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Summary  of  Recent  Results 


During  the  current  grant  period,  experimental  work  emphasized 
demonstrating  the  use  of  high  resolution  frequency  domain  nonlinear  laser 
spectroscopy  for  the  study  of  materials,  particularly  semiconductor 
heterostructures.  The  work  has  demonstrated  that  cw  four-wave  mixing 
spectroscopy  provides  new  information  regarding  the  origin  of  the  nonlinear 
optical  response  in  semiconductor  materials  as  well  as  new  information  about 
structure  and  relaxation.  We  have  also  begun  measurements  in  the  time  domain 
using  such  methods  as  three  pulse  photon  echoes,  transient  four-wave  mixing, 
and  two  wavelength  transient  absorption  in  order  to  obtain  additional 
information  as  well  as  to  enable  comparison  between  time  domain  and  frequency 
domain  measurements.  This  report  summarizes  our  experimental  and  theoretical 
results  obtained  in  GaAs/AlGaAs  quantum  well  structures  and  CdSSe  doped 
glasses. 

-  A  summary  of  the  current  progress  shows: 

First  experimental  observation  of  an  interference  effect  in  the 
nonlinear  optical  response  in  GaAs  quantum  well  structures 
corresponding  to  a  slow  component  in  the  nonlinear  optical  response. 
The  interference  effect  is  evidence  that  the  exciton  resonance 
frequency  is  shifted  in  the  presence  of  the  electron  hole  plasma 
produced  by  the  ionized  exciton. 

*  Demonstration  of  frequency  domain  nonlinear  optical  spectroscopy 
methods  for  the  direct  measurement  of  the  ambipolar  diffusion 
coefficient  and  the  electron-hole  recombination  rate  in  GaAs  quantum 
well  structures. 

°  Measurement  of  the  response  time  in  strained  quantum  well  structures 

where  the  heavy  and  light  hole  exciton  overlap  in  energy.  '<  , 

•  /  Use  of  precision  nonlinear  optical  spectroscopy  methods  for  the 

study  of  phonon  induced  spectral  diffusion  of  the  heavy  hole  exciton 
at  low  temperature  in  GaAs  quantum  well  structures.  i 

:  !  /  t 

o  Demonstration  of  a  large  cw  third  order  susceptibility  in  CdSSe 

microcrystallite  doped  glass. 

0  Demonstration  of  a  temperature  dependent  slow  component  in  the 

optical  response  of  CdSSe  microcrystallite  doped  glass  due  to  the 
presence  of  defects.  Using  precision  spectroscopy  methods,  we 
determined  the  activation  energy  for  the  trap  along  with  the 
intensity  dependence  and  the  presence  of  an  anomalous  energy 
transport  process. 

0  Demonstration  of  real  time  aberration  correction  through  cw  optical 

phase  conjugation  in  CdSSe  microcrystallite  doped  glasses. 
Demonstration  of  narrow  band,  tunable,  wide  angle  field  of  view 
filter  characteristics. 

0  Development,  of  an  analytical  description  foi  the  general 

interpretation  of  the  frequency  domain  nonlinear  optical  response  in 
terms  of  systems  which  are  quantum  mechanically  opened  or  closed. 


°  Use  of  modified  optical  Bloch  equations  to  provide  an  interpretation 

of  the  interference  profile  observed  in  the  four-wave  mixing 
spectrum  in  GaAs  quantum  well  structures. 

o  Use  of  modified  optical  Bloch  equations  to  simulate  the  nonlinear 

optical  spectrum  obtained  in  GaAs  quantum  well  structures  at  low 
temperature  showing  phonon  assisted  spectral  diffusion. 

For  a  more  complete  summary  of  the  important  results,  please  see  the  attached 
set  of  reprints  and  preprints. 


Publication  Summary 

Journal  Publications 

H.  Wang,  M.  Jiang,  and  D.G.  Steel,  "Measurement  of  phonon  assisted 
migration  of  localized  excitons  in  GaAs/AlGaAs  multiple  quantum  well 
structures,"  to  be  submitted  to  Physical  Review.  (See  attached 
preprint.  Experiments  were  done  on  this  program  but  the  data  was 
analyzed  and  the  paper  was  written  on  the  new  AFOSR  program. ) 

Hailin  Wang,  J.T.  Remillard,  M.D.  Webb,  D.G.  Steel,  J.  Pamulapati, 
J.  Oh,  and  P.K.  Bhattacharya,  "High  resolution  laser  spectroscopy  of 
relaxation  and  the  excitation  line  shape  of  excitons  in  GaAs  quantum 
well  structures,"  Surface  Science  228.  pp69-73  (1990). 

J.T.  Remillard,  H.  Wang,  M.D. Webb,  D.G.  Steel,  "Frequency  Domain 
Four-Wave  Mixing  Spectroscopy  of  Temperature  and  Optical  Intensity 
Dependent  Relaxation  in  CdSSe  Doped  Glass,"  to  be  published  in  J. 
Opt.  Soc.  Am.  B,  1990. 

J.T.  Remillard,  H.  Wang,  M.D.  Webb,  and  D.G.  Steel,  "High  Resolution 
Nonlinear  Laser  Spectroscopy  of  Room  Temperature  GaAs  Quantum  Well 
Structures:  Observation  of  Interference  Effects,"  Optics  Letters 
15,  1133  (1989). 

J.T.  Remillard,  H.  Wang,  M.D.  Webb,  D.G.  Steel,  "Optical  Phase 
Conjugation  and  Nonlinear  Optical  Band  pass  Filter  Characteristics 
in  CdSSe  Microcrystallite  Doped  Glass,"  IEEE  J.  Quant.  Elec.  QE-25 . 
p408  (1989). 

J.T.  Remillard  and  D.G.  Steel,  "Narrow  nonlinear  optical  resonances 
in  CdSSe  doped  glass,"  Optics  Letters  13,  pp30-32  (1988). 

D.G.  Steel,  J.  Remillard,  Jing  Liu,  S.C.  Rand,  "Effects  of 
population  pulsations  on  backward,  nearly  degenerate  four-wave 
mixing  spectroscopy  and  optical  phase  conjugation,"  J.  Opt.  Soc. 
Am.  B  5.  171  (1988). 

Duncan  G.  Steel  and  J.T.  Remillard,  "Resonant  nearly  degenerate 
backward  four-wave  mixing  in  open  and  closed  systems,"  Phys.  Rev.  A 
36,  pp4330-4337  (1987). 
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Hailin  Wang,  J.T.  Remillard,  M.D.  Webb,  D.G.  Steel,  "High  Resolution 
Spectroscopy  and  Photon  Echoes  of  Excitons  in  Quantum  Well 
Structures  at  Low  Temperature" ,  invited  paper  in  Laser  Spectroscopy 
IX,  1989. 

J.T.  Remillard  and  D.G.  Steel,  "Nonlinear  laser  spectroscopy  studies 
of  semiconductor  materials,"  invited  paper  in  Laser  Spectroscopy 
VIII  ,  Springer-Verlag,  W.  Persson  and  S.  Svanbe'-g,  eds.  ,  pp264-266 
(1987). 


Conference  Presentation  with  Proceedings 
Invited  Presentations 

D.G.  Steel,  J.T.  Remillard,  Hailin  Wang,  M.D.  Webb,  J.  Oh,  J. 
Pamulapati,  and  P.K.  Bhattacharya,  "High  Resolution  Nonlinear  Laser 
Spectroscopy  of  GaAs/AlGaAs  Multiple  Quantum  Wells."  OSA  Topical 
Meeting  on  Quantum  Wells  for  Optics  and  Optoelectronics,  OSA 
Technical  D’gest  10,  pp  62-65,  (1989). 

D.G.  Steel,  P.K.  Bhattacharya,  J.T.  Remillard  Hailin  Wang,  M.D. 
Webb,  J.  Pamulapati,  J.  Oh,  "High  Resolution  Nonlinear  Laser 
Spectroscopy  of  Excitons  in  GaAs  Quantum  Well  Structures,"  Quantum 
Electronics  and  Laser  Science  Conference,  OSA  Technical  Digest  12. 
MBB3  (1989). 

D.G.  Steel,  J.T.  Remillard,  H.  Wang,  and  M.D.  Webb,  "Application  of 
4-wave  mixing  spectroscopy  to  semiconductor  material  studies" , 
Conference  on  Lasers  and  Electro-Optics,  Technical  Digest  Series  7, 
pp260-261  (1988). 


Contributed  Presentations  with  Published  Proceedings 

H.  Wang,  J.T.  Remillard,  M.D.  Webb,  and  D.G.  Steel,  "High  Resolution 
Nonlinear  Laser  Spectroscopy  and  Efficient  Optical  Phase  Conjugation 
in  Semiconductor  Microcrystallite  Doped  Glasses,"  to  be  published  in 
the  Technical  Digest,  Conference  on  Lasers  and  Electrooptics,  OSA 
Technical  Digest  11.,  ppl56-157(1989) . 


H.  Wang,  J.T.  Remillard,  D.G.  Steel,  J.  Oh,  J.  Pamulapati,  and  P.K. 
Bhattacharya,  "High  Resolution  Laser  Spectroscopy  of  Relaxation  and 
the  Excitation  Line  Shape  of  Excitons  at  Low  Temperature  in  GaAs 
Quantum  Well  Structures,"  in  the  Proceedings  of  4th  International 
Conference  on  Modulated  Semiconductor  Structures,  1989.  (See 
earlier  citation  in  Surface  Science  above) 


J.T.  Remillard,  H.  Wang,  M.D.  Webb,  and  D.G.  Steel,  "Frequency 
domain  four-wave  mixing  spectroscopy  in  semiconductor  quantum  wells 
and  microcrystallites,"  International  Conference  on  Quantum 
Electronics,  Technical  Digest,  pp278-279,  IQEC  Tokyo  (1988). 
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J.T.  Remillard,  J.  Liu,  and  D.G.  Steel,  "Application  of  frequency 
domain  4-wave  mixing  spectroscopy  to  semiconductors:  studies  of 
semiconductor  doped  glasses,"  Proceedings  of  the  International 
Workshop  on  High-Speed  optical  processes  and  Optoelectronic  devices 
based  on  compound  semiconductors,  p  F9,  The  University  of  Michigan, 
1987. 


Student  Progress  Summary 

Five  students  have  been  supported  by  this  program.  Three  of  them  have 
reached  Ph.D.  candidacy  level.  One  student  received  his  Ph.D.  at  the  end  of 
winter  term,  1990.  Two  other  students  will  complete  their  Ph.D.  degree  by 
the  end  of  the  1990-1991  academic  year.  The  program  involves  students  from 
Applied  Physics,  Electrical  Engineering,  and  Physics. 


4 


PHYSICAL  REVIEW  A 


VOLUME  36.  NUMBER  9 


NOVEMBER  1,  1987 


Resonant  nearly  degenerate  backward  four-wave  mixing  in  open  and  closed  systems 

Duncan  G.  Steel  and  J.  T.  Remillard 

Departments  of  Physics  and  Electrical  Engineering,  Randall  Laboratory  of  Physics,  The  University  of  Michigan, 

Ann  Arbor,  Michigan  48109 
(Received  11  February  1987) 

This  paper  presents  perturbation  calculations  using  the  density  matrix  to  describe  the  nearly  de¬ 
generate  four-wave-mixing  spectral  response  in  resonant  systems  subject  to  state-specific  reservoir 
coupling.  The  calculations  provide  physical  insight  into  the  origin  of  narrow  resonances  observed 
in  recent  experimental  work  in  atomic  sodium  and  ultranarrow  resonances  in  crystals.  The  calcu¬ 
lations  are  for  systems  which  are  homogeneously  broadened  as  well  as  systems  which  are  inhomo- 
geneously  broadened  due  to  velocity  effects  or  random  crystal  fields.  The  calculations  show  that  it 
is  possible  to  eliminate  the  effects  of  either  kind  of  broadening  and  measure  the  decay  rate  of  the 
dipole  coherence  as  well  as  the  decay  rates  of  both  states  of  the  transition.  Furthermore,  the  cal¬ 
culations  indicate  that  the  method  of  cross-correlated  optical  fields  which  was  used  to  eliminate 
the  effects  of  laser  jitter  in  longitudinal  relaxation  rate  measurements  could  be  extended  to  the 


measurements  of  slow  dephasing  rates  in  crystals. 


INTRODUCTION 


The  study  of  line  shapes  in  nearly  degenerate  back¬ 
ward  four-wave  mixing  (NDFWM)  can  provide  impor¬ 
tant  spectroscopic  information  about  the  relaxation  of 
states  or  coherent  superposition  states  due  to  reservoir 
coupling.  A  common  example  of  this  is  spontaneous 
emission  due  to  coupling  to  the  vacuum  radiation  field. 
Spontaneous  emission  causes  decay  of  a  state  as  well  as 
decay  of  coherence  between  the  states  of  the  radiating 
transition.  There  are  also  other  important  forms  of  cou¬ 
pling  to  a  reservoir.  In  gas-phase  systems,  examples  of 
reservoir  coupling  include  state-changing,  velocity¬ 
changing,  or  phase-interrupting  collisions.  In  solids  in 
addition  to  spontaneous  emission,  examples  of  reservoir 
coupling  include  phonon  and  impurity  scattering.  De¬ 
cay  of  coherence  can  also  be  caused  by  fluctuating  fields 
in  a  crystal. 

Furthermore,  NDFWM  is  of  great  practical  impor¬ 
tance.  The  NDFWM  spectral  response  is  related  to  the 
reflectivity  bandwidth  associated  with  optical  phase  con¬ 
jugate  mirrors  produced  in  backward  four-wave  mixing. 
Such  mirrors  are  of  potential  importance  to  real-time 
aberration  correction  in  optical  systems. 

This  paper  describes  the  spectral  response  of  back¬ 
ward  NDFWM  in  a  two-level  system  in  the  presence  of 
state-specific  reservoir  coupling.  The  NDFWM  response 
is  defined  as  the  line  shape  associated  with  the  signal  as 
a  function  of  frequency  of  any  of  the  three  input  beams. 
A  number  of  line  shapes  are  discussed  for  different  cases 
of  state-specific  reservoir  coupling.  The  model  provides 
a  simple  physical  explanation  of  the  origin  of  collision- 
induced  subnatural  narrow  resonances  in  Doppler- 
broadened  gases1-2  and  the  ultranarrow  resonances  in 
solids3  observed  in  earlier  experiments.  The  physics  is 
shown  to  be  similar  to  the  origin  of  pressure-induced  ex¬ 
tra  resonances  in  FWM  (PIER4)  first  discussed  by 
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Bloembergen4-8  (see  also  review  by  Rothberg9)  and  the 
narrow  spectral  holes  and  antiholes  observed  in  modula¬ 
tion  spectroscopy  discussed  in  detail  by  Boyd. 10-14 
These  effects  arise  due  to  an  incomplete  cancellation  of 
quantum-mechanical  amplitudes.15-17  We  also  extend 
the  calculation  to  include  the  NDFWM  spectral 
response  in  a  medium  inhomogeneously  broadened  by 
random  crystal  fields.  We  show  that  in  some  cases,  it  is 
possible  to  eliminate  the  effects  of  both  inhomogeneous 
broadening  and  laser  jitter  in  the  measurement  of  the 
effective  transverse  and  longitudinal  relaxation  rate. 

MODEL 

Nearly  degenerate  four-wave  mixing  is  described  by  a 
third-order  nonlinear  optical  susceptibility  which  can  be 
evaluated  using  third-order  perturbation  theory  and 
double-sided  Feynman  diagrams.18-21  General  expres¬ 
sions  for  all  48  terms  of  the  third-order  susceptibility 
have  been  given  by  many  authors.4-22-23  In  this  calcula¬ 
tion,  we  rederive  the  third-order  nonlinear  optical  sus¬ 
ceptibility  for  a  two-level  system  for  backward  NDFWM 
in  order  to  obtain  a  physical  understanding  of  the  origin 
of  NDFWM  spectral  profiles.  The  calculation  takes  into 
account  spontaneous  emission  into  the  ground  state 
along  with  decays  of  both  states  to  the  reservoir  and 
pure  dephasing  of  the  dipole.  Because  of  its  generality 
to  multilevel  systems  and  its  extension  to  the  treatment 
of  nonexponential  decay  coupling  to  the  reservoir,  we 
use  the  density-matrix  approach  first  pioneered  by 
Bloembergen  and  Shen24  with  a  classical  description  of 
the  radiation  fields.  While  this  calculation  is  based  on 
assuming  exponential  decays  for  the  different  matrix  ele¬ 
ments,  the  density-matrix  equations  are  particularly 
powerful  since  they  can  be  generalized  to  a  more  com¬ 
plete  quantum-mechanical  transport  equation.25-28 

The  specific  form  of  the  density-matrix  equation  is 
given  as27-28 
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where  the  v-V  term  accounts  for  motion  and  v  is  the 
classical  velocity  associated  with  the  center  of  mass,  H0 
is  the  Hamiltonian  for  the  unperturbed  system,  V  is  the 
classical  electromagnetic  interaction  of  the  form  —  p-E, 
where  E  =  j'£Eae‘kaX  ‘Ua'+c.c.  is  summed  over  all 
electromagnetic  fields.  T  is  the  decay  operator;  A  is  the 
source  term  to  account  for  incoherent  pumping.  The 
other  two  terms,  (dp/dt)  |  sp  and  (dp/dt )  |  ph,  account 
for  decay  into  a  state  due  to  spontaneous  emission  from 
a  higher  state  and  decay  of  coherence  due  to  dephasing, 
respectively. 

Unlike  our  earlier  work,  we  do  not  include  any  source 
terms  from  the  reservoir.  For  example,  if  we  are 
describing  a  gas-phase  system  with  collisions,  the  reser¬ 
voir  decay  terms  account  for  loss  of  an  atom  at  velocity 
v  due  to  a  velocity-changing  collision.  However,  we  do 
not  allow  for  the  possibility  of  an  atom  in  a  given  state 
at  velocity  v'  to  undergo  a  velocity-changing  collision 
and  have  its  velocity  changed  to  v.  This  means  we  will 
fail  to  describe  such  important  behavior  as  collision- 
induced  velocity  hole  filling. 

The  polarization  is  given  by  P=Tr(/ip).  If  the  system 
is  a  gas,  then  p  is  a  function  of  velocity,  and  the  polar¬ 
ization  must  be  integrated  over  the  velocity  distribution. 
If  the  system  is  a  solid  and  inhomogeneously  broadened 
due  to  random  crystal  fields,  the  v-V  does  not  contribute 
to  the  density-matrix  equation,  but  the  polarization  must 
be  integrated  over  the  distribution  describing  the  energy 
levels.  This  integration  is  a  very  important  part  of  the 
problem  and  can  greatly  effect  the  results  and  physical 
interpretation. 

To  illustrate  the  physics  observed  in  the  experiments, 
we  consider  a  simple  two-level  system  shown  in  Fig.  1. 
The  A.,  represent  incoherent  pumping  into  level  i,  and 
the  corresponding  operator  in  the  Eq.  (1)  has  the  form 
( i  |  A  |  j  >  =  A.,  5,y .  The  T  are  exponential  decays  to  the 
reservoir  due  to  spontaneous  emission,  inelastic  col¬ 
lisions,  or  other  nonradiative-type  transitions.  The 


FIG.  1.  Energy-level  diagram  for  a  general  two-level  system. 
For  these  calculations,  the  assumed  experimental  configuration 
is  for  backward  NDFWM.  Figure  1  shows  the  appropriate  ar¬ 
rangement  of  incoming  optical  beams. 


operator  has  the  form  (i  \  T  \j  )=y,2Sy.  Note  that 
these  decays  contribute  to  the  dephasing  of  the  dipole. 
This  is  seen  by  noting  that  these  terms  contribute  to  the 
decay  of  the  off-diagonal  matrix  element,  the  optical 
coherence  pl2.  In  addition  to  this  contribution,  to  the 
decay  of  the  dipole  coherence,  there  is  additional  de¬ 
phasing.  This  term  in  the  general  equation  is  given  by 
< i  |  ( dp/dt )  1  ph  |  j )  =  —  yph(  1  —  5 jj ).  This  term  repre¬ 
sents  so-called  pure  dephasing  because  the  physical  ori¬ 
gin  for  this  term  does  not  contribute  to  any  other  de¬ 
cays.  (It  should  be  noted  that  in  the  earlier  work,1'2  ve¬ 
locity  changing  collisions  and  dephasing  collisions  arise 
from  the  same  scattering  event.  Hence,  the  effects  of 
velocity-changing  collisions  on  the  diagonal  matrix  ele¬ 
ments,  causing  loss  or  generation  of  a  state  in  a  specific 
velocity  group,  and  collision-induced  dephasing,  causing 
loss  of  dipole  coherence  and  frequency  shift  of  the  reso¬ 
nance,  is  handled  by  a  single  operator.  However,  physi¬ 
cally,  the  origin  of  the  resonant  response  does  not  de¬ 
pend  on  these  details.29)  Finally,  state  2  can  decay  by 
spontaneous  emission  back  to  1.  This  is  accounted  for 
by  a  source  term  for  state  1  given  by 


In  the  rotating-wave  approximation  (and  noting 
P'j  =Pl>  l^e  equations  for  the  matrix  elements  become 
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where  y2r=y2  +  y5p  and  rph  =  |(72r+7i)  +  7Ph-  For  level 
degeneracies,  one  must  be  more  careful  about  correctly 
including  the  effects  of  source  terms  on  coherences  be¬ 
tween  degenerate  states  and  correctly  accounting  for  the 
vector  nature  of  the  interaction  term  —  p-E. 

Figure  1  also  shows  the  optical  beam  configuration  for 
backward  NDFWM.  There  are  three  input  fields,  desig¬ 
nated  Ea,  where  a  =/,  b,  or  p  for  forward  pump,  back¬ 
ward  pump,  or  probe,  respectively.  The  nonlinear  opti¬ 
cal  interaction  of  interest  follows  from  a  third-order 
nonlinear  optical  polarization  of  the  form 
P(3l  =  7,3'£/-£(,£p*.  The  phase-matching  conditions  for 
the  wave  vectors  result  in  the  signal  field,  £, ,  counter- 
propagating  with  respect  to  the  probe  beam,  provided 
that  the  two  pump  beams  are  exactly  counterpropacat- 
ing.  By  energy  conservation,  if  the  pump  beams  are  at 
frequency  <o  and  the  probe  beam  is  at  frequencv 
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o)p  =  to  +  6,  then  the  signal  frequency  is  at  to  — 8. 

The  equation  of  motion  for  the  density  operator  is 
solved  to  third  order  in  the  electric  field.  Using  Fig.  1, 
we  can  associate  a  simple  physical  interpretation  for  the 
origin  of  the  signal  wave.  For  a  simple  system  such  as 
this  (in  the  absence  of  two-photon  resonances),  there  are 
two  important  contributions  to  the  signal.  In  the  first 
case,  the  forward  pump  and  probe  interact  to  form  an 
effective  interference  pattern  which  spatially  modulates 
the  ground-state  and  excited-state  population.  This  spa¬ 
tial  modulation  of  the  population  results  in  a  spatial 
modulation  of  the  absorption  and  dispersion.  The  “grat¬ 
ing”  spacing  is  A./[2sin(0/2)],  where  9  is  the  angle  be¬ 
tween  the  forward  pump  and  probe.  If  the  backward 
pump  is  counterpropagating  with  respect  to  the  forward 
- - 1 


pump,  then  it  Bragg  scatters  off  the  grating  in  the  direc¬ 
tion  counterpropagating  with  respect  to  the  incoming 
probe.  The  second  case  is  similar,  except  that  the  back¬ 
ward  pump  and  probe  interfere  with  a  grating  space  of 
X/[2  cost 6/2)].  In  this  case,  the  forward  pump  scatters 
from  the  spatially  modulated  population. 

DISCUSSION 

Initially,  we  consider  a  very  simple  homogeneously 
broadened  system  of  two-level  atoms  with  no  motion. 
For  a  two-level  system,  we  can  ignore  field  polarization 
effects.  The  nonlinear  polarization  for  the  case  of  fre¬ 
quency  degenerate  pump  beams  at  frequency  w,  and 
probe  beam  at  frequency  wp  =to4-6  is  given  by 


P{i)=-2Nn~XfXbX;e  iUa  iu+kp'T]\ 


8  -r'0"P' 
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where  A  is  the  pump  resonance  detuning  A=to— to0.  Xa 
is  the  Rabi  flopping  frequency  p12£a/fi  associated  with 
optical  field  Ea.  Ar0  is  the  equilibrium  population 
difference. 
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In  the  limit  of  weak  coupling  (no  pump  depletion  and  no 
coupling  of  energy  to  the  probe  beam)  the  field  equations 
decouple  and  the  frequency  dependence  (as  a  function  of 
6)  of  the  signal  intensity  is  obtained  by  looking  at  the  ab¬ 
solute  value  squared  of  the  polarization.  Equation  (5) 
shows  a  nut  her  of  denominators  with  resonances  as  a 
function  of  8.  In  general,  there  are  resonances  at  6=0 
and  ±A  with  widths  and  amplitudes  determined  by  the 
various  state  decays  and  dephasing. 

In  the  first  case,  we  examine  Eq.  (5)  in  the  limit  that 
the  system  is  absolutely  closed;  i.e.,  y,  =  y2=0  and  there 
is  no  additional  dephasing,  yph=0  and  yph  =  ysp/2.  In 
this  case  state  |  1 )  is  the  ground  state,  and  we  take 
a,  =  A2=0  and  jVo  =  <Pii  V  Figure  2(a)  shows  the 
NDFWM  spectrum  as  a  function  of  pump-probe  detun¬ 
ing,  8.  (In  this  discussion,  all  parameters  are  put  in  di¬ 
mensionless  form  by  normalizing  with  respect  to  the 
spontaneous  emission  rate.)  There  are  two  resonances. 
The  first  occurs  when  the  probe  frequency  is  coincident 
with  the  resonant  frequency  v0  at  6=—  A.  The  second 
at  6  =  A  corresponds  to  the  signal  frequency  being  reso¬ 
nant  with  to0.  The  resonances  are  associated  with  the  di¬ 
pole  coherence  and  the  width  of  the  resonances  are 
determined  by  yph.  The  resonances  are  expected  at  these 
values  of  6  because  they  correspond  to  the  frequency  of 
the  probe  where  either  the  probe  resonates  with  the  di¬ 
pole  in  absorption,  or  the  signal  resonates  with  the  di¬ 
pole  in  emission. 

The  individual  time-ordered  perturbation  sequences 


are  characterized  by  a  resonance  at  6  =  0  associated  with 
population  dynamics  (seen  in  the  second-order  contribu¬ 
tion  to  the  perturbation  sequence),  though  from  Fig.  2(a) 
there  is  obviously  no  such  resonance.  The  nonlinear 
response  in  the  vicinity  of  6  =  0  is  determined  in  part  by 
the  dynamics  of  both  the  ground  state  and  the  excited 
State.  The  dynamics  are  determined  by  the  rate  at  which 
a  specific  state  decays.  In  this  paper,  we  have  restricted 
the  discussion  to  a  simple  two-level  system  and  are  not 
considering  excited-state  or  ground-state  coherences. 
Hence,  the  tuning  of  6  can  be  interpreted  as  a  temporal 
modulation  of  the  populations.  We  first  consider  the 
origin  of  the  resonance  for  each  of  the  states  separately 


Probe  Offset  Frequency  Probe  Offset  Frequency 


FIG.  2.  NDFWM  spectrum  for  A  =  3  for  various  decays. 
All  frequencies  and  decays  are  normalized  to  the  spontaneous 
emission  rate,  (a)  The  system  is  closed  and  there  is  no  addi¬ 
tional  dephasing  (yi  =  y2  =  yPt,  =  0'.  (b>  Same  as  ta),  except  ad¬ 
ditional  dephasing  has  been  added.  yph  =  0. 5.  (cl  y ,  =  0. 1 . 
y 2  =  I ,  yPh  =  0.  (d)  yph  =  0.5,  y,=0.1.  and  y;  =  0.C5.  (See  text 
for  discussion.) 
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and  then  examine  why  in  some  cases,  the  resonance  is 
not  observed. 

The  perturbation  of  state  |  2 )  caused  by  the  optical 
field  exciting  atoms  from  state  j  1 )  decays  at  a  rate 
determined  by  the  total  decay  rate  out  of  state  |  2 ) . 
The  perturbation  of  state  j  1 )  caused  by  the  optical  field 
exciting  atoms  to  state  |  2 )  decays  with  a  rate  deter¬ 
mined  by  the  net  decay  rate  of  state  [  2 )  back  down  to 
state  [  1 )  due  to  spontaneous  emission  (a  “filling  in”  of 
state  |  1  > )  plus  the  decay  rate  out  of  state  |  1 ).  Howev¬ 
er,  for  y,  =  y2  (for  a  closed  system,  population  is  con¬ 
served  and  y1  =  y2  —  9,)  the  perturbation  of  state  |  l )  de¬ 
cays  at  exactly  the  same  rate  as  state  |  2);  hence,  there 
is  no  relaxation  of  the  ground  state  determined  by  just 
y  | :  The  resonance  width  determined  by  y,  is  multiplied 
by  an  overall  factor  of  zero  [determined  by 
1  —  ysp/(yf—  Y  !>]■  For  a  closed  system,  this  result  is 
physically  expected  since  Pn+Pn  =  ^o  for  all  time 
(The  case  for  different  state  decays  is  discussed  below  ) 

For  the  excited  state,  the  perturbation  of  state  |  2  >  is 
described  by  a  resonant  denominator  with  a  width  deter¬ 
mined  by  the  relaxation  rate  of  state  |  2  > .  The  jverall 
multiplier  of  this  factor  is  not  zero  [1  +y,p/(yf— y ,)]. 
The  resonance  occurs  at  second  order  in  the  perturba¬ 
tion  and  is  due  to  the  fact  that  at  second  order  the 
response  of  the  system  is  determined  by  the  dynamics  of 
the  diagonal  matrix  elements  or  population,  in  this  case 
state  |  2 ) .  To  understand,  however,  why  this  resonance 
is  not  observed  under  the  conditions  of  Fig.  2(a),  we  note 
that  formally,  the  numerator  obtained  by  combining  the 
two  terms  in  the  first  set  of  square  brackets  is 
—  (5  — 2/yph),  and  for  a  closed  system,  with  no  addition¬ 
al  dei  lasing,  yph  =  y,p  /2.  In  addition,  y[=y5p,  and  the 
numerator  exactly  cancels  the  resonant  denominator, 
thus  accounting  for  the  lack  of  resonance  at  6=0.  This 
represents  an  interference  of  two  different  time-ordered 
perturbation  sequences  for  the  density  matrix  operator 
such  as  the  type  represented  by  Gustafson-Yee  diagrams. 

However,  the  strength  of  this  interference  is  reduced 
in  the  presence  of  additional  dephasing,  yph,  or  nonzero 
state  decays.  Hence,  in  a  closed  system,  a  third  reso¬ 
nance  at  6  =  0  will  be  observed  in  the  presence  of  addi¬ 
tional  dephasing.  This  behavior  is  shown  in  Fig.  2(b), 
and  can  be  seen  in  the  original  equations  of  Bloember- 
gen.4  When  these  resonances  are  caused  by  the  presence 
of  pure  dephasing,  they  are  called  dephasing-induced 
coherent  emission  resonances  and  have  been  discussed 
by  Andrews  and  Hochstrasser.30  However,  the  above 
description  shows  that  these  resonances  can  also  be 
caused  by  state  decay  as  well  as  "pure  dephasing"  and  as 
the  discussion  below  will  show,  the  spectral  response  can 
be  extremely  narrow  if  yi«y2.  Finally,  for  a  closed 
system  in  the  presence  of  additional  dephasing,  the  am¬ 
plitude  of  the  6  =  0  resonance  is  determined  by  yph  and 
the  width  is  determined  by  the  spontaneous  emission 
lifetime. 

We  now  consider  the  two-level  system  with  state  cou¬ 
pling  to  the  reservoir;  i.e.,  Y\  and  Y:  arc  nonzero.  For 
small  y,  and  for  Y\  =  Yi  there  is  little  difference  in  the 
NDFWM  line  shape  with  Fig.  2(a)  except  for  some  slight 
additional  broadening  because  of  the  additional  dephas¬ 


ing  due  to  the  state  decay,  and  some  slight  filling  in 
around  6=0  due  to  the  fact  that  the  interference  be¬ 
tween  the  first-order  pump  and  probe  interactions  is  no 
longer  complete.  However,  Fig.  2(c)  shows  the 
NDFWM  response  when  yi<y2.  A  very  narrow  reso¬ 
nance  is  observed  with  a  width  determined  by  y,.  To 
understand  the  origin  of  this  behavior,  we  note  that  as 
before,  when  the  decays  are  equal,  the  resonant  enhance¬ 
ment  at  6  =  0  due  to  the  response  of  the  upper-state  dy¬ 
namics  is  reduced  by  the  interference  between  the  two 
dispersive  contributions.  In  addition,  the  ground-state 
dynamics  are  determined  by  the  decay  of  the  excited 
state  back  to  the  ground  state  in  addition  to  the  net  de¬ 
cry  out  of  the  ground  state.  Hence,  as  before,  there  is 
no  contribution  from  the  dynamics  associated  with  the 
pure  decay  of  the  ground  state  (the  factor  multiplying 
the  ground  state  6  =  0  resonance  is  still  zero).  However, 
when  y,  <y2,  the  ground-state  perturbation  does  not  de¬ 
cay  as  fast  as  the  upper  state.  Hence,  there  is  a  residual 
nonlinear  contribution  remaining  due  to  the  ground 
state.  The  dynamics  of  this  nonlinear  response  is  deter¬ 
mined  by  ylt  and  we  obtain  a  resonance  at  6  =  0  with  a 
width  given  by  y,.  In  the  case  of  an  experiment  where 
the  presence  and  magnitudes  of  the  state-specific  reser¬ 
voir  coupling  rates  is  controllable  from  zero  to  some 
finite  value,  it  would  appear  to  be  an  extra  resonance,  of 
a  type  similar  to  PIER4  resonances  referenced  earlier. 

However,  when  the  dephasing  is  large,  the  sidebands 
associated  with  resonances  at  6  =  ±A  broaden  and  disap¬ 
pear  altogether,  leaving  a  single  Lorentzian  profile  with 
a  width  given  by  y ,.  In  the  limit  of  very  large  dephasing 
<yph  »Yi'Y2-Ysp)>  this  response  is  very  similar  to  the 
work  ir.  Cr:  YA103.3  In  that  work,  the  model  was  a 
three-level  system;  however,  the  physical  result  of  decay 
to  a  metastable  state  was  to  create  a  long-lived  perturba¬ 
tion  in  the  ground  stu,*  with  the  ground-state  dynamics 
determined  by  the  decay  of  the  metastable.  We  see, 
however,  that  these  interactions  are  also  viewed  in  terms 
of  temporal  modulations  of  the  populations  at  frequency 
6  and,  hence,  are  closely  related  to  the  population  pulsa¬ 
tion  effects  which  give  rise  to  spectral  holes  in  the  non¬ 
linear  absorption  spectrum  of  homogeneously  broadened 
material. ,ou 

In  the  discussion  below,  we  will  generalize  the  calcula¬ 
tion  to  include  the  effects  of  velocity-induced  inhomo¬ 
geneous  broadening;  however,  the  physics  of  linewidth 
effects  of  reservoir  coupling  are  not  effected  by  this  addi¬ 
tion.  In  the  earlier  experiment,1,2  we  observed  collision- 
induced  line  narrowing.  For  these  experiments,  we  de¬ 
scribed  the  observations  for  A  =  0  because  of  improved 
signal-to-noise  ratio  (the  effects  were  also  observed  for 
finite  A.)  For  A  =  0,  Eq.  (5)  predicts  a  single  resonance, 
corresponding  to  the  merging  of  the  two  resonances  in 
Fig.  2(a),  with  the  width  determined  by  the  spontaneous 
emission  decay  rate  (y ,  =  y 2  =  0 ).  In  the  presence  of 
state-specific  decay,  however,  where  yi  <y;,  the  line  nar¬ 
rows  to  a  linewidth  given  by  y,  which  would  have  been 
anticipated  from  Fig.  2  in  the  limit  A— 0.  (The  physical 
origin  of  y,  for  these  experiments  is  either  velocity¬ 
changing  or  state-changing  collisions.)  We  see  from  the 
analysis,  however,  that  this  resonance  is  the  result  of  a 
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finite  contribution  from  the  g*ound-state  dynamics. 
Hence,  though  this  resonance  appears  as  a  narrowing  of 
the  spectral  response,  -the  physical  origin  is  more 
correctly  understood  as  arising  from  a  resonance  which 
is  observed  only  in  the  presence  of  state  specific  reservoir 
coupling.  This  behrv,  or  is  similar  to  the  pressure- 
induced  resonances  mentioned  earlier4-8  if  the  physical 
origin  of  y,  and  y2  is  due  to  collisions.  We  note,  howev¬ 
er,  that  there  is  really  no  fundamental  difference  between 
this  kind  of  resonance  and  the  kind  of  resonance  seen  in 
solids  such  as  Cr:YA103. 

It  should  be  noted  that  in  many  experiments,  one  de¬ 
cay  is  significant  compared  to  other  decays  resulting  in 
line  shapes  which  are  often  nearly  Lorentzian.  This  is 
not  seen  in  the  figures  shown  because  for  the  sake  of  il¬ 
lustration,  the  different  decay  rates  were  chosen  to  be 
comparable  though  diffeunt.  A  Lorentzian  line  shape  is 
unusual  for  nonlinear  spectroscopy  signatures  in  homo¬ 


geneously  broadened  systems,  but  was  seen  experimen¬ 
tally  in  the  Cr:YA103  measurements.3  This  is  in  con¬ 
trast  to  the  fully  degenerate  four-wave  mixing  response 
where  the  spectral  profile  (in  the  limit  of  weak  pumps)  is 
given  by  the  cube  of  a  Lorentzian.31 

Finally,  we  consider  the  case  where  y2<y |.  In  this 
case,  we  note  that  as  a  result  of  interference  between 
ground-state  and  excited-state  dynamics  it  is  possible  to 
produce  a  NDFWM  spectra  with  a  dip  in  the  middle  as 
hown  in  Fig.  2(d).  The  figure  shows  a  NDFWM  spec¬ 
trum  but  with  7i=0. 1,  y2  =  0.05,  and  yph  =  0.5. 

We  now  consider  the  case  where  the  atoms  are  mov¬ 
ing,  and  we  m,  t  consider  Doppler  effects  and  thermal 
washout  described  by  \V.  In  this  case,  after  velocity  in¬ 
tegration,  we  find  that  in  the  extreme  Doppler  limit32  the 
polarization  is  given  by  (assuming  a  fully  collinear 
geometry  and  a  thermal  velocity  distriburion) 
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In  the  Doppler  limit,  we  have  ignored  the  interference 
term  between  the  backward  pump  and  probe  (small  spac¬ 
ing  grating)  because  this  term  is  extremely  small  due  to 
washout  by  thermal  motion.  In  the  limit  that  5  =  0,  the 
back  grating  effects  are  discussed  very  thoroughly  by 
Lam  and  Abrams.33  Physically,  the  thermal  motion 
washout  results  when  the  transit  time  between  fringes  is 
small  compared  to  the  effective  longitudinal  relaxation 
rate.  This  motion  has  the  effect  of  lowering  the  contrast 
ratio  between  the  fringes  that  result  between  the  back¬ 
ward  pump  and  probe  that  spatially  modulate  the  popu¬ 
lation  difference,  resulting  in  a  reduced  scattering 
efficiency  for  the  forward  pump. 

In  this  case,  we  have  considered  the  problem  where 
the  reservoir  coupling  may  be  due  to  velocity-changing 
or  state-changing  collisions.  In  the  limit  that  we  do  not 
consider  velocity-changing  collisions  into  the  resonant 
velocity  group  or  consider  collision-induced  frequency 
shifts,  we  find  that  this  equation  follows  directly  from 
our  earlier  work.1  The  analytical  description  in  the  ear¬ 
lier  analysis  was  more  complete;  however,  the  physical 
implications  of  state-changing  collisions  is  not  as  easily 
seen.  To  illustrate  the  behavior,  Fig.  3(a)  shows  t>. 
NDFWM  spectral  response  as  a  function  of  pump-probe 
detuning  for  a  system  in  the  Doppler  limit  when  y,  =  y2. 
When  the  pump  beam  is  tuned  off  resonance  (with 
respect  to  the  zero  velocity  group),  we  see  that  two  reso¬ 
nances  are  observed.  The  first  resonance  occurs  at  6=0 
and  has  a  width  determined  by  the  spontaneous  emission 
rate.  The  second  resonance  occurs  at  6  =  2A  and  has  a 
width  determined  by  the  total  dephasing  rate.  In  the  ab¬ 
sence  of  state-specific  reservoir  coupling,  these  results 
are  similar  to  those  of  Yariv  and  Nilsen.  Again,  as  we 
examine  the  polarization,  we  see  that  the  resonance  at 


6=0  has  a  width  given  by  y{  but  does  not  contribute  for 
the  same  reasons  as  given  above.  However,  unlike  the 
homogeneously  broadened  system,  even  when  the  system 
is  closed,  a  resonance  occurs  at  6  =  0. 

The  physical  origin  of  the  two  resonances  is  seen  by 
consideritig  the  motion  of  the  holes  in  velocity  space  as  a 
function  of  pump-probe  detuning.  Figure  4  shows  the 
ground-state  velocity  distribution  assuming  the  pumps 
are  detuned  from  resonance  by  an  amount  A.  Assuming 
the  pumps  are  degenerate,  we  find  that  they  burn  noles 
in  the  velocity  distribution  at  velocities  given  by 
A  — k  v  =  0.  Two  holes  are  burned  corresponding  to  *  v 
locities  Doppler  shifted  into  resonance  for  the  backward 
and  forward  pump  beams.  The  probe  is  detuned  and 
burns  a  corresponding  hole,  according  to  the  amount  of 
detuning.  If  the  probe  is  given  by  w  +  6,  then  by  energy 
conservation,  the  signal  is  at  frequency  co  —  6.  A  corre¬ 
sponding  hole  is  burned  by  the  signal,  but  the  signal  k 


Probe  Offset  Frequency  Probe  Offset  Frequency 

FIG.  3.  NDFWM  spectrum  in  a  closed  system  inhomogene- 
ously  broadened  by  velocity  effects  at  A  =  2. 5.  (All  frequencies 
and  decays  normalized  to  y,p.)  (a)  Closed  system  with 

Vi  =  V:  =  yPh  =  0.  (b)  The  same  conditions  as  la)  except  now 
yph  =  0. 25.  (See  text  for  discussion.) 
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FIG.  4.  A  schematic  of  the  origin  of  NDFWM  resonances 
in  Doppler-broadened  material.  This  is  a  plot  of  the  ground- 
state  velocity  distribution  showing  four  velocity  holes  burned 
by  the  two  frequency  degenerate  pumps,  the  probe  beam,  and 
the  signal. 

vector  is  shifted  by  a  minus  sign  since  the  signal  is  prop¬ 
agating  backward  with  respect  to  the  probe  (according 
to  the  phase-matching  conditions).  Resonances  in  the 
signal  will  then  be  observed  when  holes  can  be  made  to 
overlap.  In  this  case,  the  probe  hole  can  be  adjusted. 
The  two  resonances  described  above  are  observed  to 
occur  when  6=0  and  5  =  2A.  The  first  corresponds  to 
when  the  probe  hole  aligns  with  forward  pump  hole  and 
the  second  corresponds  to  when  the  signal  hole  aligns 
with  the  forward  pump  hole.  The  remaining  resonance 
at  5=— 2A  corresponds  to  the  probe  hole  aligning  with 
backward  pump  hole,  but  is  usually  not  observed  inside 
the  Doppler  limit  because  of  thermal  washout  described 
above.  Hence,  we  see  that  the  origin  of  the  observed  res¬ 
onant  structure  can  be  interpreted  in  terms  of  velocity 
hole  burning  pictures. 

We  now  consider  the  effect  of  reservoir  coupling. 
First,  as  the  system  begins  to  experience  additional  pure 
dephasing,  we  see  in  Fig.  3(b)  that,  as  expected,  the 
second  resonance  begins  to  broaden.  To  examine  the 
case  of  state  specific  reservoir  coupling  we  consider  the 
case  for  A  =  0  (corresponding  to  experiments').  Curve 
(a)  in  Fig.  5  shows  the  NDFWM  response,  again  with  a 
width  determined  by  the  total  dephasing  for  the  case 
y,  =  y2~0  and  yph  =  0.  As  the  state  decay  rates  increase 


Probe  Offset  Frequency 

FIG.  5.  (a)  In  order  to  show  the  apparent  collisional  nar¬ 
rowing  observed  in  the  atomic  sodium  experiments  which  were 
performed  with  A  =  0,  we  show  the  predicted  NDFWM  spec¬ 
trum  under  these  conditions  with  no  state-specific  decay,  (y, 
took  on  identical  values  corresponding  to  the  transit  time 
across  the  beam.  The  value  is  0.001.)  (b)  In  the  presence  of  ve¬ 
locity  or  state-changing  collisions,  the  y,  become  different 
( y !  =  0. 1  and  y;=l)  resulting  an  apparent  narrowing  of  the 
NDFWM  spectrum.  However,  as  shown  from  the  analysis, 
this  narrowing  is  really  the  appearance  of  a  new  resonance 
with  a  width  due  to  the  ground-state  lifetime. 

in  the  presence  of  state-specific  reservoir  coupling  (due 
say  to  state-  or  velocity-changing  collisions),  we  see  in 
curve  (b)  the  line  narrows  (for  Y\<Yz^-  Nearly  identical 
behavior  is  observed  off"  resonance  (A=*=0)  in  both  the 
theory  and  experiment;  however,  in  the  experiment,  the 
signal  to  noise  becomes  poor  due  to  collision-induced  op¬ 
tical  pumping.  Agreement  between  theory  and  experi¬ 
ment  is  excellent  at  low  pressures.  At  high  pressures 
this  theory  becomes  inadequate,  and  it  is  necessary  to 
correctly  account  for  collisional  processes  that  thermal- 
ize  the  effects  of  velocity  hole  burning.117  35 

Finally,  we  consider  NDFWM  in  material  inhomo- 
geneously  broadened  by  random  crystal  fields.  In  this 
case,  we  simply  assume  a  normal  distribution  in  cu0  and 
integrate  Eq.  (5)  over  co0.  In  the  limit  of  extreme 
broadening,32  the  result  for  the  polarization  is 
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where  W  is  the  inhomogeneous  width,  and  6(/,  is  <o,  —  w;. 
In  this  polarization,  we  can  identify  two  terms  as  fp 
[ 1  /<  *>Pf  ~  Y ,  I  and  1  /( 6p/-  —  lyf']  and  two  terms 
as  bp  [  1  /< 6p/  6/fc  —  />- 1 1  and  I  /(6p/  +b,h -/>[))- 
These  terms  describe  the  dynamics  associated  with  either 
the  ground  state  (state  1)  or  the  excited  state  'state  2) 
and  are  associated  with  linewidths  y,  or  yf,  respectively. 
The  terms  identified  by  fp  are  due  to  the  interaction  by 


I 

the  forward  pump  and  probe,  and  hence  are  associated 
with  the  large  spacing  grating  (in  a  nearly  collinear  in¬ 
teraction!.  The  terms  identified  by  bp  are  due  to  the  in¬ 
teraction  by  the  backward  pump  and  probe  and  hence 
are  associated  with  the  small  spacing  grating.  In  the  line 
shape  discussions  which  follow,  we  do  not  include  con¬ 
tributions  from  the  terms  identified  by  bp  since  these 
terms  are  usually  small  even  in  solids  because  of  excita- 
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tion  diffusion.  Many  times,  these  terms  can  be  further 
reduced  experimentally,  however,  by  rotating  the  polar¬ 
ization  of  the  backwafd  pump.  The  magnitude  of  the 
backward  pump  and  probe  interaction  is  then  deter¬ 
mined  by  the  strength  of  the  alignment  term  which  is 
often  small. 

For  this  problem,  we  notice  that  there  is  no  depen¬ 
dence  on  the  absolute  resonance  frequency  (except  by 
way  of  a  Gaussian  factor  with  a  width  determined  by  W\ 
the  broadening  is  assumed  to  be  large,  and  the  value  of 
the  exponential  is  taken  as  1).  This  is  because  the 
strength  of  the  interaction  for  any  given  atom  is  the 
same  for  all  three  input  optical  fields,  independent  of  the 
field  propagation  direction.  The  only  frequency  group 
that  has  a  strong  interaction  with  the  optical  fields  is  the 
group  in  exact  resonance  with  the  fields.  Hence  there  is 
no  dependence  on  co0  except  through  the  overall  proba¬ 
bility  factor  which  determines  how  many  atoms  are  in 
the  specific  resonance  group.  These  results  are  similar  to 
those  obtained  earlier  for  polarization  spectroscopy  tech¬ 
niques  in  condensed  phases.36 

To  appreciate  the  effects  of  NDFWM  in  material  in- 
homogeneously  broadened  by  random  fields,  we  first 
note  that  in  degenerate  four-wave  mixing,  the  DFWM 
response  in  Doppler-broadened  media  is  Doppler  free. 
The  line  width  as  a  function  A  (for  8=0)  is  determined 
by  the  natural  linewidth.  Again,  the  reason  is  because  it 
is  only  for  the  optical  frequency  corresponding  to  reso¬ 
nance  for  the  zero  velocity  group  for  which  a  signal  is 
observed.  The  equal  but  opposite  Doppler  shift  for 
counter  propagating  beams  eliminates  the  signal  when 
the  laser  interacts  with  any  other  velocity  group.  This  is 
not  the  case  if  the  material  is  inhomogeneously 
broadened  by  random  fields.  The  spectral  width  of 
DFWM  would  be  determined  by  the  total  inhomogene¬ 
ous  width.  However,  if  one  of  the  optical  beams  is  de¬ 
tuned  in  frequency,  while  the  other  beams  are  held  fixed, 
it  becomes  possible  to  measure  linewidths  associated 
with  fundamental  decay  processes,  rather  than  a 
linewidth  associated  with  the  inhomogeneous  width.  In 
a  typical  system  the  probe  frequency  is  detuned  with 
respect  to  the  pumps,  then  the  linewidth  is  determined 
by  the  relaxation  rates  of  states  1  and  2  as  shown  in  Fig. 
6(a)  (6^,=0).  However,  if  the  backward  pump  beam  is 
detuned,  then  the  linewidth  is  given  by  twice  the  dipole 
coherence  relaxation  rate,  as  shown  in  Fig.  6(b)  (note  the 
scale  change  on  the  x  axis.)  This  kind  of  spectroscopy  is 
clearly  closely  related  to  hole  burning  spectroscopy.  7  It 
is  interesting  to  note  that  these  results  suggest  a  very 
powerful  aspect  to  this  kind  of  spectroscopy.  Namely, 
using  correlated  fields,  it  should  be  possible  to  eliminate 
the  effects  of  laser  jitter  on  both  transverse  and  longitu¬ 
dinal  relaxation  rate  measurements.  In  the  earlier  work, 
we  demonstrated  that  the  longitudinal  relaxation  rate 
was  measured  as  a  function  of  bfp.  Hence,  if  <oj  and  o>p 
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FIG  6.  NDFWM  spectrum  in  a  material  which  is  inhomo¬ 
geneously  broadened  by  random  crystal  fields.  In  these  curves, 
yph  was  taken  to  be  10  times  yip  and  y,  =  y,  =  0.  In  these  spec¬ 
tra,  it  is  assumed  that  the  grating  formed  by  the  backward 
pump  and  probe  is  small  either  because  of  excitation  diffusion 
or  because  the  polarization  of  the  backward  pump  is  rotated 
giving  only  a  small  alignment  contribution,  la)  The  probe 
beam  is  scanned  showing  a  width  determined  by  ylp.  (For  both 
of  these  measurements,  we  took  6/(>=0.)  (b)  The  backward 
pump  beam  is  scanned,  showing  a  width  determined  by  yph. 

are  correlated  (come  from  the  same  laser),  then  they 
have  the  same  jitter,  and  bfp  should  be  free  of  jitter,  as 
we  demonstrated.3  We  see  now  in  the  current  calcula¬ 
tion  that  the  transverse  relaxation  rate  is  measured  as  a 
function  of  8^,  which  if  we  use  correlated  fields,  would 
also  be  free  of  laser  jitter.  The  absence  of  laser  jitter  in 
these  measurements  is  because  the  nonlinear  response  is 
independent  of  the  absolute  resonance  frequency. 
Hence,  the  system  does  not  sense  changes  in  the  absolute 
laser  frequency,  and  it  should  be  possible  to  measure 
very  long  dephasing  times  without  resorting  to  the  devel¬ 
opment  of  ultrastable  lasers. 

SUMMARY 

This  paper  describes  the  physical  origin  of  many  line 
shapes  observed  in  frequency  domain  NDFWM.  While 
there  is  considerable  interest  in  using  frequency  domain 
four-wave  mixing  to  study  state  specific  reservoir  cou¬ 
pling  and  ultrafast  phenomena,38  it  is  clear  from  this 
analysis  that  it  may  become  difficult  to  interpret  any 
given  line  shape  without  having  considerable  detailed  in¬ 
formation  on  the  system  from  other  experiments. 

Note  added  in  proof.  Since  this  paper  was  first  submit¬ 
ted,  we  have  experimentally  demonstrated  the  existence 
of  narrow  spikes  and  dips  similar  to  Fig.  2  in  collision¬ 
less  systems.  These  narrow  features  are  induced  by  radi¬ 
ative  decay  and  are  present  when  either  population, 
alignment,  or  orientation  is  not  conserved.  Results  for 
systems  that  do  not  conserve  population  are  discussed  in 
Ref.  39.  Results  and  complete  theory  for  systems  that 
conserve  population  but  not  alignment  and  orientation 
will  be  discussed  elsewhere.40 
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We  describe  frequency-domain  spectroscopy  studies  of  glasses  doped  with  CdS|-,Se,  using  low-power  cw  tunable 
dye  lasers.  The  results  show  a  narrow  resonance  (4.4  kHz  at  room  temperature)  in  the  backward  nearly  degenerate 
four-wave  mixing  spectrum,  which  we  believe  is  determined  by  the  phonon-mediated  inverse  lifetime  of  a  deep  level 
trap  involved  in  the  nonlinear  response. 


Semiconductor-doped  glasses  have  been  studied  lately 
by  numerous  groups  because  of  their  large  third-order 
nonlinear-optica)  susceptibility.  Their  properties  are 
important  to  Q  switches,1  degenerate  four-wave  mix¬ 
ing  (DFWM)  and  optical  phase  conjugation,2"4  and 
optical  bistability.5  In  addition,  there  is  strong  evi¬ 
dence  that  the  semiconductor  microcrystallites  form 
by  aggregation  during  the  heating  process  of  manufac¬ 
ture  and  that  these  microcrystallites  have  a  size  distri¬ 
bution  of  the  order  of  5  nm,  which  may  make  them 
interesting  quantum -mechanical  structures.*"9 

A  number  of  laboratories  have  begun  experiments 
to  characterize  these  materials  in  terms  of  their  non¬ 
linear  response  and  the  relevant  relaxation  times.  Us¬ 
ing  pulsed  lasers,  Olbright  and  Peyghambarian10  have 
measured  n2  by  using  a  Twyman-Green  interferome¬ 
ter-  Olbright  et  al.n  followed  this  work  with  an  alter¬ 
native  approach  based  on  pulsed-laser  pump-induced 
changes  in  the  absorption  spectrum.  The  accompa¬ 
nying  theory  suggests  that  the  nonlinearity  arises  be¬ 
cause  of  band  filling  and  screening.  Similar  results 
are  obtained  by  Roussignol  et  al.12  However,  it  is 
interesting  to  note  that  many  of  the  groups  are  report¬ 
ing  different  values  for  relaxation  rates.  Jain  and 
Lind2  reported  grating  relaxation  times  in  DFWM 
shorter  than  8  nsec,  while  Yao  et  al .3  measured  grating 
decay  times  of  the  order  of  30  psec  with  a  slow  9-nsec 
tail.  In  earlier  work,  using  a  short-pulse  laser,  Rous¬ 
signol  et  al.*  measured  slow  grating  decay  times  of 
several  tens  of  nanoseconds,  and  Danielzik  et  a/.13 
measured  relaxation  times  of  the  order  of  2  nsec  in  a 
nonlinear  Fabry-Perot  configuration.  In  the  more 
recent  work,12  Roussignol  et  al.  ascribe  some  of  the 
variations  to  aging  effects  in  the  glass  filters  to  expo¬ 
sure  to  optical  radiation.  However,  it  is  interesting  to 
note  that  Cotter14  and  Nuss  et  al.1*  have  reported 
decay  times  that  are  intensity  dependent. 

Here  we  describe  nonlinear-optical  measurements 
made  based  on  four-wave  mixing  that  show  that  at 
room  temperature  at  low  intensity  the  relaxation  rate 
associated  with  the  grating  decay  time  is  of  the  order 
of  72  <2sec.  The  results  show  the  lack  of  excitation 
diffusion,  implying  the  potential  for  a  large  field  of 
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view  for  applications  to  optical  phase  conjugation, 
real-time  holography,  and  tunable  filters. 

We  used  a  standard  backward  four-wave  mixing 
configuration  for  studying  the  DFWM  and  nearly  de¬ 
generate  four-wave  mixing  (NDFWM)  response  in 
Schott  optical  filter  OG590,  which  is  doped  with 
CdSi_ISeI,  where  x  is  estimated  to  be  0.2.  In  these 
experiments,  the  probe  frequency  could  be  varied  by 
using  the  method  of  cross-correlated  optical  fields16 
for  measurements  of  bandwidths  small  compared  with 
the  laser  jitter.  A  second  tunable  dye  laser  was  used  to 
measure  the  bandwidth  associated  with  tuning  of  one 
of  the  pump  frequencies.  Phase-sensitive  detection 
was  used  to  facilitate  low-pump-power  measurements. 
Usually  the  forward  pump  beam  was  chopped. 

First  we  measured  the  optical  phase-conjugate  re¬ 
flectivity  (defined  as  the  ratio  of  the  signal  power  to 
incident  probe  power)  as  a  function  of  pump  power. 
Figure  1  shows  the  data  taken  at  \  =  589  nm.  where  the 
nonlinear  response  was  large.  For  a  weak  nonlinear 
response,  the  reflectivity  has  a  square-law  dependence 
on  the  pump  intensity.17  The  solid  line  has  an  n  =  2 
slope.  From  this  line  was  can  extract  the  nonlinear 
susceptibility  since  the  reflectivity  is  given  by  R  ~ 
[(2irw/cn)x(3)]2(8ir/cn)2L2///(,.  At  this  wavelength,  we 
infer  that  x,3)  *10""  (esu)  for  L  =  0.3  cm.  (No  local 
field  corrections  are  included.) 

A  comparison  of  the  DFWM  response  as  a  function 
of  wavelength  was  made  and  compared  with  the  ab¬ 
sorption  coefficient,  as  shown  in  the  inset  in  Fig.  1.  As 
expected,  the  nonlinear  response  decreases  at  short 
wavelength  owing  to  an  increase  in  pump  absorption. 
At  long  wavelength,  the  nonlinear  response  decreases 
owing  to  reduction  of  the  optical  coupling. 

To  measure  the  grating  relaxation  rate,  we  mea¬ 
sured  the  NDFWM  relaxation  rate  using  the  method 
of  cross-correlated  optical  fields  to  eliminate  the  ef¬ 
fects  of  laser  jitter.  We  measured  the  linewidth  at  low- 
intensity  (<  500  mW/cm2)  as  shown  in  Fig.  2.  Fitting 
a  Lorentzian  curve  to  the  data,  we  infer  a  FWHM 
linewidth  of  4400  kHz.  We  obtained  the  same 
linewidth  independently  of  the  pump  intensity  <  below 
5  W/cm2).  Furthermore,  the  linewidth  is  indepen- 
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Fig.  1.  The  FWM  reflectivity  measured  as  a  function  of 
pump  intensity.  The  straight  line  has  a  logarithmic  slope  of 
2.  The  inset  shows  DFWM  reflectivity  as  a  function  of 
wavelength  at  a  constant  pump  intensity.  The  open  points 
represent  the  absorption  spectrum  of  the  filter. 
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Fig.  2.  The  NDFWM  spectrum  observed  at  low  pump  in¬ 
tensity. 


dent  of  the  angle  between  the  pump  and  probe,  which 
shows  that  there  is  no  excitation  diffusion  on  this  time 
scale.18 

Finally,  using  a  second  dye  laser,  we  examined  the 
dependence  of  the  nonlinear  response  on  the  frequen¬ 
cy  detuning  of  the  backward  pump.  For  simple  two- 
level  resonant  systems,  this  bandwidth  of  the  detuning 
would  be  associated  with  either  the  longitudinal  or  the 
transverse  relaxation  rate.  However,  in  our  measure¬ 
ments  we  determined  that  the  nonlinear  response  was 
independent  of  the  backward  pump  detuning  as  the 
laser  was  tuned  into  the  red  (limited  by  the  laser  to  620 
nm).  As  we  tuned  the  pump  into  the  green  (X  <  590 
nm),  the  nonlinear  response  began  to  decrease,  owing 
to  absorption  of  the  backward  pump.  We  conclude 


from  this  that  it  is  unlikely  that  the  scattering  of  the 
backward  pump  beam  is  due  to  any  resonant  behavior. 

Based  on  the  above  measurements,  we  believe  that 
the  large  nonlinear  response  is  due  to  the  Filling  of 
deep  level  traps  by  the  decay  of  the  electron-hole 
excitation  produced  by  the  forward  pump  and  probe. 
Figure  3  shows  the  energy-level  diagram  believed  to 
model  the  interaction  measured  in  our  experiments. 
As  indicated  above,  for  the  case  of  a  copolarized  for¬ 
ward  pump  and  probe  and  a  cross-polarized  backward 
pump,  the  production  of  the  electron-hole  pair  is  pro¬ 
portional  to  the  product  of  the  field  amplitudes 
+  c.c.),  resulting  in  a  spatial  modulation  associated 
with  the  pair.  The  forward  pump  and  probe  are  near 
resonance  with  the  valence  band-to-conduction  band 
transition  with  energy  difference  Ec  —  Ev,  as  shown  in 
the  figure.  As  a  result  of  the  decay  to  the  traps,  the 
filled  traps  are  also  spatially  modulated. 

The  signal  results  from  a  scattering  of  the  backward 
pump  off  the  resultant  absorption  and  dispersion  grat¬ 
ing  produced  by  the  spatially  modulated  traps.  The 
backward  pump  beam  interacts  with  the  traps  through 
a  dipole  between  the  traps  and  a  state  in  the  continu¬ 
um  with  energy  Ec*.  This  interaction  is  strong  as  long 
as  the  photon  energy  exceeds  the  photoionization  en¬ 
ergy  of  the  trap.  Earlier  evidence  of  such  traps  has 
been  reported.1219 

We  take  the  nonlinear  response  to  be  due  to  trap 
Ailing  and  note  that  trap  lifetimes  are  usually  long 
compared  with  other  state  relaxation  times.  Figure  3 
shows  that  the  trap  decays  directly  to  the  ground  state, 
most  likely  by  emitting  a  photon.  If  the  trap  decay  is 
characterized  by  an  effective  net  decay  rate,  -yxp.  the 
inverse  NDFWM  bandwidth  provides  the  lifetime  of 
the  trap.  To  understand  this,  we  note  that  the  spatial 
modulation  term,  has  an  exponential  phase  de¬ 
pendence  given  by  i(w/  -  u >p)t  —  i{kf  —  kp)x.  Hence 


Fig.  3.  Energy-level  diagram  associated  with  the  NDFWM 
response.  The  decay  rate  associated  with  conduction-band 
excitation  >cv.  is  not  shown  in  this  figure. 
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the  trap  excitation  is  temporally  modulated  at  a  rate  h 
-  o)f-  up  as  the  interference  fringes  move  through  the 
medium.  When  5  <  y^p,  the  traps  can  follow  the 
modulation.  However,  when  6  >  yyp,  the  fringes  wash 
out,  lowering  the  scattering  efficiency  of  the  back 
pump. 

In  a  simple  rate-equation  approximation  using  per¬ 
turbation  theory  and  ignoring  coherences  associated 
with  band-to-band  transitions  (justified  by  the  ultra¬ 
fast  dephasing  time  of  this  transition20)  and  saturation 
effects  associated  with  trap  filling  or  trap  deexcitation 
by  photoionization,  the  nonlinear  polarization  for 
backward  NDFWM  is  given  by 

^NDFWM  ~  '  XT^tftbtp* 

X  exp[— ikp  x  -  i( u  —  6)t] 


(id  +  ycv  +  7CT)  (id  +  7TP) 

where  t*o  is  the  small-signal  absorption  coefficient  for  a 
transition  between  the  valence  band  and  the  conduc¬ 
tion  band,  xt  is  the  linear  susceptibility  associated 
with  a  filled  trap,  7cv  is  the  decay  rate  from  the  con¬ 
duction  band  to  the  valence  band,  and  7ct  is  the  decay 
rate  from  the  conduction  band  to  the  trap.  In  the 
limit  that  the  decay  rate  from  the  conduction  band  is 
large  compared  with  the  decay  from  the  trap,  we  see 
that  the  decay  rate  of  the  nonlinear  response  is  domi¬ 
nated  by  the  inverse  trap  lifetime,  7tp-  We  see  that 
the  predicted  NDFWM  spectral  response  is  Lorentz- 
ian  [the  signal  is  proportioned  to  LP(3,|2]  with  a  width 
determined  by  the  inverse  trap  lifetime.  (At  very  low 
intensity  a  NDFWM  spectrum,  such  as  that  shown  in 
Fig.  2,  shows  a  non-Lorentzian  behavior.  The  peak  is 
sharper  than  expected  for  a  Lorentzian  and  is  thought 
to  be  the  result  of  the  presence  of  a  distribution  of  trap 
lifetimes  associated  with  different  traps.  At  higher 
power,  the  curve  broadens  owing  to  photoionization  of 
the  trap,  and  the  sharp  peak  vanishes.) 

These  experiments  clearly  present  observations 
that  are  quite  distinct  from  the  earlier  results  obtained 
using  pulsed  dye  lasers.  The  relaxation  time  reported 
in  this  Letter  is  considerably  longer  than  that  reported 
earlier.  There  are  at  least  two  possible  reasons  for  this 
behavior.  The  first  is  that,  on  the  time  scale  of  the 
earlier  measurements,  the  traps  did  not  have  time  to 
fill.  That  is  to  say,  the  decay  time  from  the  conduc¬ 
tion  band  to  the  traps  is  long  compared  with  the  mea¬ 
surement  times.  The  other  possibility  is  that,  when 
high-power  pulsed  lasers  are  used,  the  traps  complete¬ 
ly  Fdl,  saturating  this  part  of  the  nonlinear-optical 
response.  Indeed,  there  is  some  evidence  for  the  lat¬ 
ter  explanation  in  the  recent  work  of  Cotter.14 

Such  long-lived  traps  suggest  that  the  trap  lifetime 
is  determined  by  a  phonon-mediated  process.  In¬ 


deed,  we  have  made  preliminary  measurements  of 
temperature  effects  and  have  found  that  as  we  cool  the 
material  the  NFDWM  linewidth  reduces  considera¬ 
bly,  implying  an  increased  trap  lifetime  and  reduced 
phonon-induced  relaxation.  Initial  measurements  of 
both  the  susceptibility  and  the  relaxation  rate  show  an 
exponential  dependence  on  the  inverse  temperature, 
as  expected  for  thermal  deexcitation  of  the  trap.  The 
activation  energy  is  estimated  to  be  at  least  93  MeV. 
However,  the  long  lifetime  suggests  that  these  traps 
may  be  associated  with  the  glass-semiconductor  inter¬ 
face. 
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We  discuss  the  effects  of  population  pulsations  on  the  spectral  response  of  backward  nearly  degenerate  four-wave 
mixing  (NDFWM).  We  include  the  effects  of  state-specific  decay  and  spontaneous  emission  on  both  two-level  and 
inverted-V-type  three-level  systems.  We  compare  the  results  with  experiments  in  both  the  gas  phase  and  the 
condensed  phase  that  show  ultranarrow  resonances.  In  addition,  we  examine  the  response  in  the  presence  of  strong 
pump  fields.  We  show  that  in  a  simple  two-level  system,  the  NDFWM  response  results  in  an  ac-Stark  splitting  of 
the  spectral  response,  whereas  in  a  simple  three-level  system  with  decay  to  a  metastable  state,  a  single  narrow 
resonance  is  observed,  which  broadens  in  the  presence  of  strong  pumps. 


1.  INTRODUCTION 

Consider  an  atom  initially  in  level  1  in  a  optical  electric  field, 
E  -  &  exp[i(k  •  x  -  ut)],  of  frequency  u>  nearly  resonant  with 
the  transition  at  frequency  o>0  between  level  1  and  an  excited 
state,  level  2.  The  probability  of  the  atom’s  being  found  in 
level  2  is  proportional  to  EE*.  If  E  is  not  monochromatic, 
but  is  instead  a  sum  of  monochromatic  fields,  E  =  £,6, 
exp[i(k,x  —  «,f)],  then  a  typical  population  term  is  propor¬ 
tional  to  the  real  part  of  &,&*  exp[i(k,  —  k,)  •  x  -  i6,;t),  where 
btJ  =  u,  —  Wj.  For  i  j,  we  then  see  in  the  rate-equation 
approximation  that  the  population  oscillates  in  time  at  fre¬ 
quency  hi,. 

Using  the  density-matrix  description,  such  population 
pulsation  effects  in  laser  systems  were  first  described  by- 
Lamb1  as  they  pertained  to  multimode  operation  in  lasers. 
The  first  observation  of  these  effects  appears  to  have  been  in 
a  saturable  absorber.2  The  first  analysis  of  these  effects 
outside  a  laser  cavity  was  by  Schwarz  and  Tan.3  These 
effects  were  later  examined  for  application  to  laser  spectros¬ 
copy  because  the  strength  of  the  pulsation  as  a  function  of  b ,, 
was  determined  by  relaxation  rates  associated  with  the 
eigenstates  of  the  transition.4  The  results  showed  that,  even 
in  homogeneously  broadened  media,  a  hole  could  be  ob¬ 
served  in  the  absorption  spectrum  of  a  probe  beam  as  a 
function  of  frequency  in  the  presence  of  a  saturating  pump 
beam. 

The  first  frequency-domain  measurement  of  the  hole  was 
made  using  modulation  spectroscopy  in  CriAljOs  (Ref.  5) 
and  was  followed  by  additional  experiments  in  alexandrite6 
and  fluorescein,7  along  with  extensive  theoretical  analysis 
pertaining  to  the  basis  of  modulation  spectroscopy.8  9  In 
these  experiments,  a  laser  beam  is  resonant  with  a  homoge¬ 
neously  broadened  system.  The  optical  beam  is  amplitude 
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modulated  at  frequency  6,  and  the  attenuation  of  the  modu¬ 
lated  signal  is  measured  using  phase-sensitive  detection. 

The  physical  origin  of  the  response  in  modulation  spec¬ 
troscopy  can  be  understood  by  considering  the  amplitude- 
modulated  beam  to  be  the  sum  of  a  central  component,  E0,  at 
frequency  w  and  two  sidebands,  E±,  at  frequency  u  ±  b. 
When  b  is  large  compared  with  the  population -difference 
relaxation  rate,  T i,  the  sidebands  interact  with  the  medium 
independently  of  the  central  component.  The  population 
terms  that  result  between  the  interference  of  the  central 
component  and  one  of  the  sidebands  oscillate  at  a  frequency 
b  that  is  large  compared  with  1/7Y  Hence  the  population 
difference  cannot  respond  to  this  interference.  However, 
for  small  6,  the  population  can  respond  to  interference  of  the 
form  EoE±*  +  E0*E±.  The  central  component  can  then 
scatter  from  this  population  difference  to  give  an  induced 
polarization  of  the  form  0]EqE±*Eo  +  dalEol2^*,  where  the 
terms  oscillate  at  frequency  u>  t  5  and  u>  ±  b,  respectively. 
Both  terms  contribute  energy  to  the  corresponding  side¬ 
band,  resulting  in  an  increased  amplitude  modulation  and 
an  apparent  reduction  in  the  attenuation  of  the  modulated 
beam  at  low  frequency.  This  effect  is  shown  dramatically  in 
Ref.  5.  From  this  discussion,  the  detected  response  in  mod¬ 
ulation  spectroscopy  can  be  seen  to  be  the  result  of  a  third - 
order  nonlinear  optical  polarization  and,  hence,  can  be  con¬ 
sidered  a  form  of  forward  four-wave  mixing. 

While  the  above  work  was  being  performed,  there  was  also 
extensive  work  in  the  area  of  resonant  backward  degenerate 
four-wave  mixing  (DFWM)  with  particular  interest  in  appli¬ 
cations  to  laser  spectroscopy  and  optical  phase  conjugation 
(OPC).10  In  OPC  by  DFWM,  a  probe  beam  interacts  in  a 
nonlinear  material  with  two  counterpropagating  pump 
beams  to  produce  a  signal  wave  that  is  phase  conjugate  to  the 
incoming  probe  wave.  For  a  weak  response,  the  signal  is  the 
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result  of  an  induced  third -order  nonlinear  optical  polariza¬ 
tion.  The  input  frequencies  of  the  three  input  beams  are  the 
same  (degenerate).  The  direction  of  propagation  of  the 
signal  is  determined  by  phase-matching  conditions  and  is 
counterpropagating  with  respect  to  the  probe  beam.  In 
resonant  DFWM,  the  forward  pump  and  probe  beam  inter¬ 
act  to  spatially  modulate  the  population  difference.  This 
spatial  modulation  or  hole  burning  produces  a  spatial  modu¬ 
lation  of  absorption  and  dispersion.  The  result  is  similar  to 
a  volume  hologram.  The  backward  pump  beam  is  counter- 
propagating  with  respect  to  the  forward  pump  and  scatters 
from  this  spatial  modulation  producing  the  signal  wave.  A 
spatial  modulation  of  the  absorption  and  dispersion  also 
results  because  of  the  interaction  between  the  backward 
pump  and  probe.  The  forward  pump  beam  scatters  from 
this  modulation,  producing  an  additional  contribution  to  the 
signal. 

To  determine  the  relative  pump-probe  frequency-detun¬ 
ing  bandwidth  for  applications  to  OPC  in  which  the  optical 
source  for  the  pump  beams  might  be  different  from  the 
optical  source  for  the  probe  beam,  the  nearly  degenerate 
four-wave  mixing  (NDFWM)  response  is  important.  In  this 
case,  as  the  frequency  of  the  probe  beam  is  detuned  from  the 
pump  beam,  the  interference  fringes  formed  by  the  forward 
pump  (backward  pump)  move  through  the  medium  in  time, 
producing  a  traveling-wave-type  excitation.  The  popula¬ 
tion  difference  becomes  temporally  modulated  (as  well  as 
spatially  modulated,  as  before).  The  temporal  modulation 
is  at  a  frequency  6  =  o>probe  -  wpunip.  In  the  limit  of  large 
dephasing  (the  total  dephasing  rate,  yphr,  is  large  compared 
with  all  state  relaxation  rates,  such  as  spontaneous  emis¬ 
sion),  effects  of  coherences  can  be  ignored,  and  the  rate- 
equation  approximation  can  be  used.  The  signal  strength  as 
a  function  of  5  is  a  maximum  at  5  =  0  and  falls  off  (as  a 
Lorentzian)  as  6  becomes  comparable  with  the  population- 
difference-relaxation  rate.  Hence,  in  contrast  to  the  work  of 
modulation  spectroscopy  in  which  a  hole  is  observed  in  a 
probe  absorption  spectrum,  NDFWM  produces  a  resonant 
emission  in  the  response.  In  both  cases,  the  width  of  the 
spectral  feature  can  be  narrower  than  the  natural  linewidth 
determined  by  dephasing.  In  some  cases,  the  width  of  the 
spectral  response  can  be  narrower  than  the  width  associated 
with  the  spontaneous  emission  rate. 

In  Section  2  we  will  see  that  this  simple  explanation  be¬ 
comes  much  more  complicated  in  the  presence  of  state- 
specific  reservoir  coupling.  Furthermore,  for  a  closed  or 
nearly  closed  quantum  system,  the  effects  of  coherence  be¬ 
come  important  and  can  result  in  interference  effects  that 
eliminate  the  observation  of  the  central  pump-probe  detun¬ 
ing  resonance  at  5  =  0  associated  with  population  pulsations. 

Perturbation  theory  is  adequate  to  describe  the  four-wave 
mixing  interaction  when  the  fields  are  weak.  However,  for 
efficient  performance  of  a  phase-conjugate  mirror,  pump 
intensities  well  in  excess  of  the  saturation  intensity  are  re¬ 
quired."  For  a  simple  two-level  system,  this  effect  corre¬ 
sponds  to  having  a  Rabi  flopping  frequency,  x  =  nE/h,  much 
larger  than  the  natural  linewidth  determined  by  the  total 
dephasing  rate,  ->phr  A  large  Rabi  flopping  frequency  re¬ 
sults  in  an  apparent  splitting  of  the  atomic  transition,  called 
the  ac-Stark  effect.  This  effect  is  accompanied  by  a  three- 
peaked  fluorescence  spectrum  in  resonant  scattering. '- 
This  effect  also  leads  to  a  splitting  of  the  NDFWM  response 


into  a  central  peak  along  with  two  symmetric  sidebands 
shifted  by  the  Rabi  flopping  frequency13-16  and  necessitates 
a  theoretical  approach  that  accounts  for  the  effect  of  the 
strong  pumps  on  the  energy  levels.  However,  when  optical 
pumping  effects  are  present,  as  they  can  be  in  multilevel 
systems,  four-wave  mixing  sidebands  disappear  and  are  re¬ 
placed  by  a  simple  broadening  of  the  NDFWM  response  that 
is  linear  in  the  pump  intensity,  as  shown  in  Section  2. 

2.  DISCUSSION 

A.  Stationary  Two-Level  Atoms  below  Saturation 
In  the  limit  of  low  intensities,  the  third-order  nonlinear 
optical  susceptibility  describing  NDFWM  can  be  evaluated 
using  third-order  perturbation  theory.  The  experimental 
geometry  is  schematically  shown  in  Fig.  1(a)  (see  Refs.  15 
and  16  for  a  detailed  discussion  of  the  experimental  setup) 
and  consists  of  two  counterpropagating  pump  beams  desig¬ 
nated  E f(uf,  kf)  and  Efc(<»b,  kfc)  and  a  probe  beam  Ep(u;p,  kp). 
For  optical  phase-conjugation  experiments,  we  usually  take 
aif  =  uh  —  ai,  k /  =  —k(„  and  up  =  u  +  5.  The  signal  wave, 
E,(w„  k,),  is  then  usually  counterpropagating  with  respect  to 
the  probe  beam  with  w,  =  w  —  5.  In  this  section,  we  examine 
the  third-order  nonlinear  optical  susceptibility  for  a  two- 
level  system  for  backward  NDFWM  to  understand  the 
NDFWM  spectral  profiles.  The  calculation  considers  spon¬ 
taneous  emission  into  the  ground  state  along  with  decays  of 
both  states  to  the  reservoir  and  pure  dephasing  of  the  dipole. 
Figure  1(b)  shows  the  energy-level  diagram.  These  calcula¬ 
tions  will  deviate  from  expectations  based  on  the  traditional 
optical  Bloch-vector  formulation  because  we  account  for  the 
specific  decay  of  each  state,  which  can  lead  to  the  noncancel¬ 
lation  of  quantum-mechanical  amplitudes.  The  results  will 
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Fig.  1.  (a)  Geometry  for  backward  four-wave  mixing,  ibl  Two- 

level  model  with  spontaneous  emission  and  decay  to  the  reservoir, 
icl  Closed  three-level  model.  The  input  optical  radiation  interacts 
only  with  the  transition  between  level  1  and  level  Level  t  is 
assumed  to  be  metastable  with  respect  to  transitions  to  the  ground 
state. 
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also  deviate  from  a  rate-equation  approach  since  we  allow 
for  the  effects  of  interference  between  coherences  induced 
by  the  different  beams.  We  then  extend  the  analysis  to  the 
simple  three-level  system  shown  in  Fig.  1(c)  and  show  how 
the  theory  for  a  two-level  system  anticipates  the  physics  in 
the  three-level  system. 

In  this  work,  we  perform  a  semiclassical  calculation  using 
classical  fields  and  the  equation  of  motion  for  the  quantum- 
mechanical  density-matrix  operator  presented  by  Bloem- 
bergen  and  Shen.17  This  calculation  is  based  on  a  simple 
model  for  reservoir  coupling,  but  the  density-matrix  de¬ 
scription  is  convenient  because  a  much  more  complete  the¬ 
ory  for  reservoir  coupling  can  be  included.18 

The  specific  form  of  the  density-matrix  equation  is  given 
as1920 

ih  (|  +  v.v)p  =  [H0,  p]  +  [  V,  p]  -  |  [r.  p)  + 

+  ih^-  - ih ^  +  ihA,  (1) 

Cu  sp  af  ph 

where  the  v  •  V  term  accounts  for  motion  and  v  is  the  classi¬ 
cal  velocity  associated  with  the  center  of  mass,  Ho  is  the 
Hamiltonian  for  the  unperturbed  system,  V  is  the  classical 
electromagnetic  interaction  of  the  form  —  p  •  E,  where  E  «= 
(l/2)£  Ea  exp(ika  •  *  —  iwa‘)  +  c.c.  is  summed  over  all 
electromagnetic  fields.  T  is  the  decay  operator  that  repre¬ 
sents  decay  of  any  state  to  the  reservoir.  Terms  of  the  form 
~(ift/2)[r,  p)+  represent  decay  of  a  specific  state  due,  for 
example,  to  spontaneous  emission.  Another  example  is  the 
effect  of  a  velocity-changing  collision  in  a  Doppler-broad¬ 
ened  system  where  the  atom  is  Doppler  shifted  out  of  reso¬ 
nance  after  the  collision.  A  is  the  source  term  to  account  for 
incoherent  pumping.  The  terms  dp/dtl,p  and  dp/dt|Ph  ac¬ 
count  for  decay  into  a  state  because  of  spontaneous  emission 
from  a  higher  state  and  decay  of  coherence  due  to  dephasing, 
respectively. 

Unlike  in  our  earlier  work,2122  we  do  not  include  any 
source  terms  from  the  reservoir,  such  as  velocity-changing 
collisions  that  shift  a  velocity  group  into  resonance.  Thus 
we  will  fail  to  describe  such  important  behavior  as  collision- 
induced  velocity  hole  filling  and  Dicke  narrowing.  Howev¬ 
er.  at  low  buffer-gas  pressures,  the  present  description  is 
quite  adequate  to  permit  an  understanding  of  the  experi¬ 
mental  observations. 

The  polarization  is  given  by  P*  Tr(*ip).  If  the  system  is  a 
gas.  then  p  is  a  function  of  velocity,  and  the  polarization 
must  be  integrated  over  the  velocity  distribution.  If  the 
system  is  a  solid  and  is  inhomogeneously  broadened  due  to 
random  crystal  fields,  the  v  •  V  does  not  contribute  to  the 
density-matrix  equation,  but  the  polarization  must  be  inte¬ 
grated  over  the  frequency  distribution  of  energy  levels. 

To  illustrate  the  physics  observed  in  the  experiments,  we 
consider  the  two-level  system  shown  in  Fig.  1(b).  The  A, 
represent  incoherent  pumping  into  level  i,  and  the  corre¬ 
sponding  operator  in  Eq.  (1)  has  the  form  (ilAlj)  ■ 

The  T  are  exponential  decays  to  the  reservoir  due  to  sponta¬ 
neous  emission,  inelastic  collisions,  or  other  types  of  transi¬ 
tions.  The  operator  has  the  form  <ilr[/)  -  y„8,,.  Note  that 
these  decays  contribute  to  the  decay  of  the  induced  optical 
coherence,  pi2.  In  addition,  some  dephasing  caused  by  in¬ 
teractions  with  the  reservoir  does  not  cause  transitions  and 


is  included  in  the  general  equation  by  setting  (lldp/dtlph  l'2) 
=  "tphPii-  This  term  represents  so-called  pure  dephasing 
because  the  physical  origin  for  this  term  does  not  contribute 
to  any  other  decays.  Finally,  state  2  can  decay  by  spontane¬ 
ous  emission  back  to  state  1.  This  decay  is  accounted  for  by 
a  source  term  for  state  1  given  by  <  lldp/<3tlspll  >  =  ysppj_>.  In 
the  rotating-wave  approximation  (and  noting  that  p ,,  =  p,,*), 
the  equations  for  the  matrix  elements  become 


^  +  V  ‘V)Pl1  =  (  V^P2\  ~  C-C.)  ”  ‘ftflPll  +  l/l75pP22  +  lhxv 

(2) 


ih  +  v  •  —  hu>gP|2  +  ( U12P00  p 1 1 | o )  ih  Tph^Pi^, 

(4) 


where  y2r  =  72  +  7sPand  >phr  =  l/2(>3r+  7i>  +  7Ph-  These 
equations  are  easily  extended  to  multilevel-type  systems, 
and  multiple-photon  effects  (e.g.,  two-photon  resonances) 
and  level  degeneracies  are  easily  incorporated.  For  level 
degeneracies,  the  effects  of  source  terms  on  coherences  be¬ 
tween  degenerate  states  must  be  included,  and  the  vector 
nature  of  the  interaction  term,  —  p  •  E,  must  be  accounted  for 
correctly. 

The  four-wave  mixing  interaction  of  interest  here  arises 
from  a  third-order  nonlinear  optical  polarization  of  the  form 
pot  -  ^EfEbEp*.  The  phase-matching  condition  for  the 
wave  vectors  results  in  the  signal  field,  E„  propagating 
counter  to  the  probe  beam,  provided  that  the  two  pump 
beams  are  exactly  counterpropagating  and  have  the  same 
frequency.  If  the  pump  beams  re  at  frequency  w  and  the 
probe  beam  is  at  frequency  uip  m  w  +  8,  then  the  signal 
frequency  is  at  w  —  8.  The  equation  of  motion  for  the 
density  operator  is  solved  to  third  order  in  the  electric  field. 
Furthermore,  in  the  limit  that  the  interaction  is  weak  (I EJ 
Ej2  «  I,  i  *  f,  b,  p ),  then  the  signal  intensity  is  simply 
proportional  to  [P31]2. 

To  understand  the  physics,  we  first  consider  a  simple 
homogeneously  broadened  system  of  two- level  atoms  with 
no  motion.  For  a  two-level  system,  we  can  ignore  field- 
polarization  effects.  The  nonlinear  polarization  for  the  case 
of  frequency-degenerate  pump  beams  at  frequency  u>  and 
probe  beam  at  frequency  o>p  is  given  by 

p3>  »  —2 N0  ~  XfXbXp *  exp[j(k,  +  k6  -  kp)  •  x  -  i(w  -  6)f  ] 

X  |(A  +  i7phrrl  +  [(-A  -  6)  +  »7phTll 

X  [(-A  +  8)  -  iy,*]-'  ^1  -  ^  ''I  y  )  <6  ~  ,>l)" 

+  ^1  + — jr*- — ^(5  -  ty2T)_lj  +  c.c.,  (5) 

where  A  is  the  pump  resonance  detuning  A  =  m  -  u.»,  is  the 
Rabi  flopping  frequency  (m ?EJh)  associated  with  optical 
field  and  No  is  the  equilibrium  population  difference: 
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PROBE  OFFSET  FREQUENCY 

Fig.  2.  (a)  NDFWM  spectrum  for  a  completely  closed  system  of 
stationary  atoms  with  no  additional  dephasing.  Interference  ef¬ 
fects  eliminate  the  resonance  at  6  *  0.  (b)  In  the  presence  of 
additional  dephasing,  the  interference  is  no  longer  complete,  and  a 
resonance  at  6  *  0  is  observed,  (c)  Both  states  decay  to  the  reser¬ 
voir.  but  >i  <  r>.  The  narrow  resonance  has  a  width  given  by  7i- 
(All  frequencies  are  normalized  to  the  spontaneous  emission  rate.) 


Equation  (5)  shows  a  number  of  denominators  with  reso¬ 
nances  as  a  function  of  J.  In  general  there  are  resonances  at 
6*0  and  ±A  with  widths  and  amplitudes  determined  by  the 
various  state  decays  and  dephasing. 

In  the  limit  that  the  system  is  absolutely  closed  (i.e..  it  * 
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>2  =  0  and  the  total  population  is  conserved)  and  there  is  no 
additional  dephasing  (7Ph  =  0  and  =  7 .p/2),  then  state 
ll)  is  the  ground  state,  and  we  take  Xi  =  X2  =  0  and  No  = 
(pn)eq-  Figure  2(a)  shows  the  NDFWM  spectrum  as  a  func¬ 
tion  of  pump-probe  detuning,  6.  (In  this  discussion,  all 
parameters  are  put  in  dimensionless  form  by  normalizing 
with  respect  to  the  spontaneous  emission  rate.)  There  are 
two  resonances.  The  first  resonance  corresponds  to  the  case 
when  the  probe  frequency  is  coincident  with  the  resonant 
frequency  wo  at  6=  -A  and  represents  the  strong  interaction 
that  occurs  when  the  input  probe  beam  resonates  with  the 
dipole.  The  second  resonance  occurs  when  the  signal  fre¬ 
quency  is  coincident  with  the  resonant  frequency  at  b  =  A. 
This  second  resonance  is  expected  because  it  is  the  result  of  a 
dipole  array’s  reradiating  a  wave  at  the  natural  frequency, 
wo.  The  resonances  are  associated  with  the  dipole  coher¬ 
ence,  and  the  width  of  the  resonances  is  determined  by  7phr. 

The  individual  time-ordered  perturbation  sequences  are 
characterized  by  a  resonance  at  6  =  0  associated  with  popula¬ 
tion  dynamics  (seen  in  the  second-order  contribution  to  the 
perturbation  sequence),  although  Fig.  2(a)  shows  no  such 
resonance.  (This  resonance  does  appear  in  a  rate-equation 
approximation.)  However,  the  resonance  is  not  observed 
because  of  interference  effects  between  two  different  coher¬ 
ent  excitations  of  the  populations.  These  interference  ef¬ 
fects  completely  eliminate  contributions  to  the  signal  from 
population  modulations.  More  precisely,  this  resonance  is 
not  observed  because  of  an  interference  of  two  different 
time-ordered  perturbation  sequences  for  the  density-matrix 
operator  such  as  the  type  represented  by  Gustafson-Yee 
diagrams.23  However,  in  the  case  when  additional  dephas¬ 
ing  is  present,  the  interference  is  no  longer  complete,  and  a 
central  resonance,  due  to  population  modulation,  is  ob¬ 
served  [see  Fig.  2(b)].  The  width  of  this  resonance  is  deter¬ 
mined  by  the  relaxation  rates  associated  with  the  popula¬ 
tions.  We  note  that  in  the  past  this  resonance  has  been 
called  a  dephasing-induced  coherence.24  The  basic  behav¬ 
ior  is  included  in  the  original  work  of  Bloembergen  et  al.'a 
However,  the  calculations  show  that  this  resonance  would 
also  be  observed  if  and  y2  were  nonzero,  even  if  yph  -  0. 

To  consider  the  origin  of  the  width  of  this  resonance  in 
more  detail,  the  nonlinear  response  in  the  vicinity  of  6  =  0  is 
determined  by  the  dynamics  of  both  the  ground  state  and 
the  excited  state.  The  dynamics  are  determined  by  the  rate 
at  which  a  specific  state  decays.  The  perturbation  of  state 
|2>  decays  at  a  rate  determined  by  the  total  decay  rate  out  of 
state  l2>.  However,  for  71  =  y>  the  perturbation  of  state  ll  > 
decays  with  a  rate  determined  by  the  net  decay  rate  of  state 
|2)  back  down  to  state  ll)  because  of  spontaneous  emission 
(a  filling  in  of  state  1 1 ) )  plus  the  decay  rate  out  of  state  1 1 ) . 
the  perturbation  of  state  ll )  decays  at  exactly  the  same  rate 
as  state  |2>;  hence  there  is  no  relaxation  of  the  ground  state 
determined  by  just  >i,  and  the  overall  linewidth  is  deter¬ 
mined  by  y>r  *  y,p  +  y>.  The  resonance  whose  width  is 
determined  by  >i  is  multiplied  by  an  overall  factor  of  zero 
[determined  by  1  -  [>,(>/(  v>r  -  >  i )])• 

However,  if  there  is  state-specific  reservoir  coupling  and 
7 1  and  7_>  are  not  equal  and  not  zero,  then  the  spectral 
response  may  be  greatly  modified.  In  particular,  if  •> ,  <  7  _. 
then  the  decay  of  the  ground-state  perturbation  by  a  filling 
in  from  the  excited  state  is  incomplete  and  is  not  offset  by 
the  ground-state  decay  to  the  reservoir,  as  it  was  earlier  In 
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Fig.  3.  NDFWM  response  in  Cr:YA103.  The  linewidth  is  9.8  Hz. 
The  inset  is  the  fluorescent  lecay  from  the  metastable  state.  The 
decay  rate  is  34  msec,  in  very  good  agreement  with  the  ground-state 
lifetime  inferred  from  the  NDFWM  linewidth. 


this  case,  we  expect  to  see  a  residual  perturbation  from  the 
ground  state  with  a  width  determined  by  7i-  This  result  is 
shown  in  Fig.  2(c).  When  the  presence  and  magnitudes  of 
the  state-specific  reservoir  coupling  rates  are  able  to  be  con¬ 
trolled  experimentally  from  zero  to  some  finite  value,  it 
could  appear  to  be  an  extra  resonance  of  a  type  similar  to 
PIER4  resonances-6-29  in  forward  four-wave  mixing  or  a 
narrowing  of  the  response  obtained  in  NDFWM  when  >i  = 
70. 21,22  The  calculations  show  that  in  fact  the  narrow  re¬ 
sponse  is  the  result  of  a  finite  contribution  of  an  existing 
resonance. 

When  the  dephasing  is  large,  the  sidebands  associated 
with  resonances  at  6  =  ±A  broaden  and  disappear  altogeth¬ 
er,  leaving  a  single  Lorentzian  profile  with  a  width  given  by 
7 1 .  We  see,  however,  that  these  interactions  are  also  viewed 
in  terms  of  temporal  modulations  of  the  populations  at  fre¬ 
quency  6  and  hence  are  closely  related  to  the  population 
pulsation  effects  that  give  rise  to  spectral  holes  in  the  nonlin¬ 
ear  absorption  spectrum  of  homogeneously  broadened  mate¬ 
rial  just  discussed. 

B.  Stationary  Three-Level  Atoms  below  Saturation 
Now  consider  the  NDFWM  response  of  the  three-level  sys¬ 
tem  in  Fig.  1(c).  The  fields  are  resonant  with  the  transition 
between  the  ground  state  and  level  2.  This  system  is  closed, 
although  the  upper  state  of  the  optical  excitation  can  decay 
to  level  3,  which  may  be  metastable  with  respect  to  transi¬ 
tions  to  the  ground  state.  If  we  consider  the  decay,  due  to 
optical  excitation,  of  the  excited-state  perturbation  (level  2), 
we  note  that  it  decays  at  the  rate  7sp  +  723.  However,  the 
ground-state  perturbation  decays  by  a  filling  in  by  spontane¬ 
ous  emission  from  level  2  and  a  long-lived  decay  from  level  3. 
If  we  take  731  to  be  small  compared  with  other  decays,  we 
expect  to  observe  a  NDFWM  spectral  response  whose  width 
is  determined  by  the  ground-state  lifetime,  731  (also  the 
lifetime  of  the  metastable). 

This  response  is  calculated  by  slightly  modifying  Eq.  (1) 
and  accounting  for  the  equation  of  motion  for  P33.  Equation 
(2)  is  modified  by  adding  a  term  +731P33  to  the  left-hand 
side.  In  Eq.  (3),  y2r  ■  7,p  +  733-  For  a  solid,  we  ignore  the 
velocity  effects.  For  P33,  we  have  P33  ”  luPn  ~  7uP33- 
Because  the  system  is  closed,  we  take  X,  ■  0  and  note  that 


probability  is  conserved  in  this  system:  p\  1  +  p> ■<  +  p.u  =  -V,,. 
For  this  system,  the  polarization  is  given  by16 

p(3)  _  -2 N0  ~  XfXbXp *  exp[i(k,  +  k6  -  kp)  •  x  -  Uw  -  6)t] 

x  r  1  +  1  (  ^  .  t,p-  Y31V 

|_6  —  i72T  >2r  “  Y31  y  ~  *Tsi 

[(-A  +  6)  -  iyphT  (a  +  iyphT  (A  +  b)  -  i7phr)_ 


+  c.c. 


(7) 


We  see  in  the  limit  that  731  is  much  less  than  the  other 
decays  and  that  the  response  is  dominated  by  the  resonant 
denominator  5  +  iy as  we  anticipated.  Using  the  method 
of  cross-correlatH  fields  to  eliminate  the  effects  of  laser 
jitter,  we  measured  the  NDFWM  response  in  Cr:YA10|.16 
Figure  3  shows  the  measurement  with  a  width  of  9.8  Hz  (7  u/ 
it).  The  solid  line  is  a  least-squares  fit  of  a  Lorentzian  curve 
showing  excellent  agreement  with  the  analysis.  The  inset 
shows  the  measured  decay  of  the  fluorescence  from  the 
metastable  state.  The  decay  is  exponential  with  a  decay 
time  (731-1)  of  34  msec,  in  good  agreement  with  the 
NDFWM  measurement. 


C.  Doppler-Broadened  Two-Level  Atoms  below 
Saturation 

We  now  consider  the  case  in  which  the  atoms  are  moving  and 
we  must  consider  Doppler  effects  and  thermal  washout  de¬ 
scribed  by  v  •  V.  In  this  case,  after  velocity  integration,  we 
find  that  in  the  extreme  Doppler  limit30  the  polarization  is 
given  by  (assuming  a  fully  collinear  geometry  and  a  thermal- 
velocity  distribution) 

f*3)  *  Mi zX/XfcXp*  expj— i((uj  -  5 )t  +  kp  •  x]| 


2  ir 


1/2 


5  —  2A  —  2i(712  +  79p/2  +  7ph> 

:[(»  +  -- - —) 

l\  >1  -  72  -  7,p/  5  -  17,  \  >1  -  >2  -  ><„/ 


1 

“  ‘(72  +  7sp)_ 


+  c.c. 


(8) 


In  the  Doppler  limit,  we  have  ignored  the  interference  term 
between  the  backward  pump  and  probe  (small  spacing  grat¬ 
ing)  because  this  term  is  extremely  small  owing  to  washout 
by  thermal  motion.  In  the  limit  that  6  ■  0,  the  backgrating 
effects  are  discussed  thoroughly  by  Lam  and  Abrams. 11 

In  this  case,  we  have  considered  the  problem  in  which  the 
reservoir  coupling  may  be  due  to  velocity-  or  state-changing 
collisions  or  to  spontaneous  emission  to  another  state.  In 
the  limit  that  we  do  not  consider  velocity-changing  collisions 
into  the  resonant  velocity  group  or  consider  collision-in¬ 
duced  frequency  shifts,  we  find  that  this  equation  is  similar 
to  our  earlier  work  on  collision  effects.21  The  analytical 
description  in  the  earlier  analysis  was  more  nearly  complete. 
However,  we  have  tried  to  emphasize  a  more  general  physi¬ 
cal  picture  in  this  presentation. 

To  illustrate  the  behavior  in  the  presence  of  atomic  mo- 
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Fig.  4.  (a)  The  NDFWM  spectrum  in  Doppler-broadened  materi¬ 
al:  A  ■  2.5-y,p.  A  fully  collinear  geometry  is  assumed,  (b)  The 
NDFWM  spectrum  measured  in  atomic  sodium  on  the  F  «  2  —  F  = 
3  transition  of  the  D2  line  at  5890  A.  (c)  The  NDFWM  spectrum  in 
the  presence  of  state-specific  reservoir  coupling  with  >i  <  >2-  (d) 
The  NDFWM  spectrum  measured  on  the  same  transition  as  (b)  but 
with  25  Torr  of  neon  added.  [The  frequencies  in  curves  (a)  and  (c) 
are  normalized  to  the  spontaneous  emission  rate.) 


tion,  Fig.  4(a)  shows  the  NDFWM  spectral  response  as  a 
function  of  pump-probe  detuning  for  a  system  in  the  Dopp¬ 
ler  limit  when  =  y-j  and  a  fully  collinear  geometry  (i.e.,  k/ 
=  kp  =  -kfc).  When  the  pump  beam  is  tuned  off  resonance 
(with  respect  to  the  zero-velocity  group),  we  see  that  two 
resonances  are  observed.  The  first  resonance  occurs  at  6  *  0 
and  has  a  width  determined  by  the  spontaneous  emission 
rate.  The  second  resonance  occurs  at  6  =  2A  and  has  a  width 
determined  by  the  total  dephasing  rate.  Figure  4(b)  is  a 
comparison  with  experiments  in  atomic  sodium  on  the  D 2 
line  (-S1/2(F  =  2)  —  -P.v^F  =  3)  transition).  (In  these 
measurements,  the  angle  between  the  pump  and  probe  re¬ 
sulted  in  a  residual  Doppler  width  of  the  order  of  a  few 
megahertz,  which  eliminated  the  observation  of  the  effects 
of  interference  due  to  coherences  between  the  magnetic  sub¬ 
states. 1J)  In  the  absence  of  state-specific  reservoir  coupling, 
the  analytical  results  are  similar  to  those  of  Nilsen  and 
Yariv. 1  ’  Again,  as  we  examine  the  polarization,  we  see  that 
the  resonance  at  6  «  0  that  has  a  width  given  by  >i  does  not 
contribute  for  the  same  reasons  as  just  given.  However, 
unlike  with  the  homogeneously  broadened  system,  even 
when  the  system  is  closed  a  resonance  occurs  at  6  =  0. 

We  now  consider  the  effect  of  reservoir  coupling.  As  the 
system  begins  to  experience  additional  pure  dephasing,  the 
second  resonance  begins  to  broaden  and  decrease  in  ampli¬ 
tude.  In  the  presence  of  state-specific  reservoir  coupling 
where  y 1  <  y>,  we  expect  the  6*0  resonance  to  narrow. 
Because  of  the  additional  dephasing,  the  second  resonance 
will  broaden  and  weaken.  Figure  4(c)  shows  the  calculated 
NDFWM  response  for  A  =  0  (A  =  0  was  used  in  the  experi¬ 
ments  to  give  the  best  signal-to-noise  ratio).  To  make  the 
comparison  experimentally,  ground-state  perturbers  (such 
as  neon  buffer  gas)  were  added  to  the  cell  of  atomic  sodium. 
In  this  case,  cross  sections  for  state-  or  velocity-changing 
collisions  between  the  sodium  atom  and  the  perturber  atom 
are  larger  for  the  3p  state  than  for  the  3,<  state.  Hence,  in 


these  experiments,  the  addition  of  buffer  gas  resulted  in  y  i  < 
7->.  Figure  4(d)  shows  the  measured  narrowing  of  the 
NDFWM  response.34 

On  the  basis  of  the  work  just  discussed,  such  resonances 
should  be  observed  if  the  upper  state  can  decay  to  another 
state  by  spontaneous  emission.  One  example  is  the  2Si/2(F 
=  2)  -*  2P|/2(F  =  2)  of  the  D1  line  in  sodium.35  Dipole- 
selection  rules  allow  the  upper  state  to  decay  to  the  F  *  1 
ground  state.  In  steady  state,  this  would  lead  to  optical 
pumping  of  the  hyperfine  split  ground  state.  (This  is  not 
the  case  on  the  previous  F  =  2  -*■  F  *  3  transition  on  the  D 2 
line).  In  this  case,  even  in  the  absence  of  any  external 
reservoir  coupling,  the  coupling  to  the  vacuum  radiation 
field,  giving  rise  to  spontaneous  emission,  is  sufficient  to 
make  the  interference  effects  incomplete,  giving  rise  to  a 
narrow  NDFWM  response  on  the  6  =  0  resonance.  Figure 
5(a)  shows  the  calculated  NDFWM  response  on  this  transi¬ 
tion  if  decay  to  the  F  =  1  ground  state  is  not  allowed.  If  the 
decay  is  allowed  (using  the  appropriate  decay  rates  for  this 
transition),  Fig.  5(b)  shows  the  expected  narrow  response; 
with  the  width  being  given  by  the  transit  time.  To  compare 
these  results  with  experiment,  it  is  necessary  to  account  for 
the  finite  angle  between  the  forward  pump  and  probe,  which 


Fig.  5.  The  calculated  NDFWM  response  in  atomic  sodium  for  the 
F  ■  2  —  F  *  2  transition  on  the  Dl  line.  For  this  calculation, 
spontaneous  emission  decay  to  the  F  ■  l  ground  state  was  not 
allowed.  (The  frequency  is  normalized  to  the  spontaneous  emission 
rate  !  (b)  The  same  calculation  as  in  (a)  except  the  F  *  upper 
level  is  allowed  to  decay  to  the  F  ■  1  ground  state.  The  solid  line  is 
for  a  fully  collinear  geometry.  The  dashed  line  is  the  calculation 
allowing  for  the  residual  Doppler  broadening  of  the  A  *  0  resonance 
owing  to  the  finite  angle  between  the  forward  pump  and  probe  id 
The  experimental  NDFWM  showing  the  narrow  A  «  0  resonance 
[This  result  should  be  compared  with  Fig  4 '  b l . | 
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results  in  a  residual  Doppler  width  associated  with  the  5  =  0 
resonance.  The  dotted  line  in  Fig.  5(b)  shows  the  NDFWM 
response  corrected  for  this  effect.  Figure  5(c)  shows  the 
experimental  measurement.  The  second  resonance  has  a 
width  given  by  twice  the  total  dephasing,  as  expected  from 
Eq.  (8).  However,  the  5  =  0  resonance  is  narrower  and 
stronger  than  the  second  resonance,  showing  that  the  inter¬ 
ference  is  incomplete  and  enabling  the  resonance  associated 
with  the  ground-state  lifetime  to  be  observed. 

D.  Stationary  Two-  and  Three-Level  Atoms  above 
Saturation 

In  a  simple  two-level  system  in  the  presence  of  strong  reso¬ 
nant  optical  excitation,  the  levels  appear  to  split  because  of 
the  ac-Stark  effect,  resulting  in  a  three-peaked  fluorescence 
spectrum.1-  The  ac-Stark  effect  is  also  directly  observable 


a.  .  set  of  equations  that  includes  the  effects  of  perturba¬ 
tion  by  a  weak  probe: 


(V  =  [H0,p1]  +  [V'V]MV\pV<'>~ 

ot 


ph 


(10) 


The  first  set  of  equations  is  solved  exactly.  Then  the  results 
are  inserted  into  the  second  set,  which  is  solved  perturba- 
tively.  From  the  solution,  we  solve  for  the  polarization 
given  by  P  =  Tr  (up),  where  p  is  the  dipole-moment  operator. 
(The  precise  set  of  terms  to  be  included  in  the  polarization  is 
determined  by  the  phase-matching  conditions  that  give  rise 
to  a  signal  counterpropagating  with  respect  to  the  input 
probe.)  We  find  that  the  polarization  is  given  by 
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(5  -  i>  „)(5  -  />2r)([A  +  5] 
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as  a  splitting  of  the  NDFWM  spectrum.1’-1'’  The  criterion 
for  observing  the  splitting  in  any  kind  of  experiment  is  that 
the  Rabi  flopping  frequency  pE/h  must  exceed  the  natural 
linewidth  determined  by  the  total  dephasing  rate,  yphr.  In  a 
simple  closed  two-level  system,  this  condition  corresponds 
exactly  to  the  requirement  that  the  optical  intensity  must 
exceed  the  system  saturation  intensity.  However,  in  a 
three-level  system  with  a  metastable  third  level,  the  satura¬ 
tion  intensity  (which  is  observed  as  a  reduction  in  the  ab¬ 
sorption  coefficient)  can  be  reduced  several  orders  of  magni- 
t"de  because  of  optical  pumping-induced  depletion  of  the 
ground  state.  In  this  case,  there  is  no  direct  relationship 
between  the  ratio  of  the  Rabi  flopping  frequency  to  the  total 
dephasing  rate  and  the  ratio  of  the  intensity  to  the  satura¬ 
tion  intensity.  Hence  we  expect  to  see  saturation  effects  in 
the  NDFWM  spectrum  without  observing  ac-Stark  split¬ 
ting. 

Assuming  strong  pumps  and  a  weak  probe  beam,  we 
solved  the  equation  of  motion  for  the  density-matrix  opera¬ 
tor  by  using  the  same  approach  discussed  in  earlier  work  on 
saturable  absorbers  (see  Refs.  11,  36,  and  37). 

We  use  Eq.  (1)  with  the  above  modifications  to  describe 
Fig.  1(c)  and  write  Vas  V°  +  V,  where  V°  is  the  interaction 
potential  due  to  the  strong  counterpropagating  pump  waves 
and  V1  is  the  interaction  potential  due  to  the  weak  probe 
wave.  (We  ignore  the  effects  of  the  signal  wave  since  we  are 
assuming  that  the  nonlinear  response  is  sufficiently  weak  to 
result  in  a  phase-conjugate  reflectivity  far  less  than  unity.) 
The  density-matrix  equations  can  now  be  written  as  a  set  of 
zero-order  equations  because  of  the  strong  pumps: 


ihp"  =  |//,„  />"|  +  [I"1,  p")  +  ih 


<V' 

dt 


‘b.  <V 


ph 


(9) 


where  x<  is  the  Rabi  flopping  frequency  (pE.  t,)  associated 
with  field  £,,  xrt  =  x/exp(ikrx)  +  xt>exp(tk(,-xi.  ,J  -  vM 
+  y»p,  -Vphr  ®  r.>r/2  +  Tph.  5  is  the  pump-probe  detuning,  and 
A  is  the  pump-resonance  detuning. 

To  predict  the  macroscopic  nonlinear  response,  however, 
we  must  insert  the  polarization  into  the  wave  equation  and 
perform  a  spatial  average  of  the  polarization  since  the  term 
xrt  is  spatially  varying  (because  of  the  standing  wave  be¬ 
tween  the  forward  and  backward  pumps).  We  defer  discus¬ 
sion  of  this  full  solution  and  simply  discuss  the  physical 
behavior  described  by  the  polarization. 

The  signal  is  proportional  to  the  square  of  the  polariza 
tion,  and.  to  simplify  the  discussion  below,  we  set  the  cross 
term  in  lx*7!-  to  zero.  (This  is  a  serious  error  if  one  is 
attempting  to  predict  absolute  phase-conjugate  efficien¬ 
cies,11  but  unless  one  is  trying  to  make  quantitative  line- 
shape  measurements  it  does  not  qualitatively  affect  the 
NDFWM  spectral  response.) 

First  we  note  that  in  the  limit  of  >ph  =  >.-i  =  0.  two 
resonances  are  expected  as  a  function  of  5  at  5  =  ±  A.  As  the 
intensity  is  increased  and  saturation  is  observed,  the  levels 
undergo  ac-Stark  splitting,  and  three  resonances  are  expect¬ 
ed  at  5  =  0  and  5  =  ±  I x/?7! •  Figure  6(a)  shows  the  expected 
behavior,  in  agreement  with  the  discussion  in  Ref.  14. 

However,  in  the  limit  of  large  dephasing  (yPh  »  ■>  :r).  finite 
> _>!,  and  a  long-lived  metastabie,  we  expect  a  single  narrow- 
resonance  with  a  width  given  by  the  metastable  lifetime  if 
the  pump  intensities  are  much  less  than  the  three-level  satu¬ 
ration  intensity.  As  the  intensity  increases  above  satura¬ 
tion.  even  though  the  ac-Stark  splitting  is  still  small  com¬ 
pared  with  the  width  of  level  2,  as  anticipated  above,  we 
expect  the  NDFWM  spectral  response  to  broaden.  This 
effect  is  the  result  of  Fig.  6(b).  Most  importantly,  we  note 
that,  from  the  analytical  description,  we  expect  the  saturat¬ 
ed  width  to  broaden  linearly  with  the  pump  intensity. 
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Fig.  6.  (a)  The  NDFWM  response  for  a  closed  two-level  system. 
The  curve  with  two  resonances  is  the  low-intensity  ( xnT  *  0)  re¬ 
sponse.  similar  to  Fig.  2(a).  The  three-peaked  structure  is  the 
strong  Field  response,  showing  the  ac-Stark  splitting  and  a  central 
resonance  (xftr  “  8).  (b)  The  NDFWM  response  for  a  closed  three- 
level  system.  The  narrow  resonance  is  the  low-intensity  response 
1  Xht  *  0).  The  broader  repsonse  shows  the  expected  power  broad¬ 
ening  (xht  “  6.71).  For  this  calculation,  y.u  “  0.005,  >23  “  2,  -yPh  * 
10-',  and  A  -  0. 

3.  SUMMARY 

In  summary,  we  have  presented  an  analytical  description  of 
backward  NDFWM  to  extract  a  physical  understanding  of 
the  origin  of  the  spectral  response  in  relatively  simple  sys¬ 
tems.  The  analysis  included  a  model  for  state-specific  res¬ 
ervoir  coupling  and  spontaneous  emission  and  showed  that 
noncancellation  of  quantum-mechanical  amplitudes  domi¬ 
nated  the  experimental  observations. 


*  Present  address,  Department  of  Electrical  Engineering, 
University  of  Michigan,  EECS  Building,  Ann  Arbor,  Michi¬ 
gan  48109. 
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Frequency  domain  nonlinear  laser  spectroscopy  is  a  powerful  method 
for  studying  resonant  optical  excitation  in  solids.  It  provides  a  means 
for  obtaining  line  shapes  related  purely  to  relaxation  of  optically  excit¬ 
ed  states  of  the  system.  In  addition,  it  often  eliminates  contributions 
from  inhomogeneous  broadening  and  enables  studies  of  the  homogeneous  line 
shape  and  the  effects  of  spectral  diffusion.  '  ’  In  this  paper,  we  dis¬ 
cuss  measurements  of  the  exciton  in  a  GaAs/AlCaAs  multiple  quantum  well 
: MQW ) . 

GaAs  MQW' s  are  semiconductor  heterostructures  of  importance  to  elec¬ 
tronic  and  optoelectronic  devices*.  They  consist  of  a  periodic  structure 
of  thin  layers  of  GaAs  and  AlCaAs.  The  larger  band  gap  of  AlGaAs  confines 
the  optical  excitation  to  the  thin  GaAs.  resulting  in  quantum  confinement 
of  the  exciton  and  an  increase  in  the  exciton  binding  energy  In  this 
work,  the  MQW  is  grown  by  molecular  beam  epitaxy  and  consists  of  65  peri¬ 
ods  of  96 A  wells  of  GaAs  and  98A  barriers  of  A1  Ga  As.  The  solid  line 

a  7 

in  Fig.  1  shows  the  absorption  spectrum  at  SIC.  The  two  resonances  are  the 
heavy  hole  ( HH1  )  and  light  hole  (LH1)  exciton  which  exhibit  strong  inhomo¬ 
geneous  broadening  due  to  fluctuations  in  the  well  thickness 


Figure  l.  The  absorption 
spectrum  and  the  degenerate 
four-wave  mixing  response  at 
5K.  The  red  shift  in  the  non¬ 
linear  response  shows  the  non¬ 
linear  response  in  this  sample 
is  due  to  localized  excitons. 


The  primary  backward  four-wave  mixing  experimental  configuration  is 

discussed  in  reference  5  Two  counter-propagating  pump  beams  (identified 

oy  field  amplitudes  &  (u  ,£  )  and  £  (w  .£  )|  Interact  with  a  third  beam 
r  f  r  eat 

.<?  fu  .it  II  through  the  third  order  susceptibility  to  produce  a  nonlinear 

p  p  p 

polarization.  The  probe  field  k-vector  it  intersects  the  forward  oump 

P 

field  k-vector.  itf,  at  an  angle  8  The  signal  field  (designated. 

£  tu  ,«  II.  arises  from  the  induced  nonlinear  response  proportional  to 

»  t 

£  £  S  Since  the  frequencies  are  nearly  degenerate,  phase  matching  en- 

f  p  b 

sures  that  the  signal  Is  neafly  counter-propagating  with  respect  to  the 

probe  field.  Physically.  S  S  form  a  traveling  wave  spatial  grating  asso- 

f  P 
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dated  with  the  excitation  of  the  system  6  scatters  from  this  grating 

to  produce  the  signal.  The  nonlinear  response  in  these  systems  is  due  to 
band  filling  and  exchange  effects 

For  £  •  £  !|2  ,  measurement  of  IG  I2  as  a  function  of  u  -u  ~i  (for 

ty-  r 11  p  s  p  f 

u  -u  =constant),  designated  the  FUMp  response,  gives  a  measure  of  the  re- 

D  f 

laxation  rates  associated  with  the  optical  field  induced  perturbation  cf 
the  excited  states  and  the  ground  state,  provided  there  is  no  oriental  icr. 
signal  (i.e.,  E  =0  if  £  # S  \£  .  which  we  have  verified  experimer tai 1 y  for 

*  b  “  f-1-  p 

low  power  cw  excitation.)  For  fixed  6,  measurement  of  |6M  as  a  function 


of  u  -u  (designated  the  the  FUMb  response)  gives  the  line  shape  associ- 

b  f 

ated  with  the  nonlinear  response  in  homogeneously  broadened  material,  and 
a  line  shape  related  to  the  homogeneous  line  shape  in  inhomogeneousiy 
broadened  material 

At  300K,  the  “xciton  is  ionized  by  LO  phonons  on  a  time  scale  of 
300f sec  The  FUMp  response  measures  dynamics  associated  with  the 

electron-hole  (e-h)  plasma.  Figure  2a  shows  a  typical  FUMp  response.  The 
solid  _urve  is  a  fit  of  a  Lorentzian.  Dynamics  are  determined  by  e-h  re¬ 
combination  and  spatial  diffusion.  The  spatial  diffusion  rate  depends  on 
the  grating  spacing  (AK«lk  -k  I)  and  hence  the  angle.  6.  Figure  2b  show 

the  angle  dependence.  The  solid  curve  is  a  fit  of  the  form  7r^*D(£K) 

The  data  gives  r  («5nsec)  and  D  ( *18cm2/sec ) .  The  measurements  of  the 
_  r«e 


Figure  2  (a)  The  FUMp 
response  at  room  temper¬ 
ature  (b)  The  angle 
dependence  of  the  FUMp 
response. 


room  temperature  nonlinear  response  show  the  presence  of  a  slow  component. 
Figure  3a  shows  a  high  resolution  display  of  the  center  of  Fig  2  The 
line  shape  shows  a  clear  effect  of  interference  in  the  scattering  channel 
giving  rise  to  the  signal.  Phenomenological  modeling  of  the  nonlinear 
response  shows  that  the  Interference  is  likely  due  to  a  shift  in  the 


Figure  3  A  high  resolution  scan  of  the  FUMp  response  (a) 

u  <<u  (b)  and  (d)  are  the  results  of  the  model, 

s  r 
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resonance.  Figure  3b  shows  the  predicted  behavior.  In  good  agreement  with 
the  experiments.  Further  confirmation  of  the  explanation  is  seen  by  com¬ 
paring  Fig  3c  and  the  corresponding  model  In  3d.  In  this  case,  the  back 
pump  beam  is  detuned  to  the  low  energy  side  of  the  resonance,  eliminating 
the  interference  effect.  The  origin  of  the  shift  may  be  due  to  a  renor¬ 
malization  of  the  exclton  energy  or  a  shift  In  the  exciton  energy  due  to 
thermal  effects. 

At  low  temperature,  the  exclton  Is  stable  against  ionization,  and  the 
dashed  line  In  F'.g  1  shows  the  nonlinear  response  (measured  by  degenerate 
FUH)  is  red  shifted  with  respect  to  the  absorption  spectrum.  The  recom¬ 
bination  time  of  the  exciton  is  1.5  nsec,  so  the  FUMp  response  would  be 
expected  to  give  a  line  shape  with  a  width  determined  by  the  inverse  life 
time.  Figure  4a  shows  a  much  more  complex  structure.  The  line  shape 
shows  the  presence  of  two  clear  decay  channels,  corresponding  to  ;oo  psec 
and  15  nsec.  In  addition,  a  high  resolution  scan  of  the  center  shows  a 
narrow  resonance  similar  to  Fig.  3a.  Figure  4b  shows  an  expanded  scan  of 


(B) 


■A 


j  Figure  4.  (a )  FUMp 
response  at  5K 
!  (b)  shows  an  ex- 
,  panded  range  along 
'with  a  curve  f  i  t 
corresponding  to 
.  the  1 . 5  nsec  recom¬ 
bination  time 


A 


Forward  Pump- Probe  Detuning  Forward  Pump-Probe  Detuning  (GHz.) 


the  Intermediate  region  along  with  a  curve  fit  corresponding  to  IS  nsec. 

*  We  believe  the  1.5  nsec  decay  Is  due  to  radiative  recombination  and  -he 
100  psec  decay  corresponds  to  the  decay  due  to  spectral  diffusion.  Solu¬ 
tions  of  the  optical  Bloch  equations  for  the  FUMp  response  modified  to 
describe  spectral  diffusion  (M08EJ3  show  the  FWMp  response  would  have  the 
form  described  by  Fig  4b  and  the  base  of  4a.  The  origin  of  the  ISnsec 
component  Is  possibly  due  to  the  presence  of  residual  electric  fields 
which  can  result  In  an  increase  in  the  exclton  lifetime. 

Further  evidence  for  spectral  diffusion  Is  seen  in  the  FUMb  response 
which  eliminates  Inhomogeneous  broadening  and  provides  a  line  shape  relat¬ 
ed  to  the  homogeneous  line  shape  along  with  contributions  from  spectral 
diffusion.  Figure  5a  shows  the  FVMb  response  obtained  28  cm  1  below  the 
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Figure  5.  The  FUMb  re¬ 
sponse  at  5K  (a) 

28cm  below  the  absorp¬ 
tion  maximum.  Ib)  5cm 
below  the  absorption 
maximum. 


absorption  maximum.  This  curve  shows  slight  asymmetry  on  the  high  energy 
side.  The  asymmetry  becomes  more  dramatic  as  the  forward  pump  and  probe 
frequency  approach  the  absorption  resonance,  as  seen  5b  It  can  be  shown 
that  in  the  absence  of  spectral  diffusion,  the  FUMb  response  is  symmetric 
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about  u  -u  (=u  1.  The  presence  of  asymmetry  is  consistent  with  excitation 

b  (  p 

diffusion  to  other  resonant  frequencies,  and  such  line  shapes  can  easily 
be  obtained  from  MOBE.  A  possible  origin  of  the  spectral  diffusion  is 
variable  range  hopping.  Recent  work  has  identified  phonon  assisted 
tunneling  and  thermal  activation  as  two  possible  mechanisms.  Figure  6 
shows  a  comparison  of  the  measured  temperature  dependence  of  the  lir.» 
width  in  Fig.  5a  to  a  theory  based  on  phonon  assisted  tunneling.  The 

curve  varies  as  TexpfB/kT'  )  and  is  in  general  agreement  with  theoret- 
wi  i - - 


Figure  6.  The  FVMb  line  width  ir. 
Fig.  5a  as  a  function  of  tempera¬ 
ture.  The  solid  curve  is  a  fit 
of  the  theory  for  phonon  ass.stea 
tunneling. 


Temperature  I K) 


ical  expectations1 * 3 * 5 * 7  Finally  the  effects  of  spectral  diffusion  which  w« 
have  observed  are  more  clearly  seen  at  higher  temperatures  Figure  'a 
shows  the  FVMb  response  at  5K  similar  to  Fig.  5b  with  an  increase  in  scar. 


Figure  7  The  FVMb  response  as  a  function  of  temperature.  (a)  5K  ibi 
15K  (c)  Solution  to  MOBE. 

At  15  K.  we  clearly  observe  the  decay  of  the  excitation  to  lower  energy 
In  these  measurements,  we  see  that  while  the  exciting  radiation  is  cen¬ 
tered  at  the  peak  of  the  absorption  spectrum,  the  excitation  clearly 
decays  to  the  region  of  the  localized  exciton.  The  large  broad  peak  is 
centered  at  the  same  wavelength  as  the  peak  of  the  nonlinear  response 
shown  in  Fig.  i  and  also  corresponds  to  the  peak  in  photolummescence 
emission  Fig.  7c  is  a  solution  of  the  MOBE  for  the  FVMb  response  in  the 
presence  of  decay  to  a  lower  lying  state. 
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We  describe  frequency-domain  four-wave  mixing  spectroscopy  measurements  near  the  fundamental  band  edge  in 
room-temperature  GaAs  multiple  quantum  wells.  The  line-shape  information  gives  a  measure  of  the  dynamical 
behavior  of  the  material  and  what  is  to  our  knowledge  the  first  observation  of  an  interference  line  shape  due  to  a 
dominant  slow  contribution  to  the  nonlinear  response. 


Using  frequency-domain  four-wave  mixing  (FWM), 
we  report  nonlinear  laser  spectroscopy  studies  of 
room-temperature  in  GaAs/AlGaAs  multiple-quan¬ 
tum-well  (MQW)  structures  in  the  vicinity  of  the  n  =  1 
exciton  features.  These  materials  have  been  studied 
in  recent  years  because  of  their  large  nonlinear  re¬ 
sponse,  which  is  of  importance  to  optoelectronic  appli¬ 
cations.  1  These  studies  are  also  of  fundamental  inter¬ 
est  since  the  optical  properties  reflect  the  effects  of 
quantum  confinement. 

Experimental  and  theoretical  research  suggests  that 
the  origin  of  the  room-temperature  nonlinear  response 
near  the  exciton  resonances  is  due  to  band  filling  and 
exchange  effects  caused  by  the  electron-hole  (e-h) 
plasma  produced  by  the  rapid  thermal  ionization  of 
the  exciton.  In  a  simple  picture  the  linear-optical 
polarization  is  proportional  to  Y.nfJ(u  -  +  iyn), 

where  the  sum  is  over  all  exciton  states.2  The  parame¬ 
ters  wn,  yn,  and  /„  represent  the  frequency,  width,  and 
oscillator  strength,  respectively.  When  the  exciton  is 
formed,  it  ionizes  on  a  time  scale  of  300  fsec,  and 
theory  has  shown  that  the  resultant  e-h  plasma  signif¬ 
icantly  modifies  the  exciton  oscillator  strength.  Ear¬ 
lier  experimental  work  on  differential  transmission2 
resulted  in  a  phenomenological  description  of  the  non¬ 
linear  response  that  described  the  experimental  mea¬ 
surements  if  the  e-h  plasma  modified  not  only  the 
oscillator  strength  but  also  the  linewidth  due  to  colli¬ 
sions.  A  small  shift  in  the  resonance  frequency  was 
also  required  to  produce  a  good  fit.  A  shift  of  this 
resonance  has  been  predicted  theoretically  as  a  result 
of  renormalization  of  the  exciton  resonance  frequency.4 

In  such  a  phenomenological  model  the  fundamental 
parameters  are  taken  to  be  the  value  in  the  absence  of 
the  optical  field  plus  a  term  that  is  linearly  dependent 
on  the  induced  e-h  plasma  density.  The  optical  re¬ 
sponse  can  then  be  expanded  about  the  zero-field  val¬ 
ue  to  obtain  a  third-order  polarization  at  frequency  u 
that  has  the  general  form  (w  =  w,  -  u,  + 

Wk)€i6*6k  exp(i(k,  -  k,  +  k*)  •  r  -  i(w,  -  +  ui*)f], 
where  Xu,\  >8  the  third-order  susceptibility.  Physical¬ 


ly  for  w,  *  Wj  near  an  exciton  resonance,  6t6*  creates 
the  e-h  plasma  modifying  the  optical  response  that  is 
sensed  by  the  third  field  6k-  Spectroscopic  measure¬ 
ments  are  made  as  a  function  of  any  of  the  three  input 
frequencies.  Information  can  be  obtained  about  spe¬ 
cific  energy-level  structure  as  well  as  relaxation  rates 
(i.e.,  dynamical  behavior)  associated  with  states  in¬ 
volved  in  the  transition.  In  this  Letter  we  describe 
measurements  of  the  line  shapes  associated  with  the 
room-temperature  susceptibility  that  provide  infor¬ 
mation  on  e-h  recombination,  ambipolar  diffusion, 
and  heavy-  and  light-hole  (hh)  and  (lh)  line  shapes. 
We  also  show  what  is  to  our  knowledge  the  first  obser¬ 
vation  of  a  narrow  interference  line  shape,  which  is 
evidence  of  a  shift  in  the  exciton  resonance  frequency. 

The  FWM  experimental  configuration  is  discussed 
in  Ref.  5.  Two  counterpropagating  pump  beams 
[identified  by  field  amplitudes  6f{w/,  k/)  and  <?6(u)b, 
kf,))  interact  with  a  third  beam  [6p(wp,  kp)]  through 
the  third -order  susceptibility  to  produce  a  nonlinear 
polarization.  The  probe  field  k  vector  kp  intersects 
the  forward  pump  field  k  vector  k/  at  an  angle  8.  The 
signal  field  [designated  6s(w„  ks)j  arises  from  the  in¬ 
duced  nonlinear  response  proportional  to  6f6p*6b- 

For  6b  -L  6 f  I!  6P,  measurement  of  |6'J2  as  a  function 
of  up  -  w/  =  <5  (for  w/  =  cot,  =  constant),  designated  the 
FWMp  response,  gives  a  measure  of  the  relaxation 
rates  associated  with  the  optical-field- induced  pertur¬ 
bation  of  the  excited  states  and  the  ground  state,6 
provided  that  there  is  no  orientation  signal"  (i.e.,  Es  - 
0  if  6b  II  Sf  _L  6P,  which  we  have  verified  experimental¬ 
ly  for  low-power  cw  excitation).  Experimentally,  one 
dye  laser  is  used  to  provide  the  pump  beams,  and  a 
second  dye  laser  provides  the  probe  beam.  For  mea¬ 
surement  of  line  shapes  narrower  than  2  MHz,  the 
interlaser  jitter  is  too  large,  and  FWMp  measurements 
are  made  with  one  laser  and  two  acousto-optic  modu¬ 
lators  using  the  method  of  correlated  optical  fields.6 
For  a  fixed  6,  measurement  of  1 6  J2  as  a  function  of  uh  - 
aif  (designated  the  FWMb  response)  gives  the  line 
shape  associated  with  the  nonlinear  response  in  homo- 
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Fig.  1.  Comparison  of  the  room-temperature  linear  absorp¬ 
tion  (solid  curve)  and  the  FWMb  response  (points). 


geneously  broadened  material  and  gives  the  line  shape 
related  to  the  homogeneous  line  shape  in  inhomogen- 
eously  broadened  material.5  Experimentally,  one  la¬ 
ser  provides  the  forward  pump  and  probe  beams  and  a 
second  laser  provides  the  back  pump  beam.  Acousto¬ 
optic  modulators  are  used  to  control  of  b. 

The  MQW  well  was  grown  by  molecular  beam  epi¬ 
taxy  on  a  GaAs  substrate  and  consists  of  alternate 
layers  of  9.6-nm  GaAs  wells  and  9.8-nm 
Alo.nGao-As  barriers  with  a  total  of  65  periods.  The 
incident  optical  beams  are  nearly  normal  to  the  growth 
direction.  All  measurements  were  done  at  low  power 
(1-5  W/cm2),  corresponding  to  an  exciton  density  in 
the  well  region  of  the  order  of  105  cm'2  and  a  carrier 
density  of  the  order  of  109  cm-2. 

In  the  first  measurements,  we  obtained  the  FWMb 
response.  As  shown  below,  by  adjusting  6,  it  is  possi¬ 
ble  to  control  the  relative  contributions  of  different 
states  involved  in  the  nonlinear  interaction  because 
these  states  are  characterized  by  different  excitation  - 
relaxation  rates.  To  reduce  contributions  from  slow 
terms  in  the  nonlinear  response,  we  set  b  -  500  kHz. 
Figure  1  shows  the  room-temperature  absorption 
spectrum  (the  solid  curve)  characterized  by  the  usual 
hh  and  lh  features  and  the  interband  transitions 
contributing  to  the  background.  Superimposed  onto 
Fig.  1  is  the  FWMb  response  (the  points).  The  spec¬ 
trum  clearly  shows  the  hh  and  lh  features  with  no 
background  contribution  from  interband  transitions. 
Such  measurements  are  similar  to  the  first  measure¬ 
ments  of  degenerate  FWM  (u>/  =  uib  =  wp)  in  a  MQW  at 
low  temperature9  and  at  room  temperature10;  howev¬ 
er.  by  tuning  only  the  back  beam,  the  relative  signal 
strengths  in  optically  thin  material  (normalized,  e.g., 
to  the  response  at  the  hh)  are  proportional  to  the  magni¬ 
tude  squared  of  the  optical  response  functions  of  the 
hh  and  lh  excitons  rather  than  the  more  complex  de¬ 
pendence  in  degenerate  FWM.  Physically,  the  origin 
of  the  behavior  is  qualitatively  understood.  The  for¬ 
ward  pump  and  probe  fields  interact  to  produce  a 
spatial  modulation  of  excitons  that  quickly  ionize, 
producing  a  spatial  modulation  (and  for  finite  b,  a 
traveling  modulation)  of  the  e-h  plasma.  Ej,  scatters 
from  this  grating,  and  scanning  of  the  backward  beam 
frequency  measures  the  basic  frequency  response  of 
the  linear  exciton  susceptibility. 


In  the  next  measurements,  we  examined  the  FWMp 
response,  which  measures  the  decay  rate  of  the  excita¬ 
tion.  The  decay  of  the  spatial  modulation  is  deter¬ 
mined  by  e-h  recombination  and  spatial  diffusion. 
The  spatial  diffusion  in  the  grating  is  determined  by 
ambipolar  diffusion  Da  =  2 Dr,  where  the  electron  mo¬ 
bility  is  much  greater  than  the  hole  mobility.”1  In 
this  limit  the  nonlinear  polarization  is  given  by 

PNL  =  i€f€p*Cb(ib  +  7  +  DjAkl2)-1 

X  exp(i6t  +  t'Ak  •  x)exp(-iuht  +  ikh  •  x)  +  c.c., 

where  f  is  a  constant,  7  is  the  decay  rate  due  to  carrier 
recombination,  lAkl2  =  Ik/  -  kp|2  =  (2k)'2  sin2(0,nl/2), 
with  0jm  the  internal  angle  between  the  forward  pump 
and  probe  beams,  and  we  have  taken  |fc/|  s  \kp\  =  k. 
Tuning  b  provides  information  on  the  relaxation  of 
these  states,  and  the  dependence  of  the  FWMp  spec¬ 
tral  response  on  Da  and  6  gives  a  way  to  measure  Da 
and  7.  A  corresponding  method  in  the  time  domain 
has  been  used  for  exciton  delocalization  studies.12 

There  are  two  regimes  for  this  experiment,  corre¬ 
sponding  to  b  -  up  -  o>  >  1  MHz  and  b  <  1  MHz  (u>  =  «/ 
=  «(,).  For  the  large  detuning  range,  we  find  line 
shapes  similar  to  that  shown  in  Fig.  2.  The  solid  curve 
is  a  fit  of  a  Lorentzian  curve,  showing  that  the  line 
shape  associated  with  the  decay  is  clearly  Lorentzian. 
As  anticipated  above,  this  line  shape  is  angle  depen¬ 
dent  as  shown  in  Fig.  3.  The  functional  form  of  the 
curve  is  given  by  7  +  Da( 2k)2  sin2(0,nt/2).  The  solid 
curve  is  a  fit  of  this  equation,  from  which  we  determine 
the  ambipolar  diffusion  coefficient  and  the  recombi¬ 
nation  time  of  the  e-h  plasma.  In  this  case,  we  find 
that  Da  -  18  cm2/sec  and  7-1  =  tr  -  5  nsec.  Da  is  in 
reasonable  agreement  with  that  expected  from  the 
measured  hole  mobility,11  and  both  Da  and  tr  are  in 
agreement  with  those  values  reported  by  Chemla  et 
al.3  using  time-domain  methods. 

For  b  <  1  MHz,  the  FWMp  spectrum  was  obtained 


Fig.  2.  FWMp  response.  The  solid  curve  is  a  least-squares 
fit  of  a  simple  Lorentzian.  Inset  (a)  is  a  high-resolution 
spectrum  of  the  tip  of  the  primary  curve  showing  an  interfer¬ 
ence  profile.  Insets  (b)  and  (c)  show  that  such  profiles  can 
be  seen  in  MOBE’s,  which  are  phenomenologially  modified 
to  account  for  the  nonlinear-optRal  response  in  semiconduc¬ 
tors.  Inset  (c)  is  an  enlargement  of  the  central  region  of 
inset  (b). 
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Fig.  3.  Linewidth  in  Fig.  2  as  a  function  of  the  angle  be¬ 
tween  the  forward  pump  and  probe  beams.  The  solid  curve 
is  a  least-squares  fit  of  the  linewidth  described  by  recombi¬ 
nation  and  ambipolar  diffusion. 

using  the  method  of  correlated  optical  fields.  Inset 
(a)  of  Fig.  2  shows  the  response,  a  profile  associated 
with  interference  effects.  The  relaxation  time  associ¬ 
ated  with  this  profile  is  of  the  order  of  10  jisec.  We 
interpret  these  data  as  being  the  result  of  interference 
between  multiple  scattering  channels  in  the  nonlinear 
interaction,  where  at  least  one  of  the  channels  is  asso¬ 
ciated  with  a  long-lived  excitation  of  the  system. 

Existing  theoretical  descriptions  have  not  described 
such  behavior,  and  the  origin  is  likely  to  be  complex. 
However,  there  have  been  recent  discussions  of  effec¬ 
tive  optical  Bloch  equations,  (EOBE’s)13  which  resem¬ 
ble  typical  density-matrix  equations  describing  opti¬ 
cal  interactions  [the  modified  optical  Bloch  equa¬ 
tions14  (MOBE’s)].  Without  solving  the  EOBE’s  it  is 
instructive  to  solve  the  MOBE’s,  where  these  equa¬ 
tions  are  further  phenomenologically  modified  to  ac¬ 
count  for  the  nonlinearity  in  semiconductors.  Recall¬ 
ing  that  the  origin  of  the  room-temperature  nonlinear 
response  is  due  to  the  presence  of  an  e-h  plasma  that 
modifies  the  fundamental  parameters  yn,  and  /„, 
we  assume  that  these  different  terms  are  modified 
according  to  the  prescription  £„  -*  £„(1  +  an£y<j,<5/), 
where  an  is  a  constant  that  describes  the  linear  depen¬ 
dence  of  a  given  parameter,  £„  (=>„,  or  /„),  on  the 
electron  and  hole  density.  If  we  assume  that  either 
the  electron  or  the  hole  is  trapped  in  a  long-lived  state 
such  as  an  interface  state,  then  there  exists  a  long- 
lived  effect  on  the  modified  parameters. 

Solving  the  corresponding  equation  for  the  FWMp 
response,  we  obtain  the  profile  shown  in  inset  (b)  of 
Fig.  2,  where  the  6  is  normalized  to  the  inverse  lifetime 
of  the  long-lived  state.  Inset  (c)  of  Fig.  2  is  an  enlarge¬ 
ment  of  the  interference  region.  The  interference  ef¬ 
fect  arises  between  the  e-h -modified  y  or  /  and  a  small 
but  necessary  modified  exciton  frequency  w.  Without 
the  shift  in  the  resonant  frequency  (owing  to  carriers 
or  some  other  laser-induced  mechanism),  the  model 
shows  just  a  narrow  resonance  with  no  interference 
dip.  The  interference  results  because  the  contribu¬ 
tion  of  the  frequency  shift  to  the  nonlinearity  is  shift¬ 
ed  in  phase  by  90°  from  the  contributions  to  the  non¬ 
linearity  due  to  shifts  in  fn  and  yn.  In  a  direct  mea¬ 
surement  of  the  FWM  response,  the  sign  of  the  shift 


cannot  be  determined  because  the  signal  is  propor¬ 
tional  to  the  square  of  the  difference  between  the  (n 
(and  >„)  term  and  the  frequency-shift  term.  Current 
research  is  based  on  use  of  a  heterodyne-detection 
measurement,  which  should  resolve  this  ambiguity. 
Furthermore,  in  this  material  (as  in  other  samples  that 
we  have  examined)  the  cw  nonlinear  response  is  domi¬ 
nated  by  the  slow  component.  The  slow  component 
gives  rise  to  a  signal  which  is  four  times  the  strength  of 
the  fast  response.  [This  was  determined  from  the  one 
laser  FWMp  measurement  shown  in  inset  (a)  of  Fig.  2 
by  taking  the  ratio  of  the  signal  at  the  peak  to  the 
signal  in  the  wings  of  the  narrow  resonance.]  This 
component  is  most  likely  electronic  in  origin  since  it 
persists  (with  the  same  width)  when  the  sample  is 
immersed  in  liquid  helium.5  Note,  however,  that  the 
room-temperature  behavior  could  also  be  the  result  of 
a  contribution  from  a  spatially  modulated  tempera¬ 
ture-induced  shift  in  the  resonance.  The  phenomeno¬ 
logical  calculation  shown  in  the  inset  of  Fig.  2  would  be 
similar  for  a  thermal  effect. 
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High  Resolution  Nonlinear  Laser  Spectroscopy  of  Excitation  and  Relaxation 
Near_  the  Band  Edge  in  CaAs  Quantum  Well  Structures 

D.G.  Steel,  P.K.  Bhattacharya,  J.T.  Remillard, 

Hallin  Wang,  M.D.  Webb,  J.  Pamulapati,  J.  Oh 

Departments  of  Electrical  Engineering  and  Physics 
Harrison  M.  Randall  Laboratory  of  Physics 
The  University  of  Michigan,  Ann  Arbor,  Michigan  48109 

The  linear  and  nonlinear  optical  properties  near  the  f undamental  band 
edge  of  GaAs  quantum  well  structures  are  important  for  applications  of  these 
materials  to  devices  as  well  as  providing  new  insight  into  the  effects  of 
quantum  confinement.  Excitons  dominate  the  optical  spectrum  in  this  region, 
however,  the  coupling  of  excitons  to  the  applied  radiation  field  is  modified 
by  dynamical  interactions  due  to  the  coupling  of  the  exciton  to  the 
surrounding  crystal  lattice  and  the  vacuum  radiation  field.  At  room 

temperature,  the  exciton  is  quickly  ionized  by  phonons  resulting  in  an 
electron-hole  plasma  which  modifies  the  optical  properties  due  to  the  effects 
of  bandfilling  and  exchange  that  lead  to  strong  optical  nonlinearities.1  At 
low  temperature,  the  exciton  is  more  stable  and  other  effects  such  as 
spontaneous  emission,  tunneling,  diffusion,  and  scattering  from  phonons  and 
defects  modify  the  optical  properties.  These  processes  can  result  in  decay  of 
the  excitation  as  well  as  decay  of  the  coherence  or  induced  polarization 

between  the  initial  and  the  excited  state.  In  this  paper,  we  describe  the  use 

of  high  resolution  nonlinear  laser  spectroscopy  based  on  four-wave  mixing 
(FWM)  to  obtain  lineshapes  associated  with  the  nonlinear  susceptibility.  The 

measurements  provide  new  understanding  of  the  physical  processes  associated 
with  the  relaxation  of  the  exciton  and  the  dynamics  of  the  optical  response. 
At  low  temperatures,  FWM  can  often  eliminate  inhomogeneous  broadening  leading 
to  a  direct  measurement  of  the  exciton  homogeneous  lineshape. 

The  experimental  configuration  for  these  measurements  is  based  on 
backward  four-wave  mixing  using  frequency  stabilized  tunable  cw  dye  lasers. 
Two  counterpropagating  pump  beams  with  fields  £  (w  )  and  2  (u  )  interact  with 

.  t  t  b  b 

a  probe  beam  E  (w  )  In  the  sample.  The  angle  between  the  forward  pump  and 
probe  beams  is  9.  The  three  beams  interact  through  the  nonlinear  response  to 
produce  a  signal  field  (£  )  arising  from  a  polarization  proportional  to 
X<3>(£f*£^)£b.  For  these  measurements,  £  |£^£  (no  signal  is  obtained  for 
eJeJ*,  for  low  power  cw  excitation),  and  the  signal  arises  from  the  coherent 
scattering  of  from  the  traveling  wave  grating  of  excitation  produced  by 
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£  •£*.  By  measuring  |E  I2  as  a  function  of  u  -w  for  fixed  d»u  -<j  (the  FWMb 

f  p  u  bf  p  I 

response,)  we  obtain  a  lineshape  associated  with  the  fundamental  excitation. 
2 

By  measuring  IE  I  as  a  function  of  5  for  fixed  u  (the  FWMp  response.)  we 

I  b 

obtain  a  lineshape  related  to  relaxation  rates  associated  with  the  excitation 

produced  by  the  forward  pump  and  probe. 

The  degenerate  FWM  response  (DFWM  where  u  ■«  )  using  cw  excitation  has 

f  b  p  23 

been  reported  for  both  low  temperature  and  room  temperature  ’  .  For 
comparison.  Fig.  1  shows  the  room  temperature  linear  absorption  and  FWMb 


Figure  1.  A  comparison 
of  the  linear  absorption 
spectrum  and  the  FWMb 
reponse. 


response  showing  two  features  corresponding  to  the  heavy  hole  (HH1)  and  light 
hole  (LH1).  The  ratio  of  the  HH1  to  LH1  FWMb  response  is  expected  to  vary  as 
the  oscillator  strength  squared,  in  contrast  to  the  more  complex  dependence  in 
DFWM.  In  these  measurememts,  6  is  fixed  at  500kHz  which  reduces  contributions 
from  slow  terms  in  the  nonlinear  response  that  (as  shown  below)  dominate  the 
cw  nonlinear  response  and  affect  the  ratio  of  the  HHI  to  LH1  response. 

At  room  temperature,  the  exciton  is  ionized  by  LO  phonons  and  the 

resultant  e-h  plasma  modifies  the  fundamental  parameters  of  the  exciton  such 

as  the  oscillator  strength,  resulting  in  the  nonlinear  response.  A 

measurement  of  the  FWMp  response  provides  a  linewidth  determined  by  the 

dynamics  of  the  e-h  plasma.  Figure  2  shows  a  typical  FWMp  spectrum,  obtained 

•2 

at  a  forward  pump-probe  angle  of  0«1O  rad.  The  solid  line  is  an  excellent 
least  squares  fit  of  a  simple  Lorentzian.  The  width  of  the  Lorentzian  is 
given  by  the  funtion  T+D#  I  £f-£  1 2  where  jr  is  the  carrier  recombination  rate 
and  D^  is  the  amblpolar  diffusion  coefficient.  Because  of  ambipolar 

diffusion,  the  linewidth  depends  on  angle,  as  seen  in  the  inset  (A].  The 
solid  line  is  a  fit  of  the  above  function  where  r  end  D  are  adjusted  to  fit 
the  data.  We  find  r*5nsec  and  D  »18cm2/sec. 

Using  the  method  of  correlated  optical  fields  to  eliminate  interlaser 
Jitter,  we  examined  the  FWMp  response  at  the  peak  of  the  main  curve  in  Fig.  2. 
The  inset  (B)  shows  the  observed  spectrum.  The  feature  is  extremely  narrow 
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(*30kHz)  corresponding  to  a  relaxation  time  of  lOpsec  and  shows  a  clear 
interference  feature  associated  with  this  profile.  We  note  further  that  the 
amplitude  associated  with  the  slow  feature  provides  a  nonlinear  signal  which 
is  a  factor  of  four  above  the  " normal "  faster  component.  Hence,  the 
importance  of  setting  6  to  be  large  in  the  FWMb  response  of  Fig.  lb  is  now 


Figure  2.  The  FWMp 
response.  The  solid  line 
is  a  Lorentzian.  Inset 
[A]  shows  the  angle 
dependence  along  with  a 
fit  of  the  expected 
response.  Inset  [B] 

shows  the  FWMp  reponse 
obtained  at  the  top  of 
the  main  curve. 


clear.  It  is  also  interesting  to  note  that  this  slow  component  persists  down 
to  liquid  helium  temperatures  and  is  clearly  electronic  in  origin.  Finally, 
we  observe  that  the  modified  optical  Bloch  equations  can  predict  such 
behavior.4  though  the  physical  origin  of  the  structure  is  not  understood. 

At  low  temperature,  the  lineshapes  and  physical  meaning  change 
significantly,  since  the  exciton  is  more  stable.  Figure  3  compares  two  FWMb 


Figure  3.  The  FWMb 
response  at  two  different 
excitation  wavelengths  at 
5  K.  (a)  is  20cm'  below 
the  HH1  absorption 
maximum,  (b)  is  at  the 
HH1  aborption  maximum. 


lineshapes  obtained  at  5K,  where  the  material  is  inhomogeneously  broadened  due 
to  static  fluctuations  in  well  thickness.  The  FWMb  response  yields  a  spectral 
response  given  by  the  convolution  of  the  homogeneous  llneshape  with  itself. 
Figure  shows  the  FWMb  response  obtained  20cm"1  below  the  HHl  absorption 
peak  (3 "500 kHz)  and  Fig.  3b  shows  the  response  at  the  peak  of  the  HHl 
absorption.  The  data  shows  the  evolution  of  a  high  energy  tail  due  to 
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contributions  from  slow  excitation  decay  terms  in  the  nonlinear  response, 
shown  below.  We  believe  this  data  represents  a  direct  measurement  of  the 
homogeneous  lineshape  of  the  exciton.  The  data  also  indicates  that  the  common 
assumption  that  the  homogeneous  lineshape  of  the  exciton  is  Lorentzian  may  not 
always  be  appropriate. 

The  FWMp  response  also  shows  complex  structure.  Figure  4  shows  the 
spectrum  is  clearly  characterized  by  two  resonances  situated  about  5*0.  In 
fact,  more  careful  study  shows  additional  structure  is  present,  consistent 
with  at  least  2  additional  resonances  within  the  bandwidth  of  the  above 
spectrum,  plus  another  resonance  which  is  observed  using  the  method  of 


correlated  optical  fields.  Based  on  other  work,  we  assign  the  broad  resonance 
at  the  base  of  Fig.  4  as  arising  from  radiative  recombination,  corresponding 
to  a  radiative  recombination  time  of  order  iOOpsec.  This  number  is  consistent 
with  other  recent  reports  obtained  at  comparable  optical  intensity5. 

This  work  is  supported  by  the  U.S.  Army  Research  Office,  the  Air  Force 
Office  of  Scientific  Research,  and  the  National  Science  Foundation. 
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Optical  Phase  Conjugation  and  Nonlinear  Optical 
Bandpass  Filter  Characteristics  in  CdSSe 
Microciystallite-Doped  Glass 

J.  T.  REMILLARD,  H.  WANG,  M.  D.  WEBB,  and  D.  G.  STEEL 


ittXrur-Thb  paper  describes  experimental  measurements  of  the 
backward  CW  four-wave  mixing  response  in  semiconductor-doped 
glass  in  a  regime  important  to  applications.  In  particular,  we  show  that 
lowering  the  operating  temperature  to  125  K  yields  a  yl3>  »  2  x  10~4 
esu.  Aberration  correction  through  optical  phase  conjugation  is  dem¬ 
onstrated  with  a  phase  conjugate  reflectivity  of  0.3  percent  along  with 
a  wide-angle  fleld-of-vlew  response.  In  addition,  we  report  on  the 
nearly -degenerate  four-wave  mixing  response  of  importance  to  tunable 
bandpass  optical  Alters. 

SEMICONDUCTOR-DOPED  glasses  (SDG)  have  been 
the  subject  of  considerable  recent  research  because  of 
their  potential  importance  to  new  nonlinear  optical  de¬ 
vices  and  because  they  provide  a  material  for  the  study  of 
quantum  confinement  effects.  These  glasses  are  doped 
with  microcrystallites  of  CdS,  _,Sex  and  are  commercially 
available  as  sharp-cut  color  glass  filters  [1].  The  com¬ 
mercial  material  is  characterized  by  a  wide  distribution  in 
microcry stallite  size  12]  (typically  100  A  in  diameter  with 
a  FWHM  size  distribution  of  50  A)  resulting  in  a  fea¬ 
tureless  absorption  edge  at  room  temperature.  (At  low 
temperature,  i.e.,  5.4  K,  these  materials  have  shown 
spectral  features  attributed  to  quantum  confinement  in  the 
luminescence  [3].)  However,  special  preparation  proce¬ 
dures  have  enabled  one  group  [4]  to  produce  structures  as 
small  as  12  A  in  CdS,  CdSe,  and  CuCl  and  another  group 
[2]  has  made  glasses  doped  with  microcrystallites  where 
the  size  distribution  is  controlled  by  a  heat  treatment  pro¬ 
cess  and  microcrystallites  as  small  as  25  A  have  been 
observed.  These  materials  show  clear  structure  in  the  room 
temperature  optical  absorption  edge,  believed  to  be  due 
to  the  effects  of  quantum  confinement  (2],  [4],  (5]. 

The  importance  of  SDG  to  device  applications  for  op¬ 
tical  switches  and  integrated  optics  has  been  recognized 
by  numerous  groups.  Very  early,  the  material  was  used 
to  demonstrate  a  Q  switch  in  a  laser  [6}.  Optical  bistabil¬ 
ity  was  demonstrated  by  McCall  and  Gibbs  [7]  and  Gibbs 
et  al.  [8);  however  the  slow  response  suggested  the  non- 
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86-0124.  The  work  of  M.  D.  Webb  was  supported  by  a  University  Re¬ 
search  Initiative  Fellowship  from  the  Army  Research  Office. 
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linearity  was  thermal  in  origin  [9].  Using  a  ^-switched 
laser,  Jain  and  Lind  [10]  demonstrated  optical  phase  con¬ 
jugation  and  showed  that  a  large  nonlinear  optical  suscep¬ 
tibility  was  observable  (x(3)  »  10"*  esu)  using  transient 
excitation.  They  determined  that  the  relaxation  time  was 
faster  than  8  ns  (determined  as  an  upper  limit.)  The  model 
for  the  nonlinear  response  proposed  by  Jain  and  Lind  was 
based  on  the  production  of  a  short-lived  electron-hole 
plasma.  This  work  was  followed  by  numerous  other  ex¬ 
perimental  [11  ]— [20]  and  theoretical  [2 1  ]— [23]  reports  of 
the  nonlinear  response  using  pulsed  excitation.  Included 
in  these  observations  have  been  several  measurements  of 
the  dynamical  behavior  [24]-[27]  of  key  importance  to 
high-speed  applications.  These  measurements  show  car¬ 
rier  relaxation  times  are  in  the  picosecond  time  domain. 
More  recently,  using  the  methods  of  ultrafast  laser  spec¬ 
troscopy,  measurements  of  transient  absorption  have 
shown  relaxation  processes  occurring  on  the  time  scale  of 
a  few  hundred  femtoseconds  [28],  [29]. 

The  large  nonlinear  optical  susceptibility  and  fast  re¬ 
sponse  times  have  been  important  for  several  recent  ap¬ 
plications.  In  particular,  25  ps  switching  in  optical  bista¬ 
bility  has  been  reported  [30].  In  addition,  the  use  of  this 
material  for  nonlinear  wave  guides  [31]-[34],  of  impor¬ 
tance  to  optical  processing  [35]— [37]  has  been  examined. 

In  this  paper,  we  examine  the  application  of  a  slower 
but  much  larger  nonthermal  contribution  to  the  nonlinear 
optical  response  obtained  under  CW  excitation.  The  first 
indication  of  the  possibility  of  a  large  CW  nonlinear  re¬ 
sponse  was  provided  by  the  observation  of  low-frequency 
dynamical  behavior  leading  to  subharmonic  generation 
and  near  chaotic  behavior  in  SDG  observed  at  low  tem¬ 
peratures  (140  K)  by  Zhedulev  et  al.  [38].  The  room  tem¬ 
perature  CW  nonlinear  response  was  measured  by  our 
group  using  DFWM  and  found  to  have  x<3>  »  I0'7  esu 
and  a  narrow  resonance  (4.4  kHz)  in  the  nearly-degener- 
ate  four-wave  mixing  (NDFWM)  response  corresponding 
to  a  72  fis  excitation  relaxation  time.  The  origin  of  the 
nonlinear  dynamics  is  attributed  to  contributions  from 
traps  and  is  the  subject  of  a  more  detailed  experimental 
study  to  be  presented  elsewhere  [40]  which  emphasizes 
the  physics  of  the  nonlinear  process  and  the  effect  on  trap 
dynamics  of  temperature,  intensity,  and  energy  transport. 
However,  in  the  present  paper,  we  show  that  reducing  the 
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temperatur  to  125  K  has  a  dramatic  effect  on  x(3>.  in¬ 
creasing  the  value  to  2  x  10-4  esu.  Based  on  this  result, 
we  demonstrate  behavior  important  to  optical  phase  con¬ 
jugation  [41],  including  wide-angle  field-of-view  reflec¬ 
tivities  as  large  as  0.3  percent,  sufficient  for  self-pumped 
ring  resonators  [42],  and  tunable  optical  filtering  [43], 
[44]  with  band  widths  (measured  below  saturation)  on  the 
order  of  a  few  hundred  hertz. 

The  experimental  configuration  consists  of  two  coun- 
terpropagating  pump  beams  £/  and  Eb  and  a  probe  beam 
Ep,  which  interact  with  the  two  pump  beams  in  the  glass. 
The  glass  is  a  standard  3  mm  thick  commercial  sample 
(RG630  from  Schott,  with  similar  results  observed  in 
Coming  material  [45]).  The  angle  between  the  forward 
pump  and  probe  is  small  enough  to  ensure  that  there  is 
complete  overlap  between  the  three  beams.  All  three 
beams  are  focused  by  a  1  m  focal  length  lens  with  a  spot 
area  of  approximately  8  x  10" 3  cm2.  The  signal  is  pro¬ 
duced  by  the  induced  nonlinear  polarization  proportional 
to  x(NL,(£)  EfEfEb  where  the  field  dependencce  of 
XNL(  £ )  is  intended  to  represent  the  possibility  of  satura¬ 
tion  of  the  four-wave  mixing  nonlinear  response  due  to 
higher  order  effects.  (Clearly,  under  saturation  condi¬ 
tions,  the  susceptibility  is  no  longer  third  order.)  As  is 
well  known,  the  phase-matching  conditions  result  in  a 
signal  which  is  counterpropagating  with  respect  to  the 
probe  beam.  The  signal  was  detected  by  either  a  PMT  or 
an  ordinary  power  meter  (e.g.,  Coherent  212).  Physi¬ 
cally,  the  signal  is  a  result  of  the  scattering  of  the  back¬ 
ward  (forward)  pump  beam  from  the  spatial  modulation 
of  the  absorption  and  dispersion  produced  by  the  forward 
(backward)  pump  and  probe  beam. 

The  effect  of  temperature  on  the  strength  of  the  nonlin¬ 
ear  response  is  shown  in  Fig.  1,  where  the  phase  conju¬ 
gate  reflectivity  is  given  as  a  function  of  incident  pump 
intensity  (below  saturation)  for  two  different  tempera¬ 
tures.  The  phase  conjugate  reflectivity  is  defined  as  the 
ratio  of  the  signal  power  to  the  probe  power.  Corrections 
have  been  made  for  the  losses  at  the  windows  of  the  dewar 
and  the  reflection  loss  of  the  sample.  The  solid  line  in  the 
figure  has  a  slope  of  two,  corresponding  to  a  quadratic 
dependence  on  pump  intensity.  The  nonlinear  response 
saturates  as  a  function  of  inverse  temperature  below  125 
K.  As  is  evident  in  Fig.  1,  the  decrease  in  temperature 
results  in  a  considerable  increase  in  the  nonlinear  re¬ 
sponse  as  is  evident  by  comparing  the  reflectivities  and 
corresponding  pump  intensities  in  Fig.  1(a)  to  the  much 
higher  reflectivities  and  lower  pump  intensities  of  Fig. 
1(b).  The  corresponding  values  for  x'3'  can  be  obtained 
using  the  expression  [46]  (in  CGS  units)  R  =  [(2 tw/cfi) 
X(3,]‘  (8ir/cfi  )2  L2lfIb,  where  R  is  the  reflectivity,  L  is  the 
interaction  length,  n  is  the  index  of  refraction,  c  is  the 
speed  of  light,  w  is  the  frequency  of  the  optical  field 
(rad/s),  and  /„  is  the  intensity  (erg/cm2/s)  of  the  a  pump 
beam  (no  local  field  corrections  are  made  and  this  suscep¬ 
tibility  value  includes  the  degeneracy  factor.)  We  find  that 
Fig.  1(a)  corresponds  to  x<3)  *  10”7  e*u  and  Fig-  1(b) 
corresponds  to  x<3'  *  10”4  esu. 


Effective  Pump  Intensity  (Wette/eq.  cm.) 


Effective  Pump  Intensity  (W»tts/sq.  cm.) 

Fig.  1.  The  OPC  reflectivity  is  a  function  of  incident  pump  intensity,  (a) 
294  K  (b)  125  K.  The  solid  lines  have  a  slope  of  2,  corresponding  to 
the  unsaturated  third-order  response. 

The  model  discussed  for  the  CW  nonlinear  response 
[39]  provides  a  physical  explanation  for  this  strong  tem¬ 
perature  dependence.  Namely,  the  forward  (backward) 
pump  and  probe  couple  a  valence  band  to  conduction  band 
transition  and  produce  a  spatial  modulation  in  the  con¬ 
duction  band.  Some  fraction  of  this  excitation  decays  to 
a  deep-level  trap,  resulting  in  a  spatial  modulation  of  the 
absorption  and  dispersion  associated  with  the  populated, 
trap.  The  strength  of  the  spatial  modulation  of  the  trap  is 
determined  by  the  rate  of  filling  of  the  trap  and  the  trap 
lifetime.  The  trap  lifetime  is  determined  by  thermal  acti¬ 
vation,  hence,  at  lower  temperatures  the  trap  lifetime  is 
longer,  resulting  in  a  stronger  spatial  modulation  of  the 
traps  and  a  higher  scattering  efficiency  of  the  backward 
(forward)  pump. 

However,  the  nonlinear  response  saturates  at  high  pump 
intensity,  as  shown  in  Fig.  2.  Currently,  we  believe  that 
the  onset  of  this  saturation  is  the  result  of  photoionization 
of  the  optically-excited  trap.  This  conclusion  is  suggested 
because  we  have  demonstrated  that  the  NDFWM  band¬ 
width,  which  is  determined  by  the  excitation  relaxation 
rate,  increases  linearly  with  pump  intensity  due  to  pho¬ 
toionization  of  the  trap.  At  low  intensity,  the  NDFWM 
bandwidth  (obtained  by  measuring  the  phase  conjugate 
reflectivity  as  a  function  of  the  detuning  of  the  probe  beam 
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Fig  2.  The  OPC  reflectivity  obtained  at  125  K  at  high  pump  intensity. 
The  data  show  the  saturation  of  the  nonlinear  response. 


frequency  with  respect  to  the  fixed  pump  beam  frequen¬ 
cies)  is  determined  by  the  intrinsic  trap  lifetime.  How¬ 
ever.  at  high  intensity,  the  trap  lifetime  is  shortened  by 
photoionization.  In  RG630,  we  see  that  the  maximum  re¬ 
flectivity  at  this  temperature  is  of  order  0.3  percent  but  is 
limited  by  this  saturation.  It  is  possible  that  engineering 
of  the  microcrystallites  may  be  able  to  decrease  the  pho- 
toio'nization  cross  section,  thus  enabling  larger  reflectiv¬ 
ities.  ^ 

Operating  at  the  high  reflectivity,  a  simple  qualitative 
demonstration  of  optical  phase  conjugation  was  made 
using  a  thermally-distorted  piece  of  glass  as  an  aberrator. 
In  the  usual  procedure,  the  far  field  of  the  probe  beam  was 
obtained  at  the  focus  of  a  2  m  focal  length  lens.  Fig.  3(a) 
shows  the  expected  near-diffraction-limited  spot.  A  single 
pass  through  the  aberrator  results  in  the  far-field  pattern 
shown  in  Fig.  3(b).  However,  Fig.  3(c)  shows  the  cor¬ 
rected  beam  obtained  by  allowing  the  phase  conjugate 
signal  produced  by  the  aberrated  probe  beam  to  retrace  its 
path  back  through  the  aberrator. 

In  addition  to  aberration  correction,  it  was  suggested 
very  early  that  based  on  optical  phase  conjugation,  the 
NDFWM  response  could  be  used  for  a  bandpass  optical 
filter  [43].  The  filtering  action  occurs  because  the  phase 
conjugate  reflectivity  depends  on  the  detuning  of  the  probe 
frequency  with  respect  to  the  pump  frequency.  In  the 
original  proposal,  the  material  was  assumed  to  have  an 
instantaneous  excitation  relaxation  rate  and  the  filter 
bandwidth  was  determined  by  phase-matching  restrictions 
on  the  detuning.  However,  the  first  demonstration  of  a 
filter  was  based  on  atomic  sodium  [44]  and  the  filter  band¬ 
width  was  determined  by  the  excitation  relaxation  rate  set 
by  spontaneous  emission.  Unfortunately,  for  CW  appli¬ 
cations,  the  field  of  view  of  devices  based  on  collisionless 
gas  phase  systems  is  limited  by  thermal  motion.  Optical 


Fig  3.  A  qualitative  demonstration  of  aberration  correction  though  CW 
OPC  in  semiconductor-doped  glass  la)  Far-field  intensity  distribution  of 
the  incident  probe  beam  (bl  Far-field  intensity  distribution  of  the  pr-be 
beam  after  a  single  pass  through  the  aberrator  (cl  Demonstration  of  the 
correction  of  the  aberration  in  ibl  after  reflection  from  the  OPC  and  pass¬ 
ing  back  through  the  aberrator 


filter  characteristics  in  materials  such  as  Cr:  YAICK  [47] 
also  have  limited  field  of  view  because  of  energy  diffu¬ 
sion.  The  field  of  view  is  determined  by  the  decrease  in 
the  signal  as  the  angle  between  the  forward  pump  and 
probe  increases. 

To  understand  why  thermal  motion  and  energy  diffu¬ 
sion  limit  the  field  of  view  of  any  device  based  on  optical 
phase  conjugation  (including  phase  conjugate  mirrors.)  we 
note  that  as  the  angle  increases,  the  period  of  spatial  mod- 
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ulation  decreases.  In  collisionless  gas  phase  systems,  if 
the  transit  time  (determined  by  the  spatial  period  spacing 
divided  by  the  thermal  velocity)  is  comparable  to  or 
smaller  than  the  spontaneous  emission  lifetime,  the  spa¬ 
tial  modulation  “washes  out,”  resulting  in  a  decrease  in 
the  scattering  efficiency.  For  small  angles,  the  snaring  is 
large  and  the  transit  time  is  long  compared  to  the  spon¬ 
taneous  emission  lifetime,  and  there  are  no  washout  ef¬ 
fects.  However,  for  large  angles,  the  transit  time  is  re¬ 
duced  considerably  and  washout  occurs,  lowering  the 
signal.  Alternatively,  for  solid-state  systems,  the  rate  of 
energy  diffusion  depends  on  the  gradient  of  the  spatial 
modulation  of  the  excitation.  As  in  the  gas  phase  system, 
if  the  energy  diffusion  time  is  comparable  to  the  state  life¬ 
times,  the  strength  of  the  spatial  modulation  is  reduced, 
again  reducing  the  scattering  efficiency.  For  small  angles, 
this  is  usually  not  important.  But  for  large  angles,  the  gra¬ 
dient  increases  and  the  energy  diffusion  time  decreases, 
again  resulting  in  washout  and  reduced  scattering  effi¬ 
ciency.  Hence,  for  optical  phase  conjugate  mirror  appli¬ 
cations,  these  effects  limit  the  field  of  view.  However,  in 
addition,  thermal  motion  and  energy  diffusion  both  con¬ 
tribute  to  the  rate  of  excitation  relaxation.  Hence,  the 
NDFWM  response  and  the  corresponding  filter  properties 
also  become  angle  dependent,  with  the  NDFWM  band¬ 
width  increasing  with  larger  angles.  It  is  interesting  to 
note  that  optical  phase  conjugation  based  on  photorefrac- 
tive  properties  in  insulators  are  not  characterized  by  large 
diffusion  coefficients  and  therefore  have  a  wide  field  of 
view  and  high  efficiency  [48],  but  are  often  limited  by 
extremely-small  bandwidths  (e.g.,  1  Hz). 

In  the  semiconductor-doped  glass,  our  earlier  room 
temperature  measurements  show  that  the  NDFWM  band¬ 
width.  obtained  by  measuring  the  reflectivity  as  a  function 
of  pump-probe  detuning,  provided  a  useful  bandwidth  on 
the  order  of  4.4  kHz,  which  is  determined  by  the  inverse 
trap  lifetime.  In  addition,  because  the  microcrystallites  are 
surrounded  by  glass  insulator  and  the  microcrystallites  are 
much  smaller  than  the  grating  spacing  associated  with  the 
nonlinear  response,  it  was  pointed  out  [10]  that  there 
should  be  no  energy  diffusion  and  no  angle  dependence  in 
the  nonlinear  response,  a  result  confirmed  by  our  room 
temperature  measurements. 

To  study  the  optical  filter  properties  of  SDG  we  first 
note  that  as  the  temperature  is  decreased,  the  NDFWM 
bandwidth  decreases,  which  as  we  indicated  earlier,  is  due 
to  the  fact  that  the  trap  lifetime  increases.  The  physical 
reason  for  this  is  that  when  the  probe  frequency  is  detuned 
from  the  pump  frequency,  the  fringes  move  through  the 
medium.  The  time  for  the  spatial  modulation  to  move  a 
half  period  is  of  order  6"'  (6  is  the  pump-probe  detuning 
frequency.)  If  that  time  is  long  compared  to  the  excitation 
relaxation  time  (trap  lifetime),  then  there  is  no  loss  of 
scattering  efficiency  for  the  pump  beam.  However,  if  6  is 
large  compared  to  the  inverse  relaxation  time,  the  spatial 
modulation  is  washed  out  by  the  moving  fringes,  thus  re¬ 
ducing  the  strength  of  the  spatial  modulation  and  the  scat¬ 
tering  efficiency.  Note  if  the  excitation  decay  rate  is  7. 


then  the  NDFWM  bandwidth  (in  hertz)  is  (7/x).  There¬ 
fore,  the  NDFWM  bandwidth  provides  a  measure  of  the 
real  time  response  capability  of  the  material. 

To  study  the  filter  properties,  we  measured  the 
NDFWM  bandwidth  using  the  method  of  correlated  op¬ 
tical  fields,  a  technique  useful  when  the  bandwidth  is  less 
than  the  interlaser  jitter  in  a  two  laser  measurement  [47], 
In  this  measurement,  the  pump  beams  and  the  probe  beam 
are  the  first-order  Bragg-scattered  beams  from  two  sepa¬ 
rate  acoustooptic  modulators  driven  by  two  phase-locked 
frequency  synthesizers.  In  this  way,  the  frequency  of  the 
probe  beam  can  be  independently  varied  with  respect  to 
the  pump  beams.  The  detuning  is  given  by  6  =  upn>tK  - 
“pump-  By  using  two  modulators  (one  for  the  pumps  and 
one  for  the  probe)  it  is  possible  to  scan  through  zero  de¬ 
tuning. 

The  low  intensity  (unsaturated)  NDFWM  bandwidth  is 
shown  in  Fig.  4  as  a  function  of  excitation  frequency.  The 
plot  shows  that  at  125  K,  narrow  bandwidths  (ranging 
from  2 1 5  to  430  Hz)  can  be  obtained  over  a  wide  range 
of  excitation  frequencies.  The  excitation  tuning  band¬ 
width  (100  A)  is  determined  by  the  slope  of  the  filter 
absorption  edge.  The  roll-off  on  the  long  wavelength  end 
is  due  to  decreased  optical  coupling  in  the  material  while 
the  roll-off  on  the  short  wavelength  end  is  due  to  pump 
absorption.  The  large  bandwidth  and  narrow  NDFWM  re¬ 
sponse  suggests  that  a  tunable  filter  based  on  this  material 
could  have  over  109  channels.  Alternatively,  such  prop¬ 
erties  may  be  important  to  phase  locking  of  narrow-band 
diode  lasers.  Larger  bandwidths  are  obtained  at  higher 
pump  powers,  due  to  the  decrease  in  trap  lifetime  by  pho¬ 
toionization.  At  this  time,  the  variation  of  the  NDFWM 
response  with  excitation  frequency  is  not  completely 
understood. 

Finally,  for  application  to  either  OPC  or  optical  filter¬ 
ing,  we  need  to  determine  the  angie  dependence  and  thus 
field  of  view  at  low  temperature.  As  discussed  above,  in 
SDG,  the  microcrystallites  are  small  compared  to  an  op¬ 
tical  wavelength  and  surrounded  by  insulator,  hence  we 
would  expect  no  excitation  diffusion.  For  reasons  which 
we  do  not  completely  understand  at  this  point,  we  do  ob¬ 
serve  an  angle  dependence  at  125  K.  Fig.  5  shows  the 
NDFWM  bandwidth  of  the  grating  produced  by  the  for¬ 
ward  pump  and  probe  as  a  function  of  the  forward  pump- 
probe  angle.  For  these  measurements,  the  backward  pump 
polarization  was  orthogonal  to  the  copolarized  forward 
pump  and  probe  in  order  to  eliminate  the  contribution  from 
the  grating  produced  by  the  backward  pump  and  probe. 
For  all  beams  copolarized,  the  angle  dependence  would 
be  weaker  since  the  back  grating  contribution  would  dom¬ 
inate  for  large  angles  with  a  bandwidth  and  efficiency  sim¬ 
ilar  to  the  forward  grating  at  small  angles.  At  low  tem¬ 
peratures,  the  data  suggest  that  even  though  the 
microcrystallites  are  surrounded  by  insulator,  there  is  still 
energy  transfer  taking  place*-  although  the  angle  depen¬ 
dence  is  far  from  the  predictions  of  simple  diffusion.  We 
should  note,  however,  that  even  though  the  bandwidth 
variation  is  a  factor  of  4  over  a  field  of  view  in  the  plane 
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Probe  Detuning  (Hz.) 

Fig.  4.  The  narrow  bandpass  optical  filter  characteristics  of  semiconduc¬ 
tor-doped  glass.  The  x  axis  is  the  pump  probe  detuning.  The  v  axis  is  the 
signal  magnitude.  The  :  axis  is  the  optical  excitation  frequency  (note  the 
pump-probe  detuning  is  very  small  compared  to  the  frequency  of  the 
optical  excitation.)  This  data  were  obtained  at  pump  intensity  below  sat¬ 
uration  of  the  response. 


External  Angle  (Dag.) 

Fig.  5.  The  NDFWM  bandwidth  of  the  forward  pump-probe  grating  as  a 
function  of  angle  between  the  forward  pump  and  probe.  The  solid  line 
is  intended  to  guide  the  eye. 

of  ir  when  the  back  pump  beam  is  orthogonally  polarized, 
this  is  quite  small  compared  to  many  diffusion-limited 
materials.  In  addition,  in  the  case  that  all  beams  are  co¬ 
polarized,  the  reflectivity  and  NDFWM  bandwidth  vari¬ 
ations  are  less  than  a  factor  of  3  over  2t. 

In  summary,  we  have  described  the  large  third-order 


nonlinear  optical  response  in  CdSSE-doped  glass.  The 
optical  phase  conjugation  reflectivities  are  limited  to  about 
0.34  percent  at  123  K  (1  percent  at  3  K)  although  the 
efficiency  is  certainly  large  enough  to  permit  laboratory 
studies  of  real  time  aberration  correction,  for  instance  in 
studies  of  imaging  through  some  forms  of  biological  sys¬ 
tems.  Furthermore,  the  high  efficiencies  are  obtainable 
with  low-power  lasers.  In  addition,  we  have  described  the 
behavior  of  this  material  as  it  might  be  important  to  tun¬ 
able  narrow -band  optical  filters.  This  latter  example  is  in¬ 
teresting  in  that  it  might  provide  a  convenient  material  to 
test  filtering  properties  when  the  probe  beam  is  a  quasi¬ 
white  light  source. 
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ABSTRACT 

A  new  class  of  measurements  on  GaAs  quantum  well  structures  based  on 
frequency  domain  nonlinear  laser  spectroscopy  is  described.  Room 
temperature  measurements  of  the  excitation  relaxation  show  an  interference 
effect  in  the  line  shape  which  is  interpreted  as  a  shift  in  exciton 
frequency.  Low  temperature  measurements  on  the  localized  exciton  provide  an 
excitation  line  shape  which  eliminates  the  effects  of  inhomogeneous 
broadening  and  shows  the  presence  of  spectral  diffusion. 

The  room  temperature  nonlinear  optical  response  near  the  band  edge  in 
GaAs  quantum  wells  has  resulted  in  considerable  interest  in  these  materials 
because  of  their  potential  importance  to  optoelectronic  devices.  Studies 
of  these  materials  have  also  resulted  in  new  understanding  of  the 
fundamental  physics  of  excitons.2'3,4  In  this  paper,  we  describe  a  new 
class  of  measurements  based  on  high  resolution  frequency  domain  nonlinear 
laser  spectroscopy  which  are  distinguished  from  earlier  measurements  because 
of  the  ability  to  report  on  excitation  line  shapes  as  weil  as  the  ability  to 
provide  greater  sensitivity  in  measurements  of  excitation  relaxation.  The 
measurements  of  the  excitation  line  shape  are  significant  because  this 
spectroscopy  eliminates  the  contribution  of  inhomogeneous  broadening  which 
characterizes  low  temperature  absorption  measurements. 

The  basic  experimental  configuration  for  these  experiments  is  based  on 
backward  four  wave  mixing.5  In  these  experiments,  two  optical  beams 
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intersect  each  other  at  a  small  angle.  One  beam  is  designated  the  probe 
beam  with  field  amplitude  f  (u  tic  )  while  the  remaining  beam  is  the  forward 

p  p  p 

pump  beam  with  field  amplitude  The  beams  are  nearly  normal  to 

the  growth  direction  of  the  quantum  well  and  produce  an  interference  pattern 

which  results  in  a  spatial  and  temporal  modulation  of  excitation 

proportional  to  ^  (u  ,1c  )«£  (u  ,ic  ).  A  coherent  signal,  £  (u  ,ic  ),  is 

fee  p  p  p  sss 

produced  by  the  scattering  of  a  back  beam  with  field  amplitude  £  (u  ,ic  ). 

b  b  b 

Spectroscopic  information  is  obtained  by  varying  any  one  frequency  with 
respect  to  the  remaining  two  frequencies.  In  particular,  the  FWMp  response 
is  obtained  by  varying  u  holding  u  =<j  while  the  FWMb  response  is  obtained 

p  f  b 

by  varying  u  ,  holding  u  -u  =constant.  The  FWMp  response  provides  a  line 

b  p  f 

shape  related  to  the  excitation  decay  rate  and  contains  contributions  such 
as  inelastic  scattering  and  carrier  or  excitonic  recombination  in  addition 
to  the  contribution  of  nonradiating  states.  The  FWMb  response  provides  a 
measure  of  the  excitation  line  shape  without  contributions  from 
inhomogeneous  broadening. 5,6 

At  room  temperature,  the  exciton  absorption  line  shape  is  dominated  by 
homogeneous  broadening  due  to  the  rapid  LO  phonon  ionization  rate.  Hence,  a 
study  of  the  FWMb  response  provides  little  additional  information.  However, 
the  nonlinear  response  arises  because  the  electron-hole  plasma  generated 
from  the  ionization  of  the  exciton  results  in  band  filling  and  exchange 
effects  which  lead  to  a  modification  of  the  optical  properties  of  the 
exciton  including  the  oscillator  strength,  dephasing  rate,  and  resonance 
frequency.  Hence,  a  grating  produced  by  the  interference  of  the  forward  pump 
and  probe  has  a  decay  rate  determined  by  the  e-h  recombination  time  and  the 
time  it  takes  carriers  to  spatially  diffuse  across  the  grating.  Figure  la 
shows  the  FWMp  response  at  room  temperature  along  with  a  fit  of  a 
Lorentzian.  The  line  width  is  given  by  j  +  DIAitl2  where  y  is  the  e-h 
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recombination  rate,  D  is  the  ambipolar  diffusion  coefficient,  and  I  Ale  I 

=  |ic  -lc  I.  Figure  lb  shows  the  angie  dependence  of  the  FWMp  line  width  as  a 
f  p 

function  of  angle  and  the  solid  line  is  a  fit  of  the  above  line  width 
expression.  From  the  shape  of  the  curve,  we  can  obtain  the  ambipolar 
diffusion  coefficient  (18  cm2/sec)  and  from  the  y-intercept,  we  can  obtain 
the  e-h  recombination  rate  (5  nsec. ) 

Using  the  method  of  correlated  optical  fields7  to  eliminate  interlaser 

jitter  contributions  to  the  line  shape,  we  examined  the  FWMp  response  near 

q  -u  =0.  Figure  2a  shows  the  line  shape.  There  are  two  important 
p  f 

observations  which  can  be  drawn  from  Fig.  2a.  First,  there  exists  a  strong 
(dominant)  contribution  to  the  room  temperature  nonlinear  response  which  is 
slow  (lOjisec).  Secondly,  the  line  shape  is  characterized  by  a  shape  which 
is  typical  of  interference  effects.8  We  believe  this  line  shape  is  due  to  a 
shift  in  the  exciton  resonance.  Solutions  to  the  optical  Bloch  type 
equations  which  have  been  phenomenologically  modified  to  account  for 
bandfilling  and  exchange  effects  show  that  this  interference  profile  results 
when  the  Bloch  equations  are  further  modified  to  account  for  a  shift  in  the 
exciton  resonance.9  The  curve  in  Fig.  2b  is  the  calculated  FWMp  response, 
showing  the  interference  line  shape.  It  is  important  to  stress  that  the 
solutions  of  the  modified  Bloch  equations  only  show  interference  effects 
when  the  exciton  resonance  is  allowed  to  shift  due  to  the  presence  of  the 
e-h  plasma  or  some  other  laser  induced  mechanism.  We  also  note  that  the 
shape  of  the  FWMp  response  is  dependent  on  the  wavelength  of  the  back  pump. 
Figure  2c  shows  the  FWMp  line  shape  measured  when  the  back  pump  is  tuned  red 
of  the  peak  of  the  heavy  hole  exciton.  (The  forward  pump  and  probe  beams 
remained  tuned  on  the  peak  of  the  heavy  hole  exciton. )  Figure  2d  shows  the 
line  shape  calculated  from  the  modified  Bloch  equations  when  the  back  pump 
is  shifted  red  of  the  exciton  resonance. 
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At  low  temperatures,  the  excitons  are  stable  against  phonon  ionization. 
Hence,  studies  of  relaxation  provide  information  on  the  exciton  dynamics. 
In  addition,  because  the  exciton  is  now  stable,  data  obtained  from  the  FWMb 
response  provide  information  on  the  excitation  line  shape.  At  5  K,  the 
absorption  line  shape  is  inhomogeneously  broadened  due  to  fluctuations  in 
the  well  thickness.  The  HH1  absorption  line  width  in  this  sample  is  3.7 
meV.  The  sample  was  further  characterized  by  performing  luminescence 
measurements  as  well  as  obtaining  a  degenerate  four  wave  mixing  (DFWM) 
spectrum.  (A  DFWM  spectrum  is  obtained  by  using  one  laser  for  all  beams  in 
the  nonlinear  interaction  and  observing  the  change  in  intensity  of  the 
signal  as  the  frequency  of  the  laser  is  tuned  across  the  exciton  resonance. ) 
The  peak  of  the  luminescence  from  the  excitons  is  shifted  2.8  meV  to  the  red 
of  the  HH1  absorption  maximum.  The  red  shifted  luminescence  is  interpreted 
as  being  due  to  the  recombination  of  excitons  which  have  become  localized 
due  to  fluctuations  in  well  thickness.  The  DFWM  measurements  show  a  peak 
which  coincides  exactly  with  the  luminescence  peak,  suggesting  the  nonlinear 
response  is  also  due  to  localized  excitons. 

In  the  first  set  of  measurements,  we  examined  the  FWMp  response  to 
determine  the  excitation  relaxation  dynamics.  In  the  simple  picture  of  the 
exciton  undergoing  a  single  relaxation  due  to  radiative  decay,  the  line 
shape  would  be  a  single  Lorentzian.  Figure  3  shows  the  FWMp  line  width 
obtained  at  5  K.  The  curve  shows  clear  evidence  of  multiple  decay  processes 
rather  than  a  single  decay.  In  particular,  three  decay  channels  are 
observed  in  Fig.  3.  The  fastest  and  slowest  channels  having  relaxation 
times  of  100  psec  and  15  nsec  are  clearly  observed  in  this  figure.  However 
a  third  component  to  the  line  shape  corresponding  to  a  relaxation  time  of 
1.5  nsec  is  observable  upon  expansion  of  the  center  portion  of  the  figure. 
The  1.5  nsec  feature  corresponds  to  the  expected  radiative  aecay  time  lor  a 
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well  of  thickness  of  98A.10  The  origin  of  the  100  psec  feature  is  most 
likely  due  to  decay  of  the  excitation  by  spectral  diffusion  via  phonon 
assisted  tunneling.  Similar  behavior  has  been  reported  in  transient 
pump-probe  measurements. 4  Further  evidence  for  this  is  given  below.  The 
origin  of  the  15  nsec  decay  rate  is  due  possibly  to  the  presence  of  electric 
fields  in  the  material  which  are  known  to  extend  the  life  time  of  the 
exciton. 

FWMb  experiments  were  then  conducted  to  determine  the  line  shape  of  the 

excitation.  Note  that  the  energies  at  which  the  following  FWMb  line  shapes 

were  recorded  are  referenced  with  respect  to  the  position  of  the  peak  of  the 

absorption.  Figure  4a  shows  the  FWMb  line  shape  obtained  3.45  meV  to  the 

low  energy  side  of  the  absorption  maximum.  The  line  shape  is  slightly 

asymmetric  on  the  high  energy  side.  This  asymmetry  becomes  much  clearer  in 

the  FWMb  spectra  recorded  at  the  absorption  line  center  as  seen  in  Fig.  4b. 

The  presence  of  asymmetry  in  the  FWMb  response  is  highly  significant,  since 

it  can  be  shown6  that  if  we  are  measuring  just  the  homogeneous  line  shape, 

the  FWMb  response  should  be  symmetric  about  u  -w  =0,  even  if  the  homogeneous 

line  shape  is  asymmetric.  However,  such  asymmetric  line  shapes  can  be 

observed  in  the  FWMb  response  in  the  presence  of  spectral  diffusion.  Hence 

we  believe  these  measurements  show  directly  that  excitons  are  spectrally 

diffusing,  as  suggested  by  earlier  work.4  Figure  4c  shows  a  comparison  of 

the  measured  temperature  dependence  of  the  line  width  in  4a  to  a  theory 

(3T1 ' 6 

based  on  phonon  assisted  tunnelling.  The  curve  varies  as  e  and  is  in 

general  agreement  with  theoretical  expectations. 11 

Further  evidence  for  spectral  diffusion  is  shown  in  Fig.  5.  Figure  5a 
is  an  extended  scan  of  the  data  shown  in  Fig.  4b  at  5  K.  Figure  5b  shows 
the  same  data  taken  at  15  K  illustrating  the  spectral  diffusion  of  excitons 
to  lower  energies,  the  peak  of  the  response  occur ing  at  the  same  energy  as 
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the  peak  of  the  luminescence.  As  in  the  case  of  the  luminescence  peak,  the 
wide  resonance  in  the  FWMb  spectra  shown  in  figure  5b  represents  the  energy 
distribution  of  the  localized  excitons.  Theoretical  results  of  the  FWMb 
line  shape  obtained  from  optical  Bloch  type  equations,  modified  to  include 
spectral  diffusion,  show  similar  behavior.  It  is  interesting  to  note  that 
while  we  obtain  clear  evidence  for  energy  diffusion,  there  is  no  evidence  in 
the  data  to  suggest  that  the  excitons  are  spatially  diffusing. 

In  summary,  we  have  used  high  resolution  frequency  domain  nonlinear 
spectroscopy  to  study  relaxation  and  excitation  line  shapes  in  GaAs  multiple 
quantum  wells.  At  room  temperature,  we  have  measured  the  ambipolar 
diffusion  coefficient  and  the  free  carrier  recombination  time.  In  addition, 
we  have  observed  an  interference  effect  in  the  FWMp  line  shape  which  we 
believe  is  due  to  a  shift  in  the  exciton  resonance.  At  low  temperatures,  we 
have  identified  contributions  to  the  relaxation  of  the  exciton  population 
due  to  spectral  diffusion  and  exciton  recombination.  Using  FWMb 
spectroscopy,  we  have  obtained  measurements  of  the  line  shape  of  the 
excitation.  The  observation  of  asymmetric  FWMb  line  shapes  and  the  measured 
temperature  dependence  of  the  line  width  provide  evidence  that  the  excitons 
are  spectrally  diffusing  via  phonon  assisted  tunnelling. 

This  work  is  supported  by  the  AFOSR,  U.S.  ARO  and  the  U.S.  ARO  URI. 
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Figure  Captions 


Figure  1. 
Figure  2. 

Figure  3. 
Figure  4. 

Figure  5. 


Figure  la  shows  the  room  temperature  FWMp  line  shape.  Figure 
2a  shows  the  angle  dependence  of  the  FWMp  line  width. 

The  FWMp  line  shape  using  the  technique  of  correlated  fields. 
Figure  2a  shows  the  interference  profile  when  the  back  beam 
wavelength  is  on  the  exciton  peak.  Figure  2c  shows  the  line 
shape  when  the  back  beam  is  tuned  red  of  the  exciton  peak. 
Figures  2b  and  2d  show  the  calculated  line  shapes  using 
modified  Bloch  equations. 

The  FWMp  line  shape  at  5  K  showing  the  presence  of  multiple 
relaxation  mechanisms. 

The  FWMb  line  shape  at  5  K.  Figure  4a  shows  the  line  shape 
red  of  the  absorption  maximum.  Figure  4b  shows  the  line 
shape  measured  at  the  peak  of  the  exciton  absorption.  Figure 
4c  shows  the  temperature  dependence  of  the  line  width  in 
figure  4a. 

FWMb  response  shown  for  2  temperatures .  Figure  5a  shows  the 
response  at  5  K.  Figure  5b  shows  the  response  at  15  K 
illustrating  the  spectral  diffusion  of  the  excitation  to 
lower  energies. 
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ABSTRACT 

This  paper  describes  cw  frequency  domain  nonlinear  laser 
spectroscopy  results  obtained  in  the  study  of  relaxation  in 
glasses  doped  with  microcrystallites  of  CdSSe.  Measurements 
are  made  as  a  function  of  temperature  and  optical  intensity 
using  low  intensity  cw  optical  excitation.  The  results  are 
interpreted  based  on  the  assumption  that  the  dynamics  of  the 
nonlinear  response  in  this  material  are  controlled  by  traps. 
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Applications  to  nonlinear  optical  problems  and  the 
possibility  of  learning  about  the  physics  of  quantum  size 
effects  has  lead  to  considerable  recent  interest  in  glasses 
doped  with  semiconductor  microcrystallites.  These  glasses  are 
doped  with  CdS^  ^Se^  and  are  commercially  available  as 
sharp-cut  color  glass  filters1. 

In  the  commercial  material,  the  microcrystallites  have  an 
average  diameter  typically  of  10 0A  with  a  FWHM  distribution  of 
50A,  as  determined  by  electron  microscopy  (obtained  for  Corning 
material2) ,  resulting  in  a  relatively  featureless  room 
temperature  absorption  edge.  At  low  temperature,  these 
materials  have  been  reported  to  exhibit  structure  in  the 
luminescence  spectrum  attributed  to  quantum  confinement3  and 
more  recent  results  have  shown  the  existence  of  structure  in 
the  modulated  absorption  spectrum4.  More  careful  material 
processing  methods  have  produced  material  with  a  uniform  size 
distribution  showing  clear  structure  in  the  room  temperature 
linear  absorption  spectrum  which  have  been  attributed  to 
quantum  size  effects5.  One  group6  has  produced  structures  as 
small  as  12A  in  CdS,  CdSe,  and  CuCl.  Another  group  has  made 
glasses  doped  with  microcrystallites  where  the  size 
distribution  is  controlled  by  a  heat  treatment  process  and 
microcrystallites  as  small  as  25A  have  been  observed.2 

Because  of  the  relative  ease  of  fabrication,  the  large 
fast  nonlinear  optical  response,  and  the  fact  that  the  host 
material  is  a  glass,  the  applications  of  these  materials  to 
optical  switches  and  integrated  optics  have  been  studied  by 
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numerous  groups.  In  addition,  there  have  been  several  studies 
discussing  the  physics  of  these  materials.  Very  early,  the 
material  was  used  to  demonstrate  a  Q-switch  in  a  laser7. 
Optical  bistability  was  demonstrated  by  McCall  and  Gibbs8  and 
Gibbs  et  al.9,  however  the  slow  response  suggested  the 
nonlinearity  was  thermal  in  origin10.  Using  a  Q-switched 
laser,  Jain  and  Lind11  demonstrated  optical  phase  conjugation 
and  showed  that  a  large  nonlinear  optical  susceptibility  was 
observable  (*  *10  esu)  using  transient  excitation.  They 
determined  that  the  relaxation  time  was  faster  than  8  nsec 
(determined  as  an  upper  limit.)  The  model  for  the  nonlinear 
response  proposed  by  Jain  and  Lind  was  based  on  the  production 
of  a  short-lived  electron-hole  plasma.  This  work  was  followed 
by  numerous  other  experimental12-22  and  theoretical23-25 
reports  of  the  nonlinear  response  using  pulsed  excitation. 
Included  in  these  observations  have  been  several  measurements 
of  the  dynamical  behavior26-29  of  key  importance  to  high  speed 
applications.  These  measurements  show  carrier  relaxation  times 
are  in  the  picosecond  time  domain.  The  range  of  fast 
relaxation  times  reported  in  the  literature  has  been  attributed 
to  defects  which  decrease  the  relaxation  time  and  which  are 
induced  by  high  power  optical  irradiation  (photo-darkening)21. 
More  recently,  using  the  methods  of  ultrafast  laser 
spectroscopy,  measurements  of  transient  absorption  have  shown 
relaxation  processes  occurring  on  the  time  scale  of  a  few 
hundred  femtoseconds30,31.  Measurements  on  CdSSe  glasses  in 
which  quantum  confinement  effects  dominate  have  demonstrated 
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the  importance  of  phonon  broadening  and  have  shown  evidence  for 
spectral  hole  burning  at  low  temperature.  Other  measurements 
of  hole  burning  in  small  (35-55A)  microcrystallites  of  CdSe 
clusters  suspended  in  a  polystyrene  film  have  resulted  in 
measurements  of  the  contributions  of  homogeneous  and 
inhomogeneous  broadening  to  the  electronic  absorption 
spectra.33  A  discussion  of  the  recent  applications  is  given  in 
reference  34. 

In  this  paper,  we  examine  the  physical  behavior  of  a 
slower  but  much  larger  nonthermal  contribution  to  the  nonlinear 
optical  response  obtained  under  cw  excitation.  The  first 
indication  of  the  possibility  of  a  large  cw  nonlinear  response 
was  provided  by  the  observation  of  low  frequency  dynamical 
behavior  leading  to  subharmonic  generation  and  near  chaotic 
behavior  in  SDG  observed  at  low  temperature  (14 OK)  by  Zheludev 
et  al.35  This  work  was  followed  by  measurements  in  our 
laboratory  of  the  room  temperature  cw  nonlinear  response  using 
DFWM  which  showed  a  nonlinear  response  much  larger  but  slower 
than  the  response  obtained  based  on  pulsed  excitation.  The 
results  showed  *^*10  7  esu.  Using  four-wave  mixing 
spectroscopy  methods,  a  narrow  resonance  (4.4kHz)  was  found  in 
the  nonlinear  susceptibility  corresponding  to  a  72*xsec 
excitation  relaxation  time.  In  that  work,  we  attributed  the 
origin  of  the  nonlinear  dynamics  to  contributions  from  traps,  a 
result  suggested  by  earlier  measurements3,13  and  confirmed  in 
more  recent  work29’37. 

This  paper  extends  the  earlier  measurements  to  show  the 
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temperature  and  optical  intensity  dependence  of  the  nonlinear 
response.  These  results  are  discussed  based  on  the  assumption 
that  the  dynamics  of  the  nonlinear  response  are  controlled  by 
traps.  As  a  result  of  this  work,  we  have  demonstrated34  that  a 
large  third  order  susceptibility  (*<3)  *  2xl0~4  esu)  can  be 
obtained  at  12 5K.  In  addition,  the  current  work  provides  a 
measure  of  an  effective  thermal  activation  energy  and  a  model 
is  discussed  which  suggests  that  the  saturation  of  the 
nonlinear  response  reported  in  our  earlier  work  is  due  to 
saturation  of  the  available  traps.  The  large  nonlinear 
response  has  been  used  to  demonstrate  efficient  optical  phase 
conjugation,  and  the  narrow  bandwidth  and  near  angle 
independent  response  has  been  used  to  demonstrate  a  tunable 
optical  filter34. 

The  experimental  configuration  for  these  measurements  is 
based  on  backward  four-wave  mixing  using  a  frequency  stabilized 
tunable  cw  dye  laser.  Two  counter  propagating  pump  beams  with 
fields  £  (cj)  and  ?  (u)  interact  with  a  probe  beam  ^  (w=<4>+5)  in 

f  b  r  p  p 

the  sample  which  is  contained  in  a  variable  temperature 


cryostat.  The  angle  between  the  forward  pump  and  probe  beams 
is  ©.  The  three  beams  interact  through  the  nonlinear  response 


to  produce  a  signal  field  (2  )  arising  from  a  polarization 
proportional  to  *<3)  (tot')t.  For  these  measurements, 

f  p  b 

(no  signal  is  obtained  for  for  low  power  cw 

excitation) ,  and  the  signal  arises  from  the  coherent  scattering 


of  £  from  the  traveling  wave  grating  of  excitation  produced  by 

b 

H  .£•.  In  the  limit  where  the  dephasing  rate  is  large  compared 
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to  the  excitation  relaxation  rates,  a  measurement  of  |E  I2  as  a 

■ 

function  of  S  for  fixed  u>  (the  FWMp  response,)  gives  a  line 
shape  related  to  relaxation  rates  associated  with  the 
excitation  grating  produced  by  the  forward  pump  and  probe.  (If 
the  dephasing  rate  is  not  large,  the  FWMp  response  may  include 
interference  effects  arising  from  the  induced  dipole.38)  A 
simple  physical  understanding  of  the  FWMp  response  can  be 
obtained  by  considering  the  response  in  the  limit  of  small  and 
large  5.  For  6  small  compared  to  the  excitation  relaxation 
rate,  the  spatial  modulation  of  the  excitation  is  in  phase  with 
the  intensity  pattern  formed  by  EE*,  and  a  strong  signal  is 

f  p 

produced  by  the  scattering  of  the  baclcward  pump  beam  from  the 
grating.  However,  for  6  large  compared  to  the  excitation 
relaxation  rates,  the  spatial  modulation  of  the  excitation  is 
washed  out  and  there  is  no  grating  to  produce  a  coherent 
signal . 39 

In  this  system,  we  assume  that  the  nonlinear  response  is 
due  to  phase  space  filling  effects  due  to  the  free  carriers 
produced  by  the  forward  pump  and  probe  which  result  in  a 
modulation  of  the  transition  oscillator  strength40.  As 
indicated  above,  because  of  the  slow  response,  it  is  believed 
that  the  dynamics  of  the  nonlinear  response  are  dominated  by 
the  presence  of  traps.41  If  an  electron  (or  hole)  decays  to  a 
long  lived  trap,  the  remaining  hole  (or  electron)  can  still 
contribute  to  the  bandfilling.  Thus  a  component  of  the 
bandfilling  relaxation  rate  will  be  determined  by  the  effective 
decay  time  of  the  carriers  from  the  traps.  The  line  width  of  a 
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FWMp  measurement  is  thus  related  to  the  effective  trap 
relaxation  rate. 

A  detailed  study  of  the  FWMp  line  width  was  performed  using 
the  experimental  setup  described  in  reference  36.  Since  the 
line  widths  were  much  smaller  than  the  interlaser  jitter 
associated  with  a  measurement  based  on  one  stabilized  dye  laser 
providing  the  pump  beams  and  a  second  dye  laser  providing  the 
probe  beam,  we  used  the  method  of  correlated  optical  fields 
using  two  acoustooptic  modulators  driven  by  two  phase  locked 
frequency  synthesizers.  In  this  approach,  the  first  order 
Bragg  scattered  beam  from  the  first  modulator  provided  the  pump 
beams  while  the  second  AO  cell  provided  the  probe  beam.  The 
first  AO  modulator  was  operated  at  a  fixed  drive  frequency  of 
40  MHz.  The  second  driver  could  be  varied  about  the  40  MHz 
drive  frequency.  In  this  way,  we  could  vary  S  over  the  range 
needed  for  the  measurement  to  obtain  line  shapes  without 
contributions  from  laser  jitter. 

Figure  1  shows  the  FWMp  line  width  obtained  at  12 5K  in 
RG630.  The  data  shows  that  the  profile  is  not  a  pure 
Lorentzian,  suggesting  that  if  trap  dynamics  are  responsible 
for  the  behavior,  then  most  likely,  the  response  is  an  average 
over  a  distribution  of  traps,  a  possibility  supported  by 
additional  data  below.  The  inset  on  the  left  shows  the  FWMp 
profile  as  a  function  of  excitation  wavelength  around  the 
filter  absorption  edge.  There  is  a  dependence  on  excitation 
wavelength  with  the  line  width  varying  from  214  Hz  to  429  Hz  at 
this  temperature  as  the  wavelength  is  varied  from  600  nm  to  614 
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run. 


The  inset  on  the  right  shows  this  dependence  with  a  linear 
fit  of  the  data.  In  addition,  the  strength  of  the  nonlinear 
response  also  varies.  On  the  blue  side  of  the  data,  the 
incident  beams  are  strongly  absorbed  by  the  material,  while  on 
the  red  side,  there  is  a  decrease  in  the  optical  coupling.  For 
the  remaining  data  below,  all  measurement  were  made  at  the  peak 
of  the  response. 

The  FWMp  response  was  studied  as  a  function  of 

temperature.  Figure  2  shows  the  log  of  the  FWMp  line  width  as 

a  function  of  1/T.  All  line  widths  were  observed  to  have  a 

linear  dependence  on  optical  intensity.  The  data  used  for  Fig. 

2  corresponds  to  the  zero  intensity  intercept.  (Further 

discussion  is  given  below  regarding  intensity  dependent 

effects.)  The  general  form  of  the  data  in  Fig.  2  suggests  a 

functional  form  for  the  effective  trap  decay  rate  given  by: 
-(<?  /kT) 

7(0)  +  re  where  7(0)  is  the  zero  temperature  decay 

rate,  and  rexp(-&b/kT)  is  the  probability  per  unit  time  for  the 
electron  or  hole  to  escape  from  a  thermally  activated  trap  (6^ 
is  the  activation  energy)  .  The  solid  line  shown  in  Fig.  2  is 
the  best  fit  of  the  functional  form  for  the  effective  trap 
decay  rate  given  above.  The  asymptotic  value  of  the  curve 
gives  the  zero  temperature  decay  rate.  The  slope  of  the  curve 
gives  the  activation  energy,  determined  to  be  137  meV.  The 
decay  parameter,  r,  is  C.7  MHz  and  the  zero  temperature  decay 
rate,  7(0)  is  43  Hz. 

It  is  most  likely  that  at  these  powers,  the  excitation  is 
in  near  equilibrium  with  the  crystal  lattice  imbedded  in  the 
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glass  material  and  the  actual  details  of  decay  are  most  likely 
far  more  complex  than  indicated  by  such  a  simple  description. 
Experiments  have  been  reported  in  colloidal  suspensions  of  CdS 
clusters42  where  luminescence  decay  is  interpreted  based  on  a 
theory  for  multiphonon  assisted  excitation  transfer  .  For  our 
measurements,  we  are  unable  to  distinguish  between  these 
different  decay  mechanisms. 

We  also  observe  that  the  FWMp  line  width  is  intensity 
dependent  as  shown  in  Fig.  3,  showing  a  linear  dependence  on 
intensity  at  low  powers.  Additionally,  as  reported  in 
reference  34,  we  observed  a  roll  off  in  the  nonlinear  response 
at  high  pump  intensities.  (At  125  K,  saturation  of  the 
nonlinear  response  occurred  around  10W/cm2)  .  A  qualitative 
understanding  of  the  saturation  behavior  reported  in  reference 
34  and  the  results  obtained  in  the  FWMp  measurements  can  be 
obtained  by  considering  a  simple  model  of  the  nonlinear 
response  based  on  the  energy  level  diagram  shown  in  figure  4. 
In  this  model,  optical  radiation  induces  a  transition  from  the 
valence  band  to  the  conduction  band,  producing  an  electron-hole 
(e-h)  pair.  The  electron  in  the  conduction  band  can  decay  back 
to  the  valence  band  with  decay  rate  7  cy  or  to  a  trap  with  decay 
rate  KNt  in  the  limit  that  the  traps  are  unsaturated  where  Nt 
is  the  total  trap  density.  (For  this  discussion,  we  consider 
only  electron  traps.)  The  trapped  electron  can  then  be 
thermally  activated  from  the  trap  and  return  to  the  conduction 
band  with  rate  7  or  recombine  with  a  hole  in  the  valence  band 

tc 

with  rate  7 
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In  the  absence  of  traps,  the  FWMp  line  shape  would  have  a 
width  determined  by  the  carrier  recombination  rate  (decay  of 
the  induced  grating  by  spatial  diffusion  of  the  carriers  is  not 
included  since  the  size  of  microcrystallites  is  small  compared 
to  the  grating  spacing  and  they  are  surrounded  by  insulating 
glass.)  However,  if  electrons  (or  holes)  can  be  trapped,  the 
remaining  electrons  and  holes  can  still  contribute  to 
bandfilling.  The  time  constant  for  the  decay  of  this  system 
would  have  a  fast  component  due  to  free  carrier  recombination 
and  a  slow  component  determined  by  the  rate  of  electron  (or 
hole)  escape  from  the  trap. 

In  the  limit  that  the  dynamics  of  the  nonlinear  response 
are  determined  by  trap  kinetics,  this  simplified  model  can  be 
used  to  obtain  a  more  quantitative  understanding  by  considering 
the  rate  equations  for  the  electron  density,  n,  and  trap 
density,  n^: 

n  =  »  -  7cyn  -  K(Nt-nt)n  +  rt(.nt  (1) 
n  =  K(N  -n  ) n  -  7  n  -  7  n  (2) 

t  '  t  t'  tc  t  rec  t  ' 

where  #  represent  the  optical  generation  rate  of  e-h  pairs. 
The  possibility  of  trap  saturation  is  included  through  the  term 
K(Nt-nt) .  The  nonlinear  response  of  this  material  is  the 
result  of  bandfilling  so  that  we  take  the  nonlinear 
susceptibility  to  be  proportional  to  the  density  of  free 
electrons  and  holes40.  However,  since  experimentally  we  detect 
only  a  slow  component  of  the  nonlinear  response  due  to  the 
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presence  of  holes  that  persist  because  of  the  presence  of 

electron  traps,  we  take  the  dominant  term  in  the  nonlinear 

susceptibility  to  be  proportional  to  the  filled  trap  density. 

The  above  set  of  rate  equations  can  be  solved  in  the  presence 

of  strong  pumps.  However,  the  essential  features  of  importance 

to  spectroscopy  can  be  seen  from  the  solution  obtained  in  the 

limit  of  perturbation  theory  and  in  the  absence  of  trap 

saturation.  (Trap  saturation  effects  are  discussed  below.) 

a  .  » 

Taking  #  =  -j—I  where  I=^/[jEicos(kix-cjit+^i)  ]2\  ,  we  obtain 

i  t 

an  expression  for  the  nonlinear  polarization  for  the  FWMp 
response  for  the  interaction  described  above  in  the  limit 
S=u-u  «r  ,  KN  : 

f  p  cv'  t 


„l  *  i(k  -k +k  )x-i(6+u  )t 
P  =£6  &  8  e  f  p  b  b  - 

f  p  b 


15+7-  +7, 

rec  tc 


r  \ 

7  +KN  ) 

CV  t  ' 


+  c.c 


i  (k  x-<*>  t) 

where  £  is  a  constant  and  E  cos(kx-ut+d  e  1  1  +c.c. 
Thus,  in  the  case  of  a  single  trap  (i.e.,  no  distribution  of 
trap  lifetimes)  the  line  width  at  zero  intensity  measured  in  a 
FWMp  experiment  is  given  by  an  effective  trap  lifetime  given  by 

r«ff  ~  *rec+ytc(r  +KN  )  * 

v  C  V  t/ 

When  the  effects  of  trap  saturation  are  included,  the 
model  predicts  a  roll  off  in  the  nonlinear  response  at  high 
pump  intensities,  as  well  as  a  linear  increase  in  the  line 
width  with  pump  intensity  in  the  low  pump  intensity  limit.  In 


11 


particular,  the  above  equations  can  be  solved  in  the  presence 
of  strong  pumps  (standing  wave  effects  are  ignored  since  the 
pump  polarizations  are  orthogonal.)  in  the  limit  that  the 
dominate  decay  channel  from  the  trap  is  7tc,the  above 
polarization  describes  the  intensity  dependent  nonlinear 
response  if  we  replace  N  with  N  -  n  and  replace  r 

t  t  t(dc)  tc 

with  7  +  Kn#J  „  where 

tc  (dc) 


n 


t  (dc) 


ytc+Kn(dc, 


(4) 


(dc) 


<X  h  1 

ab  s  o 
h(J 


(5) 


and  Jo  is  the  pump  intensity.  These  equations  do  show  the 
linear  dependence  of  the  line  width  on  intensity,  and  the  roll 
over  at  high  intensity,  however  using  the  single  trap  level 
model  presented  above,  it  is  not  possible  to  reconcile  the  low 
intensity  line  width  broadening  rate  with  the  observed 
saturation  intensity.  Estimation  of  the  saturation  intensity 
from  the  low  intensity  line  width  broadening  rate  predicts  a 
value  of  saturation  intensity  much  lower  than  what  was  observed 
experimentally.  Currently,  we  believe  the  explanation  for  this 
difference  is  that  this  material  is  characterized  by  a  range  of 
trap  lifetimes  with  correspondingly  different  concentrations 
and  positions  within  the  band-gap.  As  the  intensity  of  the 
pumps  is  increased,  the  traps  begin  to  fill  and  different 
groups  of  traps  in  the  distribution  then  contribute  to  the 
response.  Thus  at  low  pump  intensities,  the  response  most 
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likely  results  from  one  group  of  traps  and  the  values  of  trap 
relaxation  time  and  activation  energy  thus  pertain  only  to  this 
group  of  traps.  However  at  the  higher  intensities,  different 
parts  of  the  trap  distribution  contribute  to  the  signal 
resulting  in  an  effective  saturation  intensity  that  will  differ 
from  the  prediction  based  on  assuming  a  single  trap  and  using 
the  measured  low  intensity  line  width  broadening  rate.  Such  an 
intexrpretation  is  supported  by  recent  measurements  by  Tomita  et 
al37  which  have  shown  nonexponential  decays  in  luminescence 
lifetimes  with  time  constants  ranging  from  2usec  up  to 
lifetimes  beyond  20jisec  at  liquid  nitrogen  temperature.  These 
multiple  exponential  decays  are  consistent  with  our 
non-Lorentzian  line  shapes. 

In  summary,  these  measurements  provide  further  insight 
into  the  nature  of  the  nonlinear  response  obtained  in  CdSSe 
doped  glass  using  cw  excitation.  Using  frequency  domain 
nonlinear  laser  spectroscopy  methods,  we  have  determined  the 
effects  of  temperature  and  pump  intensity  on  the  dynamical 
behavior  of  the  nonlinear  response.  There  is  strong  evidence 
that  traps  are  responsible  for  the  dynamics  of  the  cw  nonlinear 
response  and  that  trap  saturation  is  responsible  for  the 
observed  roll  off  of  the  response. 

This  work  was  supported  by  AFOSR,  Grant  No.  85-0280.  In 
addition,  M.D.  Webb  was  supported  by  an  Army  Research  Office 
Fellowship  from  the  University  Research  Initiative  program. 
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FIGURE  CAPTIONS 


Figure  1.  The  FWMp  spectrum  obtained  using  the  method  of 
correlated  optical  fields.  The  inset  on  the  left  shows  the 
FWMp  line  shape  obtained  for  different  excitation  wavelengths 
at  a  constant  temperature.  The  inset  on  the  right  shows  the 
dependence  of  line  width  on  wavelength. 

Figure  2.  The  FWMp  line  width  as  a  function  of  inverse 
temperature . 

Figure  3.  The  FWMp  line  width  as  a  function  of  pump 
intensity. 

Figure  4 .  Energy  level  diagram  describing  the  optical 
interaction  believed  responsible  for  the  nonlinear  response  in 
CdSSe  doped  glass. 
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Abstract 

We  report  high  resolution  nonlinear  laser  spectroscopy  measurements  of  relaxation 
of  lowest  energy  heavy  hole  excitons  in  GaAs/AlGaAs  multiple  quantum  well  structures. 
We  show  that  excitons  below  the  absorption  line  center  are  spatially  localized,  and  migrate 
among  localized  sites  with  a  time  scale  of  order  100  ps.  The  measurement  also  gives  the 
resultant  quasi-equilibrium  energy  distribution  of  the  scattered  excitons,  and  based  on  the 
temperatur  dependence  of  the  migration  rate  confirms  the  theoretical  model  for  the  phonon 
assisted  migration. 


1 


Optical  resonant  excitation  of  excitons  with  nearly  monochromatic  light  of  energy  E 
leads  to  an  optically  induced  polarization  (coherence)  and  a  population  of  excitons  within 
AE  of  E  where  AE  ~  fin,  (rh  is  the  homogeneous  width  of  the  exciton).  The  decay  of  this 
excitation  must  be  characterized  by  decay  of  the  polarization  (often  called  dephasing)  as 
well  as  decay  of  the  population  about  energy  E.  In  a  perfect  quantum  well ,  quasi  two 
dimensional  excitons  are  described  by  a  Bloch  type  of  wave  function  and  are  free  to  move 
in  the  well  plane.  At  low  excitation  density,  decay  of  the  excitation  is  then  expected  to  be 
predominantly  due  to  exciton-phonon  scattering  along  with  exciton  recombination.  In 
practice,  however,  the  problem  becomes  more  complicated  due  to  nonideal  growth 
processes.  Recent  transport  and  chemical  lattice  imaging  measurements1^  have  shown  the 
interface  of  GaAs/AlGaAs  multiple  quantum  well  (MQW)  samples  exhibits  island-like 
structures  with  a  height  of  one  monolayer  and  a  lateral  size  of  order  50  A.  For  an  exciton 
confined  to  a  thin  GaAs  layer  in  a  MQW,  its  energy  depends  strongly  on  the  well 
thickness.  For  example,  well  width  fluctuations  of  one  monolayer  in  a  100  A 
GaAs/Alo.3Gao.7As  MQW  can  result  in  a  change  of  exciton  energy  on  the  order  of  several 
meV.  Therefore,  at  low  temperature  in  the  low  energy  region  of  the  heavy  hole  exciton 
(HH1)  absorption  spectrum,  excitons  can  be  localized  in  the  well  with  an  energy 
determined  by  the  local  environment-*  leading  to  inhomogeneous  broadening  of  the  linear 
absorption  spectrum.  Excitons  in  the  high  energy  region  may  still  be  delocalized^.  These 
excitons  are  expected  to  experience  additional  dephasing  due  to  elastic  scattering  from 
potential  fluctuations  in  addition  to  the  decay  due  to  exciton  phonon  scattering  and  exciton 
recombination. 

Localized  excitons  are  in  local  minimum  in  energy,  and  at  very  low  temperature, 
decay  of  the  localized  exciton  is  dominated  by  migration  between  localization  sites.  The 
migration  is  accompanied  by  absorption  or  emission  of  acoustic  phonons  to  compensate  for 
the  energy  difference.  Indeed,  phonon  assisted  migration  was  proposed  to  explain  the  slow 
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and  non-exponential  energy  relaxation  observed  in  time  resolved  luminescence 
measurements  in  a  GaAs  MQW4.  Furthermore,  the  migration  is  primarily  due  to  the  inersite 
dipole-dipole  interaction  when  the  intersite  distance  is  much  greater  than  the  exciton  Bohr 
radius  and  localization  length5.  Furthermore,  the  migration  is  primarily  due  to  the  inersite 
dipole-dipole  interaction  when  the  intersite  distance  is  much  greater  than  the  exciton  Bohr 
radius  and  localization  length5.  Migration  due  to  the  overlap  of  the  exciton  wave  functions 
between  different  sites  also  becomes  important  when  the  intersite  distance  is  comparable 
with  the  exciton  Bohr  radius  or  the  localization  length.  It  is  estimated5  that  the  typical 
magnitude  of  participating  phonon  wave  vectors  is  within  a  few  times  of  the  inverse  of  tne 
localization  length,  which  indicates  the  energy  of  participating  phonons  is  on  the  order  of 
0.01  to  0. 1  meV.  At  higher  temperature  (>10K),  thermal  activation  of  localized  excitons  to 
delocalized  states  becomes  important  This  process  is  associated  with  phonon  absorption, 
and  has  been  observed  in  GaAs  MQW  structures  using  resonant  Rayleigh  scattering3.  The 
estimated  activation  energy  indicates  that  the  onset  for  the  delocalized  exciton  is  near  the 
absorption  line  center. 

A  distinctive  signature  for  phonon  assisted  migration  is  the  temperature  dependence. 
Recent  work  by  Takagahara6  has  shown  the  temperature  dependence  of  the  migration  rate 
in  MQW  structures  has  a  form  given  by  expfBT01),  and  has  been  observed  in  transient  hole 
burning  experiments  in  an  InGaAs/InP  MQW7.  In  this  expression,  B  is  positive  and 
independent  of  temperature  but  is  expected  to  increase  with  the  exciton  energy,  and  a  is 
estimated  to  be  between  1.6  and  1.7.  The  predicted  temperature  dependence  is  quite 
different  from  that  of  variable  range  hopping  used  by  Mott  to  interpret  electronic  conduction 
in  the  localized  regime8,  and  is  attributed  to  the  role  of  the  long-range  dipole-dipole 
interaction  involved  in  the  migration  of  the  localized  exciton5. 

In  this  paper,  we  present  results  of  high  resolution  nonlinear  laser  spectroscopy  of 
lowest  energy  heavy  hole  excitons  in  a  GaAs/AlGaAs  MQW  at  temperatures  between  2.5  K 
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and  15  K.  Using  frequency  domain  four  wave  mixing,  we  are  able  to  obtain  relaxation 
rates  for  the  exciton  population  at  a  given  energy  E,  and  to  measure  directly  the  steady  state 
exciton  redistribution  under  narrow  band  CW  excitation.  The  measurement  shows  excitons 
below  the  absorption  line  center  are  localized,  and  migrate  among  localization  sites  in  a  time 
scale  of  order  100  ps.  The  predicated  temperature  dependence  of  the  migration  rate  is  also 
confirmed. 

The  experimental  configuration  is  discussed  elsewhere9,  and  is  based  on  the  use  of 
two  frequency  stabilized  tunable  CW  dye  lasers.  Three  incident  beams  Ef(C0f,  kf),  Ebfcot,, 
kb)  and  Ep(C0p,  kp)  (f,  b,  and  p  stand  for  forward,  backward  and  probe  respectively)  with 
EfllEpJJEb  interact  in  the  sample  through  the  third  order  susceptibility.  A  signal  beam 
Es(cos,  ks)  with  frequency  cos=(DffCOb-cop  is  generated  and  is  proportional  to  x^EfEp*Eb. 
Physically,  EpEp*  results  in  a  modulation  of  the  exciton  population10.  The  modulation 
further  modifies  the  optical  response  of  the  sample  through  exciton  phase  space  filling  and 
exchange  effects11.  The  signal  arises  from  coherent  scattering  of  the  backward  beam  from 
the  modulation.  Because  of  the  resonant  enhancement  which  occurs  when  the  emission 
dipole  is  resonant  with  the  signal  beam12,  only  excitons  that  are  resonant  with  to^tOf+GV 
(Op  contribute  significantly  to  the  nonlinear  response. 

In  general,  spectroscopic  information  related  to  the  energy  level  structure  and 
relaxation  of  the  system  is  obtained  by  measuring  the  nonlinear  response  as  a  function  of 
the  frequency  of  any  of  the  three  input  beams.  We  further  designate  the  line  shape  obtained 
as  FWMi  where  i  stands  for  f,  b,  or  p  depending  on  which  frequency  is  tuned.  The 
measurement  can  be  understood  in  the  following  simplified  picture  with  a  more  rigorous 
analysis  based  on  the  effective  optical  Bloch  equations  to  be  discussed  elsewhere13. 

In  the  FWMb  measurement,  Ef  and  Ep  (|o>f  -Oipl  «  Hi)  create  an  exciton  population 
in  a  narrow  spectral  hole  with  a  width  AE=firh  within  the  inhomogeneous  width.  Scanning 


cob  w»*.lw  keeping  cof  and  cop  fixed  effectively  changes  0)$.  The  resultant  line  shape  reflects 
a  direct  measure  of  the  steady  state  spectral  redistribution  of  excitons  that  have  a  nonzero 
emission  dipole  moment.  Hence,  in  the  case  of  localized  excitons,  scattering  to  other 
energy  states  can  then  be  directly  measured  in  the  FWMb  response.  However,  while 
delocalized  excitons  can  also  be  scattered  to  other  energies  by  inelastic  processes  such  as 
exciton-phonon  interactions,  delocalized  excitons  in  these  states  have  nonzero  momentum, 
and  as  a  result,  a  zero  emission  dipole  moment  In  this  case,  the  FWMb  line  shape  simply 
provides  a  measure  of  the  exciton  homogeneous  line  shape. 

In  the  FWMf  and  FWMp  measurement  detuning  the  probe  or  forward  pump 
generates  a  traveling  wave  modulation  oscillating  in  time  at  a  fixed  point  in  space  with  the 
detuning  frequency  8=c0f-0)p.  The  nonlinear  response  as  a  function  of  8  then  measures  the 
decay  rate  of  the  modulation  formed  by  excitons  that  are  resonant  with  the  signal  beam.  The 
rate  obtained  includes  contributions  from  exciton  recombination  as  well  as  scattering  of 
excitons  from  energy  E  to  E'  where  IE-E’I  >  IV  Spatial  diffusion  of  the  exciton  also 
contribute  to  the  decay  of  the  modulation  and  manifests  itself  as  the  dependence  of  the 
decay  rate  on  the  grating  period14.  Note  that  tuning  (Of  or  cop  will  also  change  (o§  since 
ci)s=0)ffci)b-0)p.  Complications  may  arise  when  the  exciton  energy  relaxation  rate  is 
comparable  with  Tf,.  However,  such  complications  are  reduced  by  choosing  to  tune  0)f 
instead  of  tOp1^  15,  and  introduce  only  an  overall  correction  factor  (close  to  1)  for  decay 
rates  measured  by  the  FWMf  response.  In  the  limit  that  Th  is  much  larger  than  the  energy 
relaxation  rate,  tuning  G)f  or  cop  gives  the  same  line  shape. 

The  MQW  sample  used  for  our  measurement  consists  of  65  periods  of  96  A  GaAs 
wells  and  98  A  Alo.3Gao.7As  barriers  grown  by  molecular- beam-epitaxy.  The  sample  is 
mounted  on  a  sapphire  disk  (c  axis  normal)  with  the  substrate  removed.  AH  the  experiments 
are  carried  out  on  the  HH1  exciton  with  an  exciton  density  on  the  order  of  107/cm2. 


5 


The  complex  decay  dynamics  of  the  exciton  population  is  seen  in  the  FWMf 
response.  A  typical  line  shape  below  the  absorption  line  center  is  shown  in  Fig.  la.  The 
FWHM  corresponds  to  a  relaxation  time  of  60  ps,  which  is  too  slow  to  be  due  to  phonon 
scattering  of  the  delocalized  exciton  (typically  in  a  time  scale  of  several  ps16),  and  is  over 
an  older  of  magnitude  faster  than  the  exciton  recombination  time.  Furthermore,  the  FWHM 
is  independent  of  the  grating  period,  indicating  the  contribution  from  exciton  diffusion  is 
negligible  and  that  excitons  are  localized  in  this  spectral  region.  Hence,  the  data  suggest  that 
the  decay  is  due  to  spectral  diffusion  as  a  result  of  scattering  of  localized  excitons  from 
energy  E  to  E'  (IE-E’l>rh).  In  fact,  the  decay  rate  is  in  agreement  with  the  recent  calculation 
based  on  phonon  assisted  migration  of  the  localized  exciton5.  In  addition,  excitons  that  are 
scattered  from  E  to  E'  can  later  be  scattered  back  and  establish  a  quasi-equilibrium  exciton 
population  about  E  before  they  eventually  recombine.  This  quasi-equilibrium  population 
contributes  to  the  FWMf  response  with  a  decay  rate  characterized  by  the  exciton 
recombination13.  The  FWMf  response  shown  in  Fig.  la  indeed  shows  a  small  and  narrow 
feature  at  the  top  of  the  line  shape.  The  narrow  feature  becomes  more  pronounced  at  higher 
temperature  due  to  faster  exciton  migration  as  shown  in  Fig.  lb.  The  width  associated  with 
the  feature  corresponds  to  a  decay  time  of  1.2  ns,  consistent  with  the  exciton  recombination 
rate17. 

We  further  examine  the  relaxation  mechanism  of  the  localized  exciton  by  studying  the 
temperature  dependence  of  the  exciton  relaxation  rate.  Fig.  2a  shows  the  temperature 
dependence  of  the  exciton  relaxation  rate  obtained  at  0.6  meV  and  1.5  meV  below  the 
absorption  line  center  using  the  FWMf  response.  The  data  is  in  good  agreement  with  the 
theory  of  phonon  assisted  migration  discussed  above  with  a=1.6.  The  measurement 
indicates  that  the  dominant  contribution  to  relaxation  of  the  localized  exciton  is  phonon 
assisted  migration  up  to  a  temperature  of  15  K.  Note  that  earlier  measurements  using 
resonant  Rayleigh  scattering  have  reported  observations  of  activation  type  of  temperature 
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dependence  fen*  the  localized  exciton  at  temperatures  between  7  and  20  K,  indicating  that  in 
this  temperature  region  relaxation  for  the  localized  exciton  is  dominated  by  thermal 
activation  to  delocalized  states.  It  has  been  suggested  that  depending  on  details  of  quantum 
well  interface  roughness,  the  effective  activation  energy  can  be  much  higher  than  simply  the 
energy  difference  between  the  mobility  edge  and  the  localized  exciton*®.  Therefore  it  is 
possible  that  at  relatively  low  temperature  (e.g.  below  15  K),  the  thermal  activation  rate  is 
much  smaller  than  the  phonon  assisted  migration  rate. 

The  energy  dependence  of  the  decay  of  the  exciton  population  below  the  absorption 
line  center  is  shown  in  Fig.  2b.  For  exciton  energies  less  than  1.5  meV  below  the 
absorption  line  center,  the  decay  rate  depends  very  weakly  on  the  energy  .  The  rate 
increases  rapidly  when  the  exciton  energy  approaches  the  absorption  line  center,  suggesting 
a  transition  from  localized  to  delocalized  excitons3. 

If  we  imagine  that  localized  excitons  are  optically  excited  at  energy  E  and  then 
migrate  among  localized  sites  to  different  energies,  a  quasi-equilibrium  exciton  population 
over  a  broad  spectral  range  can  be  established  assuming  the  exciton  migration  rate  is  large 
compared  with  the  recombination  rate.  As  indicated  earlier,  indeed,  it  is  the  decay  of  this 
quasi-equilibrium  distribution  that  gives  rise  to  the  narrow  feature  in  Fig.  la,  b.  The 
spectral  distribution  of  the  population  can  be  directly  measured  in  the  FWMb  response  by 
scanning  0%  while  keeping  0)f  and  CDp  fixed  at  E.  Fig.  3  shows  a  FWMb  line  shape  where 
excitons  are  optically  excited  at  1.5  meV  below  the  absorption  line  center.  The  nonlinear 
response  is  corrected  for  sample  absorption,  and  is  proportional  to  the  quasi-equilibrium 
exciton  population  assuming  all  excitons  in  the  spectral  region  concerned  give  rise  to  the 
same  CW  nonlinear  response. 

To  quantitatively  describe  the  exciton  migration  process,  we  introduce  a  distribution 
kernel  f(E,  E'),  the  rate  for  a  exciton  to  migrate  from  energy  E  to  E'.  f(E,  E')  is  analogous 
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to  collision  kernels  which  are  used  to  describe  velocity  changing  collisions  in  atomic 
vapor19.  Decay  of  the  exciton  population  can  then  be  described  by  the  following  transport 
equation  p(E)  =  -(Y  +  rg)p(E)  +  / f(E',  E)p(E')dE',  where  p(E)  is  the  density  of  the 
exciton  at  E,  y  is  the  recombination  rate,  =  J  f  (E,  E'  )dE'  is  the  overall  rate  at  which 
excitons  migrate  out  of  states  with  energy  E.  In  principle,  the  transport  equation,  together 
with  the  equation  for  the  nonlinear  polarization11-20,  can  provide  a  description  for  the 
FWMf  and  FWMb  line  shape,  though  in  general,  f(E,  E')  is  unknown  and  has  to  be 
determined  by  experiments.  However,  by  assuming  f(E,  E’)  to  be  independent  of  the  initial 
state  energy,  analytical  solutions13  can  be  obtained  for  standard  distributions  and  show  that 
the  FWMf  response  has  an  additional  component  with  a  width  given  by  y  as  anticipated 
based  on  physical  arguments.  The  FWMb  line  shape  in  Fig.  3  can  be  described  by  a  simple 
model  which  neglects  migration  to  states  above  the  excitation  energy,  and  is  based  on  the 
understanding  of  the  solution  to  the  above  transport  equation.  The  model  assumes  a 
Gaussian  distribution  for  the  quasi-equilibrium  population  of  excitons  that  have  migrated  to 
states  below  the  excitation  energy.  The  result  is  plotted  as  the  solid  line  in  Fig.  3. 

Finally,  we  note  that  earlier  measurements  have  suggested  that  excitons  above  line 
center  are  weakly  delocalized3.  Indeed,  our  measurements  of  the  FWMf  response  show  a 
quadratic  dependence14  of  the  modulation  decay  rate  on  the  inverse  of  the  excitation  grating 
period  yielding  a  diffusion  coeffificient  on  the  order  of  4  cm2/s.  More  strikingly,  for  large 
grating  spacing  where  there  is  negligible  contribution  from  exciton  spatial  diffusion,  the 
FWMf  line  shape  is  completely  dominated  by  the  recombination  component  due  to  the  rapid 
exciton- phonon  scattering  (on  the  order  of  several  ps16).  A  typical  FWMf  line  shape  in  this 
case  is  shown  in  Fig.  lc.  The  FWHM  of  the  line  shape  corresponds  to  an  exciton 
recombination  time  of  1.2  ns.  Using  the  FWMb  response,  we  also  determined  the 
homogeneous  width  to  be  on  the  order  0.5  meV  (corresponding  to  a  dephasing  time  of  1.3 
ps )  at  a  temperature  of  5  K,  consistent  with  what  is  expected  for  the  delocalized  exciton16. 
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In  summary,  using  high  resolution  nonlinear  laser  spectroscopy,  we  have  shown  that 
lowest  energy  excitons  are  localized  and  migrate  in  energy  by  absorption  or  emission  of 
acoustic  phonons.  The  line  shapes  have  enabled  us  to  determine  the  migrating  rate  along 
with  the  temperature  dependence  and  the  resultant  quasi-equilibrium  energy  distribution. 

This  work  was  supported  by  the  Air  force  Office  of  Scientific  Research. 
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Figure  Captions: 


Fig.  1.  The  FWMf  responses,  (a)  at  1.5  meV  below  line  center  and  at  2.5  K.  (b)  at  0.6 
meV  below  line  center  and  at  10  K.  (c)  at  2  meV  above  line  center  and  at  10  K. 

Fig.  2.  (a)  The  temperature  dependence  of  the  exciton  decay  rate.  Circles  and  crosses  are 
data  obtained  at  1.5  meV  and  at  0.6  meV  below  line  center  respectively.  Dashed  lines  are  fit 
to  the  theory  of  phonon  assisted  migration,  (b)  The  energy  dependence  of  the  exciton 
decay  rate  at  10  K.  The  dash  line  is  only  a  guide  to  the  eye. 

Fig.  3.  The  FWMb  response  with  the  frequencies  of  the  forward  and  probe  beams  at  1.5 
meV  below  line  center. 
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